CN102339478A - Method and device for generating digital elevation model from contour map - Google Patents

Method and device for generating digital elevation model from contour map Download PDF

Info

Publication number
CN102339478A
CN102339478A CN2011102724447A CN201110272444A CN102339478A CN 102339478 A CN102339478 A CN 102339478A CN 2011102724447 A CN2011102724447 A CN 2011102724447A CN 201110272444 A CN201110272444 A CN 201110272444A CN 102339478 A CN102339478 A CN 102339478A
Authority
CN
China
Prior art keywords
grid
height value
sampled point
elevation
contour map
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2011102724447A
Other languages
Chinese (zh)
Other versions
CN102339478B (en
Inventor
史明昌
黄兆伟
李团宏
曹刚
孙成宝
李嵩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Datum Science & Technology Development Co Ltd
Original Assignee
Beijing Datum Science & Technology Development Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Datum Science & Technology Development Co Ltd filed Critical Beijing Datum Science & Technology Development Co Ltd
Priority to CN2011102724447A priority Critical patent/CN102339478B/en
Publication of CN102339478A publication Critical patent/CN102339478A/en
Application granted granted Critical
Publication of CN102339478B publication Critical patent/CN102339478B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Processing Or Creating Images (AREA)
  • Instructional Devices (AREA)

Abstract

The application provides a method and a device for generating a digital elevation model from a contour map, which relate to the field of geographic information. The method comprises the steps that: a rectangular enclosing box of the contour map is divided into a grid according to the longitudinal and horizontal resolution of the required digital elevation model; the elevation values of the grid in which a sampling point is in and a plurality of grids around the grid are obtained according to the elevation changes of each sampling point on each contour line in the contour map and the elevation value of the sampling point; the whole grid is browsed, the elevation values of the treated grids which are most adjacent to the grid in a plurality of directions are searched for the grid which is not treated, and the average value of the values of the treated grids is adopted as the elevation value of the grid which is not treated; and the whole digital elevation model is output according to the elevation value of the whole grid. When the elevation value of the grid is obtained through the application, a dynamic effect model (DEM) is generated by direct calculation without selecting a reference point set, so that the computing speed and efficiency are obviously improved.

Description

A kind of method and apparatus that generates digital elevation model from contour map
Technical field
The application relates to the geography information field, particularly relates to a kind of method and apparatus that generates digital elevation model from contour map.
Background technology
Digital elevation model (Digital Elevation Model; Abbreviation DEM) be the planimetric coordinates (X of regular grid point in the certain limit; Y) and the data set of elevation (Z); It mainly is a space distribution of describing regional landform shape, is to carry out data acquisition (comprising sampling and measurement) through level line or similar stereoscopic model, carries out interpolation of data then and forms.DEM is the virtual representation to landform shape, can derive information such as level line, slope map, also can superpose with remote sensing image or other thematic data, is used for the analytical applications relevant with landform, the basic data that the while itself is still made remote sensing image.Digital elevation model DEM is a kind of basic data of GIS-Geographic Information System, can carry out 3-D display and carry out terrain analysis on this basis.Contour map has crucial research and using value as a kind of important and cheap Data Source that generates DEM.
In the prior art; The common practice that generates DEM from contour map is to adopt the interpolation algorithm of DEM; Interior slotting total thought of DEM is to obtain the height value of point to be calculated according to the elevation of several known reference point, and concrete interpolation algorithm develops into a variety of in practical application again; Concrete, generate DEM from contour map, for the grid (being pixel) of height value to be calculated, according to prior art; Need from a series of known level lines, select the reference point set S1 of influence, S2; ..., Sn, and then carry out interior inserting with these RPs and ask height value.Different algorithms is selected with reference to the method for point set difference to some extent, but when DEM to be generated is mass data, when the contour map data volume is very big again, this select with reference to the process of point set consuming time very high.
Summary of the invention
The application's technical matters to be solved provides a kind of method and apparatus that generates digital elevation model from contour map; Utilize on the contour map basis; To calculating the grid of height value, need not select with reference to the point set process, directly calculate and generate DEM; Under the situation that guarantees the good quality of data, computing velocity efficient is significantly improved.
In order to address the above problem, the application discloses the method that a kind of contour map generates digital elevation model, comprising:
Grid generates step, and the vertical and horizontal resolution of digital elevation model as required is divided into grid with the quadrilateral bounding box of contour map;
The first height value obtaining step; To each sampled point on each level line in the said contour map; According to the elevation change rate of said sampled point and the height value of said sampled point, obtain the grid at said sampled point place and the height value of a plurality of grids around this grid;
The second height value obtaining step; Travel through whole grid; To being in the grid of the state of being untreated, to search and the said grid height value of the grid of having handled of contiguous a plurality of directions, the average of value of getting said a plurality of grids of having handled is as the height value of the said grid that is untreated;
Digital elevation model output step is exported whole digital elevation model according to the height value of the whole grid that is obtained.
Preferably, carry out following steps for each sampled point:
Step 121 to the sampled point on the level line in the said contour map, according to isocontour spacing and the difference of elevation nearest with said sampled point, obtains the elevation change rate of said sampled point;
Step 122 according to the grid at said sampled point place, is confirmed a plurality of grids that said grid is contiguous on every side;
Step 123, obtain said sampled point and said around the distance of each grid element center point in contiguous a plurality of grids; And whether judge said distance greater than threshold value, if less than threshold value, then change step 124 over to, if, change step 125 over to then greater than threshold value;
Step 124 is with the height value of this sampled point height value as this grid;
Step 125 according to the elevation change rate of said distance and corresponding sampled point, obtains the elevation change amount of corresponding said grid, and judges that said elevation change increases progressively or successively decreases, if increase progressively, then changes step 126 over to, if successively decrease, then changes step 127 over to;
Step 126 adds result that the elevation change amount of corresponding said grid the obtains height value as said grid with the height value of said corresponding sampled point;
Step 127 deducts result that the elevation change amount of corresponding said grid the obtains height value as said grid with the height value of said corresponding sampled point.
Preferably, when a grid is handled, carry out following steps:
Judge whether this grid is handled: when this grid is untreated, with the height value of this grid that obtains directly as the height value of this grid; When this grid has been handled, the height value and the existing height value of this grid of this grid that obtains are got the height value of average as this grid.
Preferably, judge in the step 124 that elevation change increases progressively or the process of successively decreasing comprises:
Steps A 1 confirms that the central spot of said grid is inner still outside in this sampled point place level line;
Steps A 2 judges that according to said contour map distribution trend said central point increases progressively or successively decreases.
Preferably, in step 123, also comprise: for sampled point place grid, directly with the height value of sampled point height value as this grid.
Preferably, said steps A 2 specifically comprises:
When contour map from innermost layer when outermost height raises, if said central point is inner at said high line, then elevation change is successively decreased;
When contour map from innermost layer when outermost height reduces, if said central point is inner at said high line, then elevation change increases progressively.
Corresponding the application also discloses a kind of device from contour map generation digital elevation model, it is characterized in that, comprising:
The grid generation module, digital elevation model vertical and horizontal resolution as required is divided into grid with the minimum quadrilateral of contour map;
The first height value acquisition module; Be used for to each sampled point on each level line of said contour map; According to the elevation change rate of said sampled point and the height value of said sampled point, obtain the grid at said sampled point place and the height value of a plurality of grids around this grid;
The second height value acquisition module; Be used to travel through whole grid; To being in the grid of the state of being untreated, to search and the said grid height value of the grid of having handled of contiguous a plurality of directions, the average of value of getting said a plurality of grids of having handled is as the height value of the said grid that is untreated;
The digital elevation model output module is used for exporting whole digital elevation model according to the height value of whole grid.
Preferably, the said first height value acquisition module specifically comprises:
The elevation change rate is obtained submodule, is used for according to isocontour spacing and the difference of elevation nearest with said sampled point, obtaining the elevation change rate of said sampled point to the sampled point on the level line of said contour map;
Grid obtains submodule, is used for the grid according to said sampled point place, confirms a plurality of grids that said grid is contiguous on every side;
Distance is obtained judge module, is used for obtaining said sampled point and the said distance of each grid element center point of contiguous a plurality of grids on every side; And whether judge said distance greater than threshold value, and if less than threshold value, then change height value over to and obtain submodule one, if greater than threshold value, change the elevation change amount over to and obtain the judgement submodule;
Height value is obtained submodule one, with the height value of this sampled point height value as this grid;
The elevation change amount is obtained the judgement submodule; Be used for elevation change rate, obtain the elevation change amount of corresponding said grid, and judge that said elevation change increases progressively or successively decreases according to said distance and corresponding sampled point; If increase progressively; Then change height value over to and obtain submodule two,, then change height value over to and obtain submodule three if successively decrease;
Height value is obtained submodule two, is used for result that height value with said corresponding sampled point adds that the elevation change amount of corresponding said grid the obtains height value as said grid;
Height value is obtained submodule three, is used for result that elevation change amount that height value with said corresponding sampled point deducts corresponding said grid the obtains height value as said grid.
Preferably, also comprise overlapping processing module, be used for when this grid is untreated, with the height value of this grid that obtains directly as the height value of this grid; When this grid has been handled, the height value and the existing height value of this grid of this grid that obtains are got the height value of average as this grid.
Preferably, the said second height value acquisition module also comprises:
First submodule confirms that the central spot of said grid is inner still outside in this sampled point place level line;
Second submodule judges that according to said contour map distribution trend said central point increases progressively or successively decreases.
Compared with prior art, the application comprises following advantage:
The application passes through at first to the elevation change rate of each known sampled point and the height value of this sampled point; Obtain the said sampled point place grid and the height value of a plurality of grids on every side thereof; And then travel through whole grid; To being in the grid of the state of being untreated; Search and the said grid height value of the grid of having handled of contiguous a plurality of directions, the average of value of getting said a plurality of grids of having handled is as the height value of the said grid that is untreated, and the process that obtains whole grid height value so need not selected directly to calculate with reference to the process of point set and generate DEM; And can be under the situation that guarantees the good quality of data, computing velocity efficient is significantly improved.
Description of drawings
Fig. 1 is a kind of schematic flow sheet that generates the method for digital elevation model from contour map of the application;
Fig. 2 is the actual level line legend of of the application;
Fig. 3 is a kind of schematic flow sheet that generates the preferred first height value obtaining step of method of digital elevation model from contour map of the application;
Fig. 4 is a kind of preferred schematic flow sheet of further handling a grid of method that generates digital elevation model from contour map of the application;
Fig. 5 is the digital elevation model figure that the application is generated by actual level line legend;
Fig. 6 is a kind of structural representation that generates the device of digital elevation model from contour map of the application;
Fig. 7 is a kind of structural representation that generates the preferred first height value acquisition module of device of digital elevation model from contour map of the application.
Embodiment
For above-mentioned purpose, the feature and advantage that make the application can be more obviously understandable, the application is done further detailed explanation below in conjunction with accompanying drawing and embodiment.
The application passes through at first to the elevation change rate of each known sampled point and the height value of this sampled point; Obtain the said sampled point place grid and the height value of a plurality of grids on every side thereof; And then travel through whole grid; To being in the grid of the state of being untreated; Search and the said grid height value of the grid of having handled of contiguous a plurality of directions, the average of value of getting said a plurality of grids of having handled is as the height value of the said grid that is untreated, and the process that obtains whole grid height value so need not selected directly to calculate with reference to the process of point set and generate DEM; And can be under the situation that guarantees the good quality of data, computing velocity efficient is significantly improved.
With reference to Fig. 1, show that the application is a kind of to generate the schematic flow sheet of the method for digital elevation model from contour map, its step comprises:
Grid generates step 110, and the vertical and horizontal resolution of digital elevation model as required is divided into grid with the quadrilateral bounding box of contour map.
After obtaining contour map; General at first through calculating the bounding box of contour map; Obtain comprising the minimum quadrilateral of these high line charts; And then the resolution division grid number vertical, horizontal according to DEM, wherein when dividing grid, be the basis with minimum quadrilateral, progressively adjustment meets the quadrilateral of the encirclement contour map of an integer grid.After this step is being divided the grid number, can these trellis state be set to " being untreated ".
With reference to Fig. 2; It is the application's an actual level line legend, and these high line charts are contour maps of a hillside type, and its isocontour value is successively decreased to outermost layer by innermost layer gradually; Wherein the level line difference of elevation is 100 meters, and the isocontour height value of innermost layer is 1000 meters.
Comprising all isocontour minimum quadrilaterals through calculating, also is the frame rectangle among Fig. 2, obtains the long a and the wide b of this rectangle simultaneously, the demand of DEM resolution according to actual needs; The resolution of DEM is confirmed as 50 meters, and the sizing grid that promptly obtains DEM is a/50, b/50, wherein; This result calculated just is an integer not necessarily, and such as a=110, but number of grid is integer certainly, so; The smallest positive integral 2 of getting a/50=2.2 again+1, promptly we get 3, readjust the long a value of rectangle like this; Be a=3*50=150 rice now, carried out aforesaid operations adjustment b value with same principle, finally made a/50 and b/50 be integer for b/50.Here the original state of each grid is " being untreated ".
The first height value obtaining step 120; To each sampled point on each level line in the said contour map; According to the elevation change rate of said sampled point and the height value of said sampled point, obtain the grid at said sampled point place and the height value of a plurality of grids around this grid.
Have known a plurality of sampled points in the reality on each level line of contour map.To level line CL all in the contour map 1, CL 2..., CL nCarry out circular treatment, for level line CL i(n>=i>=1) is to all the sampled point P on it 1, P 2..., P mCircular treatment.Wherein, P no matter iWhere (m>=i>=1) is, its must be in a certain grid of DEM or the limit on, be the center with this grid; It has a plurality of grids on every side; Preferably, getting with a grid among the application is the center, adjacent 8 grids around it (upper and lower, left and right, upper left, a left side is following, upper right, totally 8 of bottom rights); Can obtain one 3 * 3 macrolattice A, calculate the height value of 9 grids in macrolattice A the inside then according to the height value of the elevation change rate of this sampled point and said sampled point.
Further, with reference to Fig. 3, being that the application is a kind of generates the schematic flow sheet of the preferred first height value obtaining step of method of digital elevation model from contour map, specifically comprises:
Step 121 to the sampled point on the level line in the said contour map, according to isocontour spacing and the difference of elevation nearest with said sampled point, obtains the elevation change rate of said sampled point.
For the sampled point (P on certain level line x, P y); On two level lines adjacent, find that level line nearest, and calculate described minimum distance L, in conjunction with the difference of elevation d between the level line with this sampled point with this sampled point place level line; Obtain the elevation change rate in this shop, generally obtain elevation change rate: T through following formula j=d/L, wherein T jThe elevation change rate of representing this sampled point.
Step 122 according to the grid at said sampled point place, is confirmed a plurality of grids that said grid is contiguous on every side.
After whole contour map is vertical according to DEM in step 110, lateral resolution is divided into a plurality of grids; Each sampled point on every level line all can be at a grid on the inner or limit; When its when a grid is inner, be f grid around the center grid is chosen directly with this grid; When its on the limit of a grid time, picked at random one of them grid of the grid of limit or concurrent is altogether chosen the grid around it as the center grid; The application preferably gets adjacent 8 grids on every side of center grid, can obtain one 3 * 3 macrolattice A.Such as among Fig. 2, have 35 sampled points in the level line of innermost layer, to each sampled point P, (x y) calculates (P in the position of DEM grid by its locus P x, P y), the grid that comprises of macrolattice A is respectively like this: (P X+1, P Y+1), (P x, P Y+1), (P X+1, P y),
(P x?1,P y?1)、(P x+1,P y?1)、(P x?1,P y+1)、(P x,P y?1)、(P x?1,P y)、(P x,P y)。
Step 123, obtain said sampled point and said around the distance of each grid element center point in contiguous a plurality of grids; And whether judge said distance greater than threshold value, if less than threshold value, then change step 124 over to, if, change step 125 over to then greater than threshold value.
After having chosen grid, the sampled point P in computing center's grid jTo around the distance of grid element center point, after having chosen 3 * 3 macrolattice A, computing grid arrives the distance of this sampled point.For each said grid to this sampled point apart from Q, whether judge it greater than threshold value, if less than threshold value, then change step 124 over to, if greater than threshold value then change step 125 over to; Wherein said threshold value can be set to the size of a grid or the size of half grid.Such as (P X+1, P Y+1) to sampled point P (x, distance y) is Q, but this Q value is less than the length of a grid, the time, change step 124 so over to; But (P x, P Y+1) to sampled point apart from the length of Q value greater than a grid, then change step 125 over to.
Wherein, generally speaking, the center grid at sampled point place is directly with the height value of the sampled point height value as this center grid.
Step 124 is with the height value of this sampled point height value as this grid.
Because this grid is to the sampled point distance grid less than threshold value, so directly with the height value of this sampled point height value as this grid.(x, height value y) is 1000, and (P such as a sampled point P on the aforementioned innermost layer level line X+1, P Y+1) to sampled point P (x, distance y) is the length less than a grid, then with 1000 directly as grid (P X+1, P Y+1) height value.
Step 125 according to the elevation change rate of said distance and corresponding sampled point, obtains the elevation change amount of corresponding said grid, and judges that said elevation change increases progressively or successively decreases, if increase progressively, then changes step 126 over to, if successively decrease, then changes step 127 over to.
For to the grid of sampled point distance,, obtain the elevation change amount of corresponding said grid according to the elevation change rate of said distance and corresponding sampled point greater than threshold value; And judge that said elevation change amount increases progressively or successively decreases, if increase progressively, then changes step 126 over to; If successively decrease, then change step 127 over to.Wherein the elevation change amount can be calculated through following mode: W=Q * T j, wherein W representes elevation change amount that should grid, T jFor aforesaid expression to elevation change rate that should sampled point, Q is the distance that aforesaid this grid element center is put this sampled point.Such as aforementioned (P x, P Y+1) to sampled point apart from the length of Q value greater than a grid, then calculate elevation change amount W that should grid, and judge that said elevation change increases progressively or successively decreases.
Further, can judge that the elevation change amount increases progressively or the process of successively decreasing through following steps:
Steps A 1 confirms that the central spot of said grid is inner still outside in this sampled point place level line;
Steps A 2 judges that according to said contour map distribution trend said central point increases progressively or successively decreases.
Further, said steps A 2 specifically comprises:
When contour map from innermost layer when outermost height raises, if said central point is inner at said high line, then elevation change is successively decreased;
When contour map from innermost layer when outermost height reduces, if said central point is inner at said high line, then elevation change increases progressively.
This process can be to increase progressively the limit or successively decrease the limit at this isocontour elevation according to this grid; Isocontour elevation increases progressively limit definition: if level line internal point elevation greater than on the level line, then the level line internal point belongs to elevation and increases progressively the limit; If level line internal point elevation is less than on the level line, then the level line external point belongs to elevation and increases progressively the limit.
Step 126 adds result that the elevation change amount of corresponding said grid the obtains height value as said grid with the height value of said corresponding sampled point.
Elevation change amount such as a grid is W, and the height value of this sampled point is S, and the elevation change that this grid is corresponding increases progressively, then the height value Z=S+W of this grid.
Step 127 deducts result that the elevation change amount of corresponding said grid the obtains height value as said grid with the height value of said corresponding sampled point.Elevation change amount such as a grid is W, and the height value of this sampled point is S, and the elevation change that this grid is corresponding increases progressively, then the height value Z=S-W of this grid.
Can obtain the center grid at sampled point place and the height value of grid on every side thereof like this.
So round-robin carries out above-mentioned each process to all sampled points and grid, can obtain the height value of each grid of the correspondence of all sampled points on all level lines in all contour maps.
Wherein, with reference to Fig. 4, when a grid is handled, carry out following steps:
Judge whether this grid is handled: when this grid is untreated, with the height value of this grid that obtains directly as the height value of this grid; When this grid has been handled, the height value and the existing height value of this grid of this grid that obtains are got the height value of average as this grid.Such as, to sampled point P 1The height value that obtains this grid after corresponding grid d handles for the first time is g, and sampled point P 3Also comprise grid d in the corresponding grid, obtained P by abovementioned steps again 3Obtain after should another height value h of grid, with the mean value of g and h promptly (g+h)/2 as the new height value of grid d.
The 3rd height value obtaining step 130; Travel through whole grid; To being in the grid of the state of being untreated, to search and the said grid height value of the grid of having handled of contiguous a plurality of directions, the average of value of getting said a plurality of grids of having handled is as the height value of the said grid that is untreated.
In the reality, after the corresponding a plurality of grids of sampled point are handled, also can there be a plurality of untreated grids, be the grid that obtains height value.In order to obtain the height value of these grids, travel through whole grid, for untreated grid, search handling of a plurality of directions that this grid closes on most the height value of grid, get the height value of the mean value of these a plurality of height value as this grid.The application preferably searches 8 directions (maximum 8 are respectively upper and lower, left and right, upper left, upper right, left following, bottom right that this grid closes on most.If on the border, then get in esse several direction) the grid of having handled, the mean value of height value of asking these grids is as the height value of this grid.
Digital elevation model output step 140 is exported whole digital elevation model according to the height value of the whole grid that is obtained.
With reference to Fig. 5, be the final DEM that obtains through aforementioned process by contour map among Fig. 2.As scheme innermost layer and represent that height value is bigger, belonging to the summit, hillside, outermost layer is represented the subregion that elevation is lower.
The application's method has been removed selection and has directly been calculated generation DEM with reference to the process of point set, and can be under the situation that guarantees the good quality of data, and computing velocity efficient is significantly improved.
Disclosed herein as well is a kind of device that generates digital elevation model from contour map accordingly.
With reference to Fig. 6, being that the application is a kind of generates the structural representation of the device of digital elevation model from contour map, and said device comprises:
The grid generation module, digital elevation model vertical and horizontal resolution as required is divided into grid with the quadrilateral bounding box of contour map.
The first height value acquisition module; Be used for to each sampled point on each level line of said contour map; According to the elevation change rate of said sampled point and the height value of said sampled point, obtain the grid at said sampled point place and the height value of a plurality of grids around this grid.
The second height value acquisition module; Be used to travel through whole grid; To being in the grid of the state of being untreated, to search and the said grid height value of the grid of having handled of contiguous a plurality of directions, the average of value of getting said a plurality of grids of having handled is as the height value of the said grid that is untreated.
The digital elevation model output module is used for exporting whole digital elevation model according to the height value of whole grid.
Further, with reference to Fig. 7, be a kind of structural representation that generates the preferred first height value acquisition module of device of digital elevation model from contour map of the application.The said first height value acquisition module specifically comprises:
The elevation change rate is obtained submodule, is used for according to isocontour spacing and the difference of elevation nearest with said sampled point, obtaining the elevation change rate of said sampled point to the sampled point on the level line of said contour map.
Grid obtains submodule, is used for the grid according to said sampled point place, confirms a plurality of grids that said grid is contiguous on every side.
Distance is obtained judge module, is used for obtaining said sampled point and the said distance of each grid element center point of contiguous a plurality of grids on every side; And whether judge said distance greater than threshold value, and if less than threshold value, then change height value over to and obtain submodule one, if greater than threshold value, change the elevation change amount over to and obtain the judgement submodule.
Height value is obtained submodule one, with the height value of this sampled point height value as this grid;
The elevation change amount is obtained the judgement submodule; Be used for elevation change rate, obtain the elevation change amount of corresponding said grid, and judge that said elevation change increases progressively or successively decreases according to said distance and corresponding sampled point; If increase progressively; Then change height value over to and obtain submodule two,, then change height value over to and obtain submodule three if successively decrease.
Further, the elevation change amount is obtained and is judged that submodule comprises:
Judge submodule one, be used for confirming that the central spot of said grid is inner still outside in this sampled point place level line.
Judge submodule two, be used for judging that according to said contour map distribution trend said central point increases progressively or successively decreases.
Height value is obtained submodule two, is used for result that height value with said corresponding sampled point adds that the elevation change amount of corresponding said grid the obtains height value as said grid.
Height value is obtained submodule three, is used for result that elevation change amount that height value with said corresponding sampled point deducts corresponding said grid the obtains height value as said grid.
Further, also comprise overlapping processing module, be used for when this grid is untreated, with the height value of this grid that obtains directly as the height value of this grid; When this grid has been handled, the height value and the existing height value of this grid of this grid that obtains are got the height value of average as this grid.
Further, the said second height value acquisition module also comprises:
First submodule confirms that the central spot of said grid is inner still outside in this sampled point place level line.
Second submodule judges that according to said contour map distribution trend said central point increases progressively or successively decreases.
For system embodiment, because it is similar basically with method embodiment, so description is fairly simple, relevant part gets final product referring to the part explanation of method embodiment.
Each embodiment in this instructions all adopts the mode of going forward one by one to describe, and what each embodiment stressed all is and the difference of other embodiment that identical similar part is mutually referring to getting final product between each embodiment.
More than to the application provided a kind of method and apparatus that generates digital elevation model from contour map; Carried out detailed introduction; Used concrete example among this paper the application's principle and embodiment are set forth, the explanation of above embodiment just is used to help to understand the application's method and core concept thereof; Simultaneously, for one of ordinary skill in the art, according to the application's thought, the part that on embodiment and range of application, all can change, in sum, this description should not be construed as the restriction to the application.

Claims (10)

1. the method from contour map generation digital elevation model is characterized in that, comprising:
Grid generates step, and the vertical and horizontal resolution of digital elevation model as required is divided into grid with the quadrilateral bounding box of contour map;
The first height value obtaining step; To each sampled point on each level line in the said contour map; According to the elevation change rate of said sampled point and the height value of said sampled point, obtain the grid at said sampled point place and the height value of a plurality of grids around this grid;
The second height value obtaining step; Travel through whole grid; To being in the grid of the state of being untreated, to search and the said grid height value of the grid of having handled of contiguous a plurality of directions, the average of value of getting said a plurality of grids of having handled is as the height value of the said grid that is untreated;
Digital elevation model output step is exported whole digital elevation model according to the height value of the whole grid that is obtained.
2. method according to claim 1 is characterized in that, carries out following steps for each sampled point:
Step 121 to the sampled point on the level line in the said contour map, according to isocontour spacing and the difference of elevation nearest with said sampled point, obtains the elevation change rate of said sampled point;
Step 122 according to the grid at said sampled point place, is confirmed a plurality of grids that said grid is contiguous on every side;
Step 123, obtain said sampled point and said around the distance of each grid element center point in contiguous a plurality of grids; And whether judge said distance greater than threshold value, if less than threshold value, then change step 124 over to, if, change step 125 over to then greater than threshold value;
Step 124 is with the height value of this sampled point height value as this grid;
Step 125 according to the elevation change rate of said distance and corresponding sampled point, obtains the elevation change amount of corresponding said grid, and judges that said elevation change increases progressively or successively decreases, if increase progressively, then changes step 126 over to, if successively decrease, then changes step 127 over to;
Step 126 adds result that the elevation change amount of corresponding said grid the obtains height value as said grid with the height value of said corresponding sampled point;
Step 127 deducts result that the elevation change amount of corresponding said grid the obtains height value as said grid with the height value of said corresponding sampled point.
3. method according to claim 1 and 2 is characterized in that, when a grid is handled, carries out following steps:
Judge whether this grid is handled: when this grid is untreated, with the height value of this grid that obtains directly as the height value of this grid; When this grid has been handled, the height value and the existing height value of this grid of this grid that obtains are got the height value of average as this grid.
4. method according to claim 2 is characterized in that:
Judge in the step 124 that elevation change increases progressively or the process of successively decreasing comprises:
Steps A 1 confirms that the central spot of said grid is inner still outside in this sampled point place level line;
Steps A 2 judges that according to said contour map distribution trend said central point increases progressively or successively decreases.
5. method according to claim 2 is characterized in that:
In step 123, also comprise: for sampled point place grid, directly with the height value of sampled point height value as this grid.
6. method according to claim 4 is characterized in that:
Said steps A 2 specifically comprises:
When contour map from innermost layer when outermost height raises, if said central point is inner at said high line, then elevation change is successively decreased;
When contour map from innermost layer when outermost height reduces, if said central point is inner at said high line, then elevation change increases progressively.
7. the device from contour map generation digital elevation model is characterized in that, comprising:
The grid generation module, digital elevation model vertical and horizontal resolution as required is divided into grid with the minimum quadrilateral of contour map;
The first height value acquisition module; Be used for to each sampled point on each level line of said contour map; According to the elevation change rate of said sampled point and the height value of said sampled point, obtain the grid at said sampled point place and the height value of a plurality of grids around this grid;
The second height value acquisition module; Be used to travel through whole grid; To being in the grid of the state of being untreated, to search and the said grid height value of the grid of having handled of contiguous a plurality of directions, the average of value of getting said a plurality of grids of having handled is as the height value of the said grid that is untreated;
The digital elevation model output module is used for exporting whole digital elevation model according to the height value of whole grid.
8. device according to claim 7 is characterized in that:
The said first height value acquisition module specifically comprises:
The elevation change rate is obtained submodule, is used for according to isocontour spacing and the difference of elevation nearest with said sampled point, obtaining the elevation change rate of said sampled point to the sampled point on the level line of said contour map;
Grid obtains submodule, is used for the grid according to said sampled point place, confirms a plurality of grids that said grid is contiguous on every side;
Distance is obtained judge module, is used for obtaining said sampled point and the said distance of each grid element center point of contiguous a plurality of grids on every side; And whether judge said distance greater than threshold value, and if less than threshold value, then change height value over to and obtain submodule one, if greater than threshold value, change the elevation change amount over to and obtain the judgement submodule;
Height value is obtained submodule one, with the height value of this sampled point height value as this grid;
The elevation change amount is obtained the judgement submodule; Be used for elevation change rate, obtain the elevation change amount of corresponding said grid, and judge that said elevation change increases progressively or successively decreases according to said distance and corresponding sampled point; If increase progressively; Then change height value over to and obtain submodule two,, then change height value over to and obtain submodule three if successively decrease;
Height value is obtained submodule two, is used for result that height value with said corresponding sampled point adds that the elevation change amount of corresponding said grid the obtains height value as said grid;
Height value is obtained submodule three, is used for result that elevation change amount that height value with said corresponding sampled point deducts corresponding said grid the obtains height value as said grid.
9. according to claim 7 or 8 described devices, it is characterized in that:
Also comprise overlapping processing module, be used for when this grid is untreated, with the height value of this grid that obtains directly as the height value of this grid; When this grid has been handled, the height value and the existing height value of this grid of this grid that obtains are got the height value of average as this grid.
10. device according to claim 8 is characterized in that:
The said second height value acquisition module also comprises:
First submodule confirms that the central spot of said grid is inner still outside in this sampled point place level line;
Second submodule judges that according to said contour map distribution trend said central point increases progressively or successively decreases.
CN2011102724447A 2011-09-14 2011-09-14 Method and device for generating digital elevation model from contour map Active CN102339478B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011102724447A CN102339478B (en) 2011-09-14 2011-09-14 Method and device for generating digital elevation model from contour map

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011102724447A CN102339478B (en) 2011-09-14 2011-09-14 Method and device for generating digital elevation model from contour map

Publications (2)

Publication Number Publication Date
CN102339478A true CN102339478A (en) 2012-02-01
CN102339478B CN102339478B (en) 2013-11-13

Family

ID=45515190

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011102724447A Active CN102339478B (en) 2011-09-14 2011-09-14 Method and device for generating digital elevation model from contour map

Country Status (1)

Country Link
CN (1) CN102339478B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103778191A (en) * 2014-01-03 2014-05-07 南京师范大学 Vector contour line data partitioning method with space proximity relation considered
CN104331928A (en) * 2014-10-28 2015-02-04 中国航空工业集团公司洛阳电光设备研究所 Automatic contour elevation value assignment method based on triangulation network
CN105046267A (en) * 2015-05-19 2015-11-11 广州杰赛科技股份有限公司 Dynamic grid generating method, generating apparatus, clustering method and clustering apparatus
CN107393002A (en) * 2017-07-28 2017-11-24 北京金风科创风电设备有限公司 Apparatus and method for extracting terrain boundary
CN108280880A (en) * 2018-01-24 2018-07-13 长春工程学院 A method of improving the digital elevation data resolution of massif using remote sensing image
CN108898681A (en) * 2018-06-26 2018-11-27 中煤航测遥感集团有限公司 Digital elevation model processing method and processing device
CN109191579A (en) * 2018-07-02 2019-01-11 广东省地质测绘院 Method, apparatus, computer equipment and the storage medium of coordinate conversion
CN112419495A (en) * 2020-10-26 2021-02-26 天津大学 Elevation point automatic extraction method based on multi-scale DEM space model

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101477711A (en) * 2009-01-20 2009-07-08 武汉数字媒体工程技术有限公司 Method for automatically generating tri-dimensional topography in Torque engine

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101477711A (en) * 2009-01-20 2009-07-08 武汉数字媒体工程技术有限公司 Method for automatically generating tri-dimensional topography in Torque engine

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
傅炜: "GIS数字地形模型的建立及应用研究", 《陕西师范大学学报(自然科学版)》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103778191B (en) * 2014-01-03 2017-02-15 南京师范大学 Vector contour line data partitioning method with space proximity relation considered
CN103778191A (en) * 2014-01-03 2014-05-07 南京师范大学 Vector contour line data partitioning method with space proximity relation considered
CN104331928A (en) * 2014-10-28 2015-02-04 中国航空工业集团公司洛阳电光设备研究所 Automatic contour elevation value assignment method based on triangulation network
CN104331928B (en) * 2014-10-28 2017-10-03 中国航空工业集团公司洛阳电光设备研究所 A kind of contour automatic labeling method based on the triangulation network
CN105046267A (en) * 2015-05-19 2015-11-11 广州杰赛科技股份有限公司 Dynamic grid generating method, generating apparatus, clustering method and clustering apparatus
CN105046267B (en) * 2015-05-19 2018-08-28 广州杰赛科技股份有限公司 Generation method, generating means, polymerization and the polyplant of dynamic grid
US11227433B2 (en) 2017-07-28 2022-01-18 Beijing Goldwind Science & Creation Windpower Equipment Co., Ltd. Device and method for extracting terrain boundary
CN107393002A (en) * 2017-07-28 2017-11-24 北京金风科创风电设备有限公司 Apparatus and method for extracting terrain boundary
CN108280880A (en) * 2018-01-24 2018-07-13 长春工程学院 A method of improving the digital elevation data resolution of massif using remote sensing image
CN108898681A (en) * 2018-06-26 2018-11-27 中煤航测遥感集团有限公司 Digital elevation model processing method and processing device
CN108898681B (en) * 2018-06-26 2022-07-05 中煤航测遥感集团有限公司 Digital elevation model processing method and device
CN109191579A (en) * 2018-07-02 2019-01-11 广东省地质测绘院 Method, apparatus, computer equipment and the storage medium of coordinate conversion
CN112419495A (en) * 2020-10-26 2021-02-26 天津大学 Elevation point automatic extraction method based on multi-scale DEM space model
CN112419495B (en) * 2020-10-26 2022-11-15 天津大学 Elevation point automatic extraction method based on multi-scale DEM space model

Also Published As

Publication number Publication date
CN102339478B (en) 2013-11-13

Similar Documents

Publication Publication Date Title
CN102339478A (en) Method and device for generating digital elevation model from contour map
CN105336003B (en) The method for drawing out three-dimensional terrain model with reference to the real-time smoothness of GPU technologies
CN105869202B (en) The digital terrain model data creation method and device of electronic map
CN102930509B (en) Intelligent filtering method for airborne laser point cloud data
CN108629231A (en) Obstacle detection method, device, equipment and storage medium
CN108648271B (en) Interpolation method for generating complex terrain grid model based on GIS data
KR20100136604A (en) Real-time visualization system of 3 dimension terrain image
CN104821013A (en) Method and system for specific surface area extraction based on geodetic coordinate system digital elevation model
CN105957134A (en) Systems and methods for 3-d scene acceleration structure creation and updating
CN111047682B (en) Three-dimensional lane model generation method and system
CN108053477B (en) Numerical processing method for deformation in pipeline
CN107507127B (en) Global matching method and system for multi-viewpoint three-dimensional point cloud
US20110050689A1 (en) Method and Apparatus for Creating of 3D Direction Displaying
CN112084280B (en) Multi-scale terrain cutting and splicing method
CN103077549A (en) Real-time large-scale terrain visualization implementation method based on kd tree
CN104751733B (en) The region method for drafting and device of map, path distance sorting technique and system
CN107564098A (en) A kind of Fast rendering method of big Local Area Network three-dimensional noise map
CN102270236A (en) Rasterized geographic information system (GIS)-based spatial relationship judging method and system
CN105844224A (en) Point cloud fast ordering method for on-vehicle LiDAR road points
CN113286982A (en) System and method for generating, updating and enhancing large-scale high-precision 3D road map and multi-level road map
CN105631939B (en) A kind of three-dimensional point cloud distortion correction method and its system based on curvature filtering
CN105005964A (en) Video sequence image based method for rapidly generating panorama of geographic scene
CN102117494A (en) Method and device for delaminating and blocking data in three-dimensional modeling
CN104572924A (en) Multiscale expression information generating method for GIS (geographic information system) vector building polygon
CN107992588B (en) Terrain display system based on elevation tile data

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant