CN106097450A - A kind of contour lines creation method and device - Google Patents
A kind of contour lines creation method and device Download PDFInfo
- Publication number
- CN106097450A CN106097450A CN201610454545.9A CN201610454545A CN106097450A CN 106097450 A CN106097450 A CN 106097450A CN 201610454545 A CN201610454545 A CN 201610454545A CN 106097450 A CN106097450 A CN 106097450A
- Authority
- CN
- China
- Prior art keywords
- equal pitch
- contour
- pitch contour
- elevation
- line
- 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/05—Geographic models
-
- 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
-
- 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/04—Architectural design, interior design
Abstract
The present invention relates to a kind of contour lines creation method and device, belong to GIS-Geographic Information System field.The present invention has initially set up digital elevation model based on Delaunay triangulation network, improves the form of net form.Then according to various form lines, geosutures and cultural symbol, net form is modified, thus ensure that DTM farthest reduces landform original appearance according to measured data.Follow the trail of equal pitch contour further according to this accurate DTM and then ensure that isocontour precision.Finally, the Curvilinear polygon border of cultural symbol is turned off by equal pitch contour, preferably solve equal pitch contour and meet atural object, particularly overlapping atural object and the disconnection process problem of large area atural object.The DTM that the present invention sets up can be used not only for following the trail of equal pitch contour, is additionally operable to the cubic metre of earth and stone and calculates and the field such as section process and engineering design.So, the present invention also establishes elevation model at slope and road etc. without the region drawn a contour line.
Description
Technical field
The present invention relates to a kind of contour lines creation method and device, belong to GIS-Geographic Information System field.
Background technology
In topography, represent that a kind of important means of morphologic prominence is equal pitch contour, in conventional panel mapping, equal pitch contour
By artificial hand drawing in the wild.In mapping automatically, computer isocontour algorithm of automatically sketching is considerably complicated
, also it is one of the difficult point of automatic cartography system.It first has to set up digital elevation model, i.e. field acquisition arrive mutually the most not
Relevant discrete point carries out networking, i.e. sets up the contact between corresponding discrete point, then directly carries out in network model
Interpolation and equal pitch contour of sketching out (also referred to as contour).
Equal pitch contour in map, is different from other general isogram, and its precision prescribed is higher, is setting up simultaneously
During network model, various form line and geosutures must be considered, such as ridge, mountain valley and steep bank etc..It addition, be accustomed to according to map
Requirement, equal pitch contour also needs to when running into atural object, annotation etc. disconnect, these, bring relatively for computer automatic drafting equal pitch contour
Big difficulty.
Summary of the invention
The present invention mainly solves the above-mentioned technical problem existing for prior art, it is provided that a kind of contour lines creation side
Method and device.The method and device have used optimisation technique as advanced in plug-in type constructor algorithm, dynamic increment index etc. in a large number,
Improve isocontour formation efficiency, and preferably solve equal pitch contour and meet atural object, particularly overlapping atural object and in large area
The disconnection process problem of thing.
The above-mentioned technical problem of the present invention is mainly addressed by following technical proposals:
A kind of contour lines creation method, including:
Step 1, generates equal pitch contour, confirms scale and contour interval, reads in the dot file of pending figure, according in figure
Having elevational point to carry out automatic network forming, in net type, selecting blueness be common triangle selvage, and redness and green are limited deckle, wherein,
Limiting limit and be restricted limit, it can not swap limit exactly with the difference on common limit, and it is necessary in order to carry out for limiting limit
The needs that disconnect of atural object and the limit specified；Isocontour tracking is just carried out after network forming；To the triangulation network before following the trail of equal pitch contour
Carry out suitable finishing, specifically include:
Step 101, when carrying out triangle network forming, if certain angle of triangle >=135 degree, then easily cause this
Place's equal pitch contour distortion, thus affect data precision, for avoiding this kind of situation, at the tetragon that adjacent two obtuse triangles are constituted
In, this tetragon diagonal is swapped, thus in the case of not deleting triangle, the triangle at bigger obtuse angle is become
Acute triangle；
Step 102, will be not required to network forming or the underproof limit of net type is left out；
Step 103, adds net type limit and is maybe again recovered on the limit by mistake deleted；
Step 104, additions and deletions border, additions and deletions border has delimitation equal pitch contour forbidden zone and generates the double function in equal pitch contour region,
Delimiting equal pitch contour forbidden zone is after overall network forming, marks additions and deletions border, is automatically deleted by the triangulation network in the range of this, superposition etc.
Equal pitch contour in the range of being then somebody's turn to do during high line disconnects automatically；Delimiting equal pitch contour forbidden zone is before network forming, first delimits the border of additions and deletions,
Just network forming in the region then only delimited；
Step 2, isocontour superposition and graphics process, open and need the isocontour figure of superposition, self-defined contour number
Spacing and decimal digits, scatterplot annotation spacing and decimal digits thereof, choose whether equal pitch contour is carried out polishing process；Wherein,
If choosing this, then the equal pitch contour generated smooths attractive in appearance, but degree of accuracy may decrease；If not choosing, then equal pitch contour is relatively
Truly, but the most smooth；Then, high line can be added on figure automatically, topography now is carried out following process:
Step 201, is rule to the direction of end elevation by the direction of elevation on figure, adopts what line intersected with equal pitch contour at screen
Local generation annotation, try again equal pitch contour disconnection process, then at annotation, equal pitch contour disconnects automatically；
Step 202, as required to discrete point annotation in addition in self-defining scope, specific practice is: first draw certain
One region, carries out annotation to the scatterplot in the range of this, if not drawing the most any scope, then carries out annotation in the range of full figure.
Step 203, uses equal pitch contour fairing mode to adjust equal pitch contour smoothness, including: tensioned spline interpolation, H sample
Bar interpolation and B-spline matching and merge closely three kinds of modes；
Step 3, uses Triangular meshes method to carry out Earthwork calculation, including step chosen below:
Step 301, uses horizontal reference surface model to carry out Earthwork calculation, wherein, horizontal reference surface model
The relative plane of reference referring to fill out excavation is the situation of horizontal plane, first reads in the eb literary composition of topography participating in calculating
Part, then carries out region network forming, then input reference face elevation, and wherein in this dialog box, the elevation of acquiescence is the highest of current figure
Elevation and the meansigma methods of lowest elevation, determine the region participating in Earthwork Calculation the most again, then calculate volume；
Step 302, uses Tilt Datum Plane surface model to carry out Earthwork calculation, and wherein, Tilt Datum Plane surface model refers to fill out excavation
The relative plane of reference be the situation of plane tilted, its relative level has an inclination angle；Tilt Datum Plane surface model parameter whole
Body meaning is: inclined-plane is a little determined by outside datum line and line, and the elevation of datum line is the elevation of the horizontal plane that datum line is passed through,
And the angle of some inclination angle i.e. inclined plane and horizontal plane of (i.e. orientation point) relative level outside line；Wherein, Tilt Datum Plane face mould
Formula needs to select Tilt Datum Plane face, then inclination angle, input reference face while the elevation of input reference face in parameter is arranged；
Step 302, uses any base level model to carry out Earthwork calculation, and wherein, any base level model refers to carry out soil
The relative datum face that side's amount calculates is the face of arbitrary shape；Including: structure bilayer model, first reads in the first model and i.e. fills out and dig
After object module, then network forming, reselection constructs the second model, reads in the second model, i.e. master mould, then network forming, then selects
Select "current" model, return to the first model, delimit works scope, next can carry out Earthwork Calculation.
Optimize, above-mentioned contour lines creation method, use equal pitch contour fairing mode to adjust equal pitch contour smoothness, bag
Include: tensioned spline interpolation, B-spline matching, H spline interpolation.
Optimize, above-mentioned contour lines creation method, also include step 4, read in the file of the topography participating in calculating, structure
Build the triangulation network, section part line need to increased, generating section according to line.
Optimize, above-mentioned contour lines creation method, also include step 5, after reading in equal pitch contour process file network forming, screen
Display is divided into two parts: areal model edit segment and 3-D effect show part.
A kind of contour lines creation device, including:
Contour lines creation module, confirms scale and contour interval, reads in the dot file of pending figure, according to having in figure
Elevational point carries out automatic network forming, and in net type, selecting blueness be common triangle selvage, and red and green is limited deckle, wherein, limits
Deckle is restricted limit, and it can not swap limit exactly with the difference on common limit, and it is necessary in order to carry out for limiting limit
Needs that atural object disconnects and the limit specified；Isocontour tracking is just carried out after network forming；Before following the trail of equal pitch contour, the triangulation network is entered
The finishing that row is suitable, specifically includes:
Opposite side crosspoint, when carrying out triangle network forming, if certain angle of triangle >=135 degree, the most easily make
Become equal pitch contour distortion at this, thus affect data precision, for avoiding this kind of situation, on four limits that adjacent two obtuse triangles are constituted
In shape, this tetragon diagonal is swapped, thus by the triangularity at bigger obtuse angle in the case of not deleting triangle
Acutangulate triangle；
Straight line deletes net unit, will be not required to network forming or the underproof limit of net type is left out；
Unit is inserted on limit, adds net type limit and is maybe again recovered on the limit by mistake deleted；
Border additions and deletions unit, additions and deletions border has delimitation equal pitch contour forbidden zone and generates the double function in equal pitch contour region, delimiting
Equal pitch contour forbidden zone is after overall network forming, marks additions and deletions border, is automatically deleted by the triangulation network in the range of this, superposition equal pitch contour
Time then should in the range of equal pitch contour automatically disconnect；Delimiting equal pitch contour forbidden zone is before network forming, first delimits the border of additions and deletions, then only
There is just network forming in the region of delimitation；
Superposition and pattern process module, open and need the isocontour figure of superposition, self-defined contour number spacing and little
Numerical digit number, scatterplot annotation spacing and decimal digits thereof, choose whether equal pitch contour is carried out polishing process；Wherein, if choosing this,
The equal pitch contour then generated smooths attractive in appearance, but degree of accuracy may decrease；If not choosing, then equal pitch contour is relatively true, but differs
Fixed the most smooth；Then, high line can be added on figure automatically, topography now is carried out following process:
Locally equal pitch contour mark unit, on figure by the direction of elevation to the direction of end elevation rule, screen adopt line with
The local generation annotation that equal pitch contour intersects, try again equal pitch contour disconnection process, then at annotation, equal pitch contour disconnects automatically；
Local Discrete point annotation unit, as required to discrete point annotation, specific practice in addition in self-defining scope
Being: first draw a certain region, the scatterplot in the range of this is carried out annotation, if not drawing the most any scope, then entering in the range of full figure
Row annotation.
Fairing mode selects unit, uses equal pitch contour fairing mode to adjust equal pitch contour smoothness, including: tensioned batten
Interpolation, H spline interpolation and B-spline matching and merging closely three kinds of modes；
Earthwork calculation module, uses Triangular meshes method to carry out Earthwork calculation, including step chosen below:
Horizontal pattern unit, uses horizontal reference surface model to carry out Earthwork calculation, and wherein, horizontal reference surface model refers to fill out
The relative plane of reference of excavation is the situation of horizontal plane, first reads in the eb file of the topography participating in calculating, then carries out region
Network forming, then input reference face elevation, wherein in this dialog box, the elevation of acquiescence is elevation and the lowest elevation of current figure
Meansigma methods, the most again determine participate in Earthwork Calculation region, then calculate volume；
Tilt mode unit, uses Tilt Datum Plane surface model to carry out Earthwork calculation, and wherein, Tilt Datum Plane surface model refers to fill out
The relative plane of reference of excavation is the situation of the plane tilted, and its relative level has an inclination angle；Tilt Datum Plane surface model parameter
Whole meaning be: inclined-plane is a little determined by outside datum line and line, and the elevation of datum line is the horizontal plane that datum line is passed through
Elevation, and the angle of some inclination angle i.e. inclined plane and horizontal plane of (i.e. orientation point) relative level outside line；Wherein, Tilt Datum Plane
Surface model needs to select Tilt Datum Plane face in parameter is arranged, and then while the elevation of input reference face, incline in input reference face
Angle；
Arbitrary patterns unit, uses any base level model to carry out Earthwork calculation, and wherein, any base level model is fingering
The relative datum face of row Earthwork Calculation is the face of arbitrary shape；Including: structure bilayer model, first reads in the first model i.e.
Filling out the object module after digging, then network forming, reselection constructs the second model, reads in the second model, i.e. master mould, then network forming,
Reselection "current" model, returns to the first model, delimit works scope, next can carry out Earthwork Calculation.
Optimize, above-mentioned contour lines creation device, use equal pitch contour fairing mode to adjust equal pitch contour smoothness, bag
Include: tensioned spline interpolation, B-spline matching, H spline interpolation.
Optimize, above-mentioned contour lines creation device, also include section generation module, read in the topography of participation calculating
File, constructs the triangulation network, need to increase section part line, generate section according to line.
Optimize, above-mentioned contour lines creation device, also include display module, after reading in equal pitch contour process file network forming,
Screen shows and is divided into two parts: areal model edit segment and 3-D effect display part.
Therefore, present invention have the advantage that the method and device have used such as plug-in type constructor algorithm in a large number, dynamically increased
The advanced optimisation techniques such as amount index, improve isocontour formation efficiency, and preferably solve equal pitch contour chance atural object, special
It not overlapping atural object and the disconnection process problem of large area atural object.
Accompanying drawing explanation
Accompanying drawing 1 is the schematic diagram before exchange opposite side；
Accompanying drawing 2 is the schematic diagram after exchange opposite side；
Accompanying drawing 3 is the schematic diagram before additions and deletions border；
Accompanying drawing 4 is the schematic diagram behind additions and deletions border.
Detailed description of the invention
Below by embodiment, and combine accompanying drawing, technical scheme is described in further detail.
In Fig. 1, the opposite side that 1-is exchanged.
Embodiment:
A kind of contour lines creation method, including:
Step 1, generates equal pitch contour, confirms scale and contour interval, reads in the dot file of pending figure, according in figure
Having elevational point to carry out automatic network forming, in net type, selecting blueness be common triangle selvage, and redness and green are limited deckle, wherein,
Limiting limit and be restricted limit, it can not swap limit exactly with the difference on common limit, and it is necessary in order to carry out for limiting limit
The needs that disconnect of atural object and the limit specified；Isocontour tracking is just carried out after network forming；To the triangulation network before following the trail of equal pitch contour
Carry out suitable finishing, specifically include:
Step 101, when carrying out triangle network forming, if certain angle of triangle >=135 degree, then easily cause this
Place's equal pitch contour distortion, thus affect data precision, for avoiding this kind of situation, at the tetragon that adjacent two obtuse triangles are constituted
In, this tetragon diagonal is swapped, thus in the case of not deleting triangle, the triangle at bigger obtuse angle is become
Acute triangle；
Step 102, will be not required to network forming or the underproof limit of net type is left out；
Step 103, adds net type limit and is maybe again recovered on the limit by mistake deleted；
Step 104, additions and deletions border, additions and deletions border has delimitation equal pitch contour forbidden zone and generates the double function in equal pitch contour region,
Delimiting equal pitch contour forbidden zone is after overall network forming, marks additions and deletions border, is automatically deleted by the triangulation network in the range of this, superposition etc.
Equal pitch contour in the range of being then somebody's turn to do during high line disconnects automatically；Delimiting equal pitch contour forbidden zone is before network forming, first delimits the border of additions and deletions,
Just network forming in the region then only delimited；
Step 2, isocontour superposition and graphics process, open and need the isocontour figure of superposition, self-defined contour number
Spacing and decimal digits, scatterplot annotation spacing and decimal digits thereof, choose whether equal pitch contour is carried out polishing process；Wherein,
If choosing this, then the equal pitch contour generated smooths attractive in appearance, but degree of accuracy may decrease；If not choosing, then equal pitch contour is relatively
Truly, but the most smooth；Then, high line can be added on figure automatically, topography now is carried out following process:
Step 201, is rule to the direction of end elevation by the direction of elevation on figure, adopts what line intersected with equal pitch contour at screen
Local generation annotation, try again equal pitch contour disconnection process, then at annotation, equal pitch contour disconnects automatically；
Step 202, as required to discrete point annotation in addition in self-defining scope, specific practice is: first draw certain
One region, carries out annotation to the scatterplot in the range of this, if not drawing the most any scope, then carries out annotation in the range of full figure.
Step 203, uses equal pitch contour fairing mode to adjust equal pitch contour smoothness, including: tensioned spline interpolation, H sample
Bar interpolation and B-spline matching and merge closely three kinds of modes；
Step 3, uses Triangular meshes method to carry out Earthwork calculation, including step chosen below:
Step 301, uses horizontal reference surface model to carry out Earthwork calculation, and wherein, horizontal reference surface model refers to fill out excavation
The relative plane of reference be the situation of horizontal plane, first read in the eb file of topography participating in calculating, then carry out region network forming,
Input reference face elevation again, wherein in this dialog box, the elevation of acquiescence is the average of elevation and the lowest elevation of current figure
Value, determines the region participating in Earthwork Calculation the most again, then calculates volume；
Step 302, uses Tilt Datum Plane surface model to carry out Earthwork calculation, and wherein, Tilt Datum Plane surface model refers to fill out excavation
The relative plane of reference be the situation of plane tilted, its relative level has an inclination angle；Tilt Datum Plane surface model parameter whole
Body meaning is: inclined-plane is a little determined by outside datum line and line, and the elevation of datum line is the elevation of the horizontal plane that datum line is passed through,
And the angle of some inclination angle i.e. inclined plane and horizontal plane of (i.e. orientation point) relative level outside line；Wherein, Tilt Datum Plane face mould
Formula needs to select Tilt Datum Plane face, then inclination angle, input reference face while the elevation of input reference face in parameter is arranged；
Step 302, uses any base level model to carry out Earthwork calculation, and wherein, any base level model refers to carry out soil
The relative datum face that side's amount calculates is the face of arbitrary shape；Including: structure bilayer model, first reads in the first model and i.e. fills out and dig
After object module, then network forming, reselection constructs the second model, reads in the second model, i.e. master mould, then network forming, then selects
Select "current" model, return to the first model, delimit works scope, next can carry out Earthwork Calculation；
Step 4, reads in the file of the topography participating in calculating, constructs the triangulation network, need to increase section part line, according to drawing
Line generates section；
Step 5, reads in equal pitch contour and processes after file network forming, and screen show and is divided into two parts: areal model edit segment with
3-D effect display part.
Wherein, equal pitch contour fairing mode is used to adjust equal pitch contour smoothness, including: tensioned spline interpolation, B-spline are intended
Conjunction, H spline interpolation.
Tensioned spline interpolation, tension-adjustable coefficient, its coefficient value is from 1 to 10, and coefficient is the least, and equal pitch contour is the most smooth,
Coefficient is the biggest, and equal pitch contour is the sharpest.When the difference that the equal pitch contour coefficient of tension and polishing process is by polishing process, equal pitch contour passes through
Summit change, adjust the equal pitch contour coefficient of tension and the most only adjust isocontour smoothness, the summit that equal pitch contour passes through is not
Become；
B-spline matching, can make equal pitch contour relatively smooth, reaches degree the most attractive in appearance, but owing to it is the most contour
The approximate fits of line, so its error is relatively large, can select according to practical situation oneself this user, as wanted pattern precision
Ask the highest, B-spline matching can be carried out and polishing processes, if but higher to required precision, it is proposed that only adjust tensile splines interpolation i.e.
Can.
H spline interpolation, is that equal pitch contour is carried out segment processing, compared with B-spline matching, its refresh rate faster, precision
Also higher than B-spline matching.
After using said method, the present invention has initially set up digital elevation model based on Delaunay triangulation network, improves
The form of net form.Then according to various form lines, geosutures and cultural symbol, net form is modified, thus ensure that DTM
Landform original appearance is farthest reduced according to measured data.Follow the trail of equal pitch contour further according to this accurate DTM then to ensure that
The precision of high line.Finally, the Curvilinear polygon border of cultural symbol is turned off by equal pitch contour, preferably solves equal pitch contour
Meet atural object, particularly overlapping atural object and the disconnection process problem of large area atural object.The DTM that the present invention sets up can be used not only for following the trail of
Equal pitch contour, is additionally operable to the cubic metre of earth and stone and calculates and the field such as section process and engineering design.So, the present invention is in the nothing such as slope and road
The region that need to draw a contour line also establishes elevation model.
In above-mentioned processing procedure, the present invention has used in a large number such as advanced persons such as plug-in type constructor algorithm, dynamic increment indexes
Optimisation technique, improve contour lines creation efficiency.
Specific embodiment described herein is only to present invention spirit explanation for example.Technology neck belonging to the present invention
Described specific embodiment can be made various amendment or supplements or use similar mode to replace by the technical staff in territory
Generation, but without departing from the spirit of the present invention or surmount scope defined in appended claims.
Claims (8)
1. a contour lines creation method, it is characterised in that including:
Step 1, generates equal pitch contour, confirms scale and contour interval, reads in the dot file of pending figure, according to the existing height in figure
Journey point carries out automatic network forming, and in net type, selecting blueness be common triangle selvage, and red and green is limited deckle, wherein, restriction
While be restricted limit, its difference with common limit can not swap limit exactly, and it is necessary in order to carry out for limiting limit
Needs that thing disconnects and the limit specified；Isocontour tracking is just carried out after network forming；Before following the trail of equal pitch contour, the triangulation network is carried out
Suitable finishing, specifically includes:
Step 101, when carrying out triangle network forming, if certain angle of triangle >=135 degree, then easily cause at this etc.
High line distortion, thus affect data precision, for avoiding this kind of situation, in the tetragon that adjacent two obtuse triangles are constituted, will
This tetragon diagonal swaps, thus in the case of not deleting triangle, the triangularity at bigger obtuse angle is acutangulated three
Dihedral；
Step 102, will be not required to network forming or the underproof limit of net type is left out；
Step 103, adds net type limit and is maybe again recovered on the limit by mistake deleted；
Step 104, additions and deletions border, additions and deletions border has delimitation equal pitch contour forbidden zone and generates the double function in equal pitch contour region, delimiting
Equal pitch contour forbidden zone is after overall network forming, marks additions and deletions border, is automatically deleted by the triangulation network in the range of this, superposition equal pitch contour
Time then should in the range of equal pitch contour automatically disconnect；Delimiting equal pitch contour forbidden zone is before network forming, first delimits the border of additions and deletions, then only
There is just network forming in the region of delimitation；
Step 2, isocontour superposition and graphics process, open and need the isocontour figure of superposition, self-defined contour number spacing
And decimal digits, scatterplot annotation spacing and decimal digits thereof, choose whether equal pitch contour is carried out polishing process；Wherein, if selecting
In this, then the equal pitch contour generated smooths attractive in appearance, but degree of accuracy may decrease；If not choosing, then equal pitch contour is relatively true,
But it is the most smooth；Then, high line can be added on figure automatically, topography now is carried out following process:
Step 201, is rule to the direction of end elevation by the direction of elevation on figure, adopts, at screen, the place that line intersects with equal pitch contour
Generating annotation, try again equal pitch contour disconnection process, then at annotation, equal pitch contour disconnects automatically；
Step 202, as required to discrete point annotation in addition in self-defining scope, specific practice is: first draw a certain district
Territory, carries out annotation to the scatterplot in the range of this, if not drawing the most any scope, then carries out annotation in the range of full figure.
Step 203, uses equal pitch contour fairing mode to adjust equal pitch contour smoothness, including: tensioned spline interpolation, H batten are inserted
Value and B-spline matching and merging closely three kinds of modes；
Step 3, uses Triangular meshes method to carry out Earthwork calculation, including step chosen below:
Step 301, uses horizontal reference surface model to carry out Earthwork calculation, and wherein, horizontal reference surface model refers to fill out the phase of excavation
It is the situation of horizontal plane to the plane of reference, first reads in the eb file of the topography participating in calculating, then carry out region network forming, more defeated
Entering datum elevation, wherein in dialog box, the elevation of acquiescence is the elevation meansigma methods with lowest elevation of current figure, connects
Get off and determine the region participating in Earthwork Calculation again, then calculate volume；
Step 302, uses Tilt Datum Plane surface model to carry out Earthwork calculation, and wherein, Tilt Datum Plane surface model refers to fill out the phase of excavation
Being the situation of the plane tilted to the plane of reference, its relative level has an inclination angle；The overall meaning of Tilt Datum Plane surface model parameter
Justice is: inclined-plane is a little determined by outside datum line and line, and the elevation of datum line is the elevation of the horizontal plane that datum line is passed through, and line
The inclination angle i.e. inclined plane of outer a bit (i.e. orientation point) relative level and the angle of horizontal plane；Wherein, Tilt Datum Plane surface model needs
Tilt Datum Plane face, then inclination angle, input reference face while the elevation of input reference face is selected in parameter is arranged；
Step 302, uses any base level model to carry out Earthwork calculation, and wherein, any base level model refers to carry out earth volume
The relative datum face calculated is the face of arbitrary shape；Including structure bilayer model, first read in after the first model i.e. fills out and dig
Object module, then network forming, reselection constructs the second model, reads in the second model, i.e. master mould, then network forming, and reselection is worked as
Front model, returns to the first model, delimit works scope, next can carry out Earthwork Calculation.
Contour lines creation method the most according to claim 1, it is characterised in that use equal pitch contour fairing mode to adjust contour
Line smoothness, including: tensioned spline interpolation, B-spline matching, H spline interpolation.
Contour lines creation method the most according to claim 1, it is characterised in that also include step 4, reads in and participates in calculating
The file of topography, constructs the triangulation network, need to increase section part line, generate section according to line.
Contour lines creation method the most according to claim 1, it is characterised in that also include step 5, reads in equal pitch contour and processes
After file network forming, screen shows and is divided into two parts: areal model edit segment and 3-D effect display part.
5. a contour lines creation device, it is characterised in that including:
Contour lines creation module, confirms scale and contour interval, reads in the dot file of pending figure, according to the existing elevation in figure
Point carries out automatic network forming, and in net type, selecting blueness be common triangle selvage, and red and green is limited deckle, wherein, restriction limit
Being restricted limit, it can not swap limit exactly with the difference on common limit, and limiting limit is to carry out necessary atural object
The needs that disconnect and the limit specified；Isocontour tracking is just carried out after network forming；Before following the trail of equal pitch contour, the triangulation network is fitted
When finishing, specifically include:
Opposite side crosspoint, when carrying out triangle network forming, if certain angle of triangle >=135 degree, then easily cause this
Place's equal pitch contour distortion, thus affect data precision, for avoiding this kind of situation, at the tetragon that adjacent two obtuse triangles are constituted
In, this tetragon diagonal is swapped, thus in the case of not deleting triangle, the triangle at bigger obtuse angle is become
Acute triangle；
Straight line deletes net unit, will be not required to network forming or the underproof limit of net type is left out；
Unit is inserted on limit, adds net type limit and is maybe again recovered on the limit by mistake deleted；
Border additions and deletions unit, additions and deletions border has delimitation equal pitch contour forbidden zone and generates the double function in equal pitch contour region, delimiting contour
Line forbidden zone is after overall network forming, marks additions and deletions border, is automatically deleted by the triangulation network in the range of this, during superposition equal pitch contour then
Equal pitch contour in the range of Gai disconnects automatically；Delimiting equal pitch contour forbidden zone is before network forming, first delimit the border of additions and deletions, then only draws
Just network forming in fixed region；
Superposition and pattern process module, open and need the isocontour figure of superposition, self-defined contour number spacing and decimal place thereof
Number, scatterplot annotation spacing and decimal digits thereof, choose whether equal pitch contour is carried out polishing process；Wherein, if choosing this, then give birth to
The equal pitch contour become smooths attractive in appearance, but degree of accuracy may decrease；If not choosing, then equal pitch contour is relatively true, but the most very
Smooth；Then, high line can be added on figure automatically, topography now is carried out following process:
Locally equal pitch contour mark unit, is rule to the direction of end elevation by the direction of elevation on figure, adopts line with contour at screen
The local generation annotation that line intersects, try again equal pitch contour disconnection process, then at annotation, equal pitch contour disconnects automatically；
Local Discrete point annotation unit, as required to discrete point in addition annotation in self-defining scope, specific practice is: first
Drawing a certain region, the scatterplot in the range of this is carried out annotation, if not drawing the most any scope, then noting in the range of full figure
Note.
Fairing mode selects unit, uses equal pitch contour fairing mode to adjust equal pitch contour smoothness, including: tensioned batten is inserted
Value, H spline interpolation and B-spline matching and merging closely three kinds of modes；
Earthwork calculation module, uses Triangular meshes method to carry out Earthwork calculation, including step chosen below:
Horizontal pattern unit, uses horizontal reference surface model to carry out Earthwork calculation, and wherein, horizontal reference surface model refers to fill out excavation
The relative plane of reference be the situation of horizontal plane, first read in the eb file of topography participating in calculating, then carry out region network forming,
Input reference face elevation again, wherein in this dialog box, the elevation of acquiescence is the average of elevation and the lowest elevation of current figure
Value, determines the region participating in Earthwork Calculation the most again, then calculates volume；
Tilt mode unit, uses Tilt Datum Plane surface model to carry out Earthwork calculation, and wherein, Tilt Datum Plane surface model refers to fill out excavation
The relative plane of reference be the situation of plane tilted, its relative level has an inclination angle；Tilt Datum Plane surface model parameter whole
Body meaning is: inclined-plane is a little determined by outside datum line and line, and the elevation of datum line is the elevation of the horizontal plane that datum line is passed through,
And the angle of some inclination angle i.e. inclined plane and horizontal plane of (i.e. orientation point) relative level outside line；Wherein, Tilt Datum Plane face mould
Formula needs to select Tilt Datum Plane face, then inclination angle, input reference face while the elevation of input reference face in parameter is arranged；
Arbitrary patterns unit, uses any base level model to carry out Earthwork calculation, and wherein, any base level model refers to carry out soil
The relative datum face that side's amount calculates is the face of arbitrary shape；Including: structure bilayer model, first reads in the first model and i.e. fills out and dig
After object module, then network forming, reselection constructs the second model, reads in the second model, i.e. master mould, then network forming, then selects
Select "current" model, return to the first model, delimit works scope, next can carry out Earthwork Calculation.
Contour lines creation device the most according to claim 5, it is characterised in that use equal pitch contour fairing mode to adjust contour
Line smoothness, including: tensioned spline interpolation, B-spline matching, H spline interpolation.
Contour lines creation device the most according to claim 5, it is characterised in that also include section generation module, reads in ginseng
With the file of the topography calculated, construct the triangulation network, section part line need to increased, generating section according to line.
Contour lines creation device the most according to claim 5, it is characterised in that also include display module, reads in equal pitch contour
Processing after file network forming, screen show and is divided into two parts: areal model edit segment and 3-D effect show part.
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811364179.3A CN109523631B (en) | 2016-06-21 | 2016-06-21 | Contour line generation method and device based on B spline fitting |
| CN201811364177.4A CN109523630B (en) | 2016-06-21 | 2016-06-21 | Contour line generation method and device based on tension spline interpolation |
| CN201610454545.9A CN106097450B (en) | 2016-06-21 | 2016-06-21 | A kind of contour lines creation method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610454545.9A CN106097450B (en) | 2016-06-21 | 2016-06-21 | A kind of contour lines creation method and device |
Related Child Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811364177.4A Division CN109523630B (en) | 2016-06-21 | 2016-06-21 | Contour line generation method and device based on tension spline interpolation |
| CN201811364179.3A Division CN109523631B (en) | 2016-06-21 | 2016-06-21 | Contour line generation method and device based on B spline fitting |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106097450A true CN106097450A (en) | 2016-11-09 |
| CN106097450B CN106097450B (en) | 2019-01-22 |
Family
ID=57238810
Family Applications (3)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811364177.4A Active CN109523630B (en) | 2016-06-21 | 2016-06-21 | Contour line generation method and device based on tension spline interpolation |
| CN201811364179.3A Active CN109523631B (en) | 2016-06-21 | 2016-06-21 | Contour line generation method and device based on B spline fitting |
| CN201610454545.9A Active CN106097450B (en) | 2016-06-21 | 2016-06-21 | A kind of contour lines creation method and device |
Family Applications Before (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811364177.4A Active CN109523630B (en) | 2016-06-21 | 2016-06-21 | Contour line generation method and device based on tension spline interpolation |
| CN201811364179.3A Active CN109523631B (en) | 2016-06-21 | 2016-06-21 | Contour line generation method and device based on B spline fitting |
Country Status (1)
| Country | Link |
|---|---|
| CN (3) | CN109523630B (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106557573A (en) * | 2016-11-24 | 2017-04-05 | 中铁第五勘察设计院集团有限公司 | The definition of TIN data memory format and data processing method based on directory index structure |
| CN106592357A (en) * | 2016-12-14 | 2017-04-26 | 中国人民解放军空军工程大学 | Airfield pavement curve transitional surface elevation designing method |
| CN106909740A (en) * | 2017-03-02 | 2017-06-30 | 中国石油天然气集团公司 | The computational methods and computing device of a kind of Oil ＆ Gas Storage yard total figure volume of earthwork |
| CN107067446A (en) * | 2017-01-19 | 2017-08-18 | 中国测绘科学研究院 | A kind of contour proximity relations judges and direction automatic adjusting method |
| CN107145669A (en) * | 2017-05-08 | 2017-09-08 | 攀钢集团矿业有限公司 | A kind of ore-rock amount computational methods and its device |
| CN109255008A (en) * | 2018-09-28 | 2019-01-22 | 贵州省水利水电勘测设计研究院 | A kind of data processing method, system and electronic equipment |
| CN112419495A (en) * | 2020-10-26 | 2021-02-26 | 天津大学 | Elevation point automatic extraction method based on multi-scale DEM space model |
| CN114046776A (en) * | 2021-09-22 | 2022-02-15 | 北京洛斯达科技发展有限公司 | Power transmission engineering water and soil conservation measure implementation checking system |
| US11618477B2 (en) * | 2020-07-16 | 2023-04-04 | Waymo Llc | Service area maps for autonomous vehicles |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110414095B (en) * | 2019-07-11 | 2023-06-27 | 上海交通大学 | Data preprocessing method in fluid-solid load spline interpolation conversion |
| CN110489511B (en) * | 2019-08-23 | 2021-07-30 | 西南大学 | Contour line edge-contact elevation error correction method and system, electronic device and medium |
| CN111724402B (en) * | 2020-06-18 | 2021-07-20 | 北京小白世纪网络科技有限公司 | Medical image labeling method, system and device |
| CN114329575B (en) * | 2022-03-15 | 2022-06-10 | 武汉峰岭科技有限公司 | Contour line encryption method, system, medium, and device based on interpolated ridge line |
| CN115329622A (en) * | 2022-07-16 | 2022-11-11 | 南京尚水数能科技有限公司 | Method for generating isoline based on channel terrain sounding point construction grid |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0562449A2 (en) * | 1992-03-19 | 1993-09-29 | Matsushita Electric Industrial Co., Ltd. | Object shape expression apparatus |
| WO1996013812A1 (en) * | 1994-10-29 | 1996-05-09 | Masaharu Ishii | Apparatus and method for terrain model reproduction |
| JP2005234787A (en) * | 2004-02-18 | 2005-09-02 | Yoshiaki Goto | Visualization technique using cg |
| KR100898617B1 (en) * | 2008-11-20 | 2009-05-27 | 주식회사 범아엔지니어링 | Construction method for digital elevation model of area coexisting the ground and water through verification of tin data of lidar and mbes measure value |
| CN102314711A (en) * | 2010-07-01 | 2012-01-11 | 中国地质科学院矿产资源研究所 | Three-dimensional visualization method and device for mineral resource evaluation information |
| CN102682479A (en) * | 2012-04-13 | 2012-09-19 | 国家基础地理信息中心 | Method for generating three-dimensional terrain feature points on irregular triangulation network |
| CN103226845A (en) * | 2013-03-18 | 2013-07-31 | 北京市测绘设计研究院 | TIN-based earth surface refinement processing method and system |
| CN103530901A (en) * | 2013-10-17 | 2014-01-22 | 中国测绘科学研究院 | Method and system for visualizing dynamic water area by matching with real three-dimensional environment |
| CN103871102A (en) * | 2014-03-28 | 2014-06-18 | 南京大学 | Road three-dimensional fine modeling method based on elevation points and road outline face |
| CN104331389A (en) * | 2014-11-21 | 2015-02-04 | 南京南瑞集团公司 | Eight-point method-based contour-line tracing algorithm |
| CN105303590A (en) * | 2015-11-25 | 2016-02-03 | 中国地质大学（武汉） | Method for cutting block model of three-dimensional ore body with attribute to generate two-dimensional profile map |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101750037A (en) * | 2008-11-28 | 2010-06-23 | 中国二十冶建设有限公司 | Accurate measuring method of large-scale irregular earthwork excavated volume |
| KR100916474B1 (en) * | 2009-04-28 | 2009-09-08 | 최준수 | Methods of partitioning a region represented by contours into smaller polygonal zones and calculating data for digital elevation model and data for constructing geographical surfaces of the region, and a recording medium in which programs for implementing the methods are recorded |
| CN102129712A (en) * | 2011-03-10 | 2011-07-20 | 东南大学 | Triangular net model constructing method based on multi-stratum and three-dimensional quantity of earth work |
| WO2013057760A1 (en) * | 2011-10-19 | 2013-04-25 | Hitachi, Ltd. | Digital elevation model generating system and method of generating a digital elevation model |
| CN103399990A (en) * | 2013-07-18 | 2013-11-20 | 北京工业大学 | Method of constructing fine discrete road grid in urban drainage simulation system |
-
2016
- 2016-06-21 CN CN201811364177.4A patent/CN109523630B/en active Active
- 2016-06-21 CN CN201811364179.3A patent/CN109523631B/en active Active
- 2016-06-21 CN CN201610454545.9A patent/CN106097450B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0562449A2 (en) * | 1992-03-19 | 1993-09-29 | Matsushita Electric Industrial Co., Ltd. | Object shape expression apparatus |
| WO1996013812A1 (en) * | 1994-10-29 | 1996-05-09 | Masaharu Ishii | Apparatus and method for terrain model reproduction |
| JP2005234787A (en) * | 2004-02-18 | 2005-09-02 | Yoshiaki Goto | Visualization technique using cg |
| KR100898617B1 (en) * | 2008-11-20 | 2009-05-27 | 주식회사 범아엔지니어링 | Construction method for digital elevation model of area coexisting the ground and water through verification of tin data of lidar and mbes measure value |
| CN102314711A (en) * | 2010-07-01 | 2012-01-11 | 中国地质科学院矿产资源研究所 | Three-dimensional visualization method and device for mineral resource evaluation information |
| CN102682479A (en) * | 2012-04-13 | 2012-09-19 | 国家基础地理信息中心 | Method for generating three-dimensional terrain feature points on irregular triangulation network |
| CN103226845A (en) * | 2013-03-18 | 2013-07-31 | 北京市测绘设计研究院 | TIN-based earth surface refinement processing method and system |
| CN103530901A (en) * | 2013-10-17 | 2014-01-22 | 中国测绘科学研究院 | Method and system for visualizing dynamic water area by matching with real three-dimensional environment |
| CN103871102A (en) * | 2014-03-28 | 2014-06-18 | 南京大学 | Road three-dimensional fine modeling method based on elevation points and road outline face |
| CN104331389A (en) * | 2014-11-21 | 2015-02-04 | 南京南瑞集团公司 | Eight-point method-based contour-line tracing algorithm |
| CN105303590A (en) * | 2015-11-25 | 2016-02-03 | 中国地质大学（武汉） | Method for cutting block model of three-dimensional ore body with attribute to generate two-dimensional profile map |
Non-Patent Citations (3)
| Title |
|---|
| QIAO CHAOFEI等: "Study of the Spatial Relations of Topographic Contour Lines", 《GEOMATICS AND SPATIAL GEOGRAPHIC INFORMATION》 * |
| 张尧等: "基于Delaunay三角网的等高线树生成方法", 《测绘学报》 * |
| 张瑞军: "遥感图像三维可视化在康定机场建设中的应用", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106557573B (en) * | 2016-11-24 | 2019-11-19 | 中铁第五勘察设计院集团有限公司 | The definition of irregular triangle network data memory format and data processing method based on directory index structure |
| CN106557573A (en) * | 2016-11-24 | 2017-04-05 | 中铁第五勘察设计院集团有限公司 | The definition of TIN data memory format and data processing method based on directory index structure |
| CN106592357A (en) * | 2016-12-14 | 2017-04-26 | 中国人民解放军空军工程大学 | Airfield pavement curve transitional surface elevation designing method |
| CN106592357B (en) * | 2016-12-14 | 2018-10-23 | 中国人民解放军空军工程大学 | A kind of airfield pavement bend transition face elevation design method |
| CN107067446A (en) * | 2017-01-19 | 2017-08-18 | 中国测绘科学研究院 | A kind of contour proximity relations judges and direction automatic adjusting method |
| CN106909740A (en) * | 2017-03-02 | 2017-06-30 | 中国石油天然气集团公司 | The computational methods and computing device of a kind of Oil ＆ Gas Storage yard total figure volume of earthwork |
| CN106909740B (en) * | 2017-03-02 | 2020-04-17 | 中国石油天然气集团有限公司 | Method and device for calculating total map earthwork amount of oil and gas storage and transportation station |
| CN107145669B (en) * | 2017-05-08 | 2020-07-03 | 攀钢集团矿业有限公司 | Rock mass calculating method and device |
| CN107145669A (en) * | 2017-05-08 | 2017-09-08 | 攀钢集团矿业有限公司 | A kind of ore-rock amount computational methods and its device |
| CN109255008A (en) * | 2018-09-28 | 2019-01-22 | 贵州省水利水电勘测设计研究院 | A kind of data processing method, system and electronic equipment |
| CN109255008B (en) * | 2018-09-28 | 2020-09-25 | 贵州省水利水电勘测设计研究院有限公司 | Data processing method and system and electronic equipment |
| US11618477B2 (en) * | 2020-07-16 | 2023-04-04 | Waymo Llc | Service area maps for autonomous vehicles |
| CN112419495A (en) * | 2020-10-26 | 2021-02-26 | 天津大学 | Elevation point automatic extraction method based on multi-scale DEM space model |
| CN114046776A (en) * | 2021-09-22 | 2022-02-15 | 北京洛斯达科技发展有限公司 | Power transmission engineering water and soil conservation measure implementation checking system |
| CN114046776B (en) * | 2021-09-22 | 2023-04-21 | 北京洛斯达科技发展有限公司 | Verification system for implementing water and soil conservation measures of power transmission engineering |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106097450B (en) | 2019-01-22 |
| CN109523630B (en) | 2022-07-29 |
| CN109523631A (en) | 2019-03-26 |
| CN109523630A (en) | 2019-03-26 |
| CN109523631B (en) | 2022-08-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106097450A (en) | A kind of contour lines creation method and device | |
| CN107025685B (en) | Airborne building roof point cloud modeling method under topology perception | |
| CN106446910B (en) | Complex geological curved surface feature extraction and reconstruction method | |
| De Floriani et al. | Applications of computational geometry to geographic information systems. | |
| CN100514369C (en) | Non-homogeneous space partition based scene visibility cutting method | |
| CN101872376A (en) | True three-dimensional road intelligent design method and system | |
| CN104966317B (en) | A kind of three-dimensional method for automatic modeling based on ore body contour line | |
| CN104766366A (en) | Method for establishing three-dimensional virtual reality demonstration | |
| CN104835202A (en) | Quick three-dimensional virtual scene constructing method | |
| CN102663800A (en) | City building complex and rendering method considering city image | |
| CN102663811A (en) | Real-time three-dimensional visualization method of large-scale scene city building based on spatial cognition | |
| CN102538754B (en) | Method and equipment for calculating road gradient information | |
| CN105184867A (en) | Three-dimensional stratum modeling method based on regular grids and corner-point grid technology | |
| CN104978763B (en) | A kind of simultaneous synthesis map emulation mode of network of waterways key element and DEM based on three-dimensional Douglas Peucker algorithms | |
| CN111858810B (en) | Modeling elevation point screening method for road DEM construction | |
| CN109035364A (en) | Method for rapidly drawing profile map based on CAD topographic map | |
| CN102651141A (en) | Method for automatically superposing and rendering vector map and high-resolution digital elevation models in large scene | |
| CN108198143A (en) | A kind of dimensional topography avoiding cracks method that LoD levels is gone to constrain | |
| CN106777779A (en) | A kind of railway culvert method for designing based on BIM | |
| CN109118588A (en) | A kind of colored LOD model automatic forming method decomposed based on block | |
| CN107918957A (en) | A kind of three-dimensional building model simplifying method of holding structure and textural characteristics | |
| CN107220372A (en) | A kind of automatic laying method of three-dimensional map line feature annotation | |
| Wang et al. | A hybrid GIS for 3-D city models | |
| CN101694726A (en) | Fusing and drawing method based on multi-source terrain data | |
| CN102567594B (en) | Method for simulation modeling of offshore island reef type artificial fish reef cluster flow field |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |