CN102339478B - 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 PDFInfo
- Publication number
- CN102339478B CN102339478B CN2011102724447A CN201110272444A CN102339478B CN 102339478 B CN102339478 B CN 102339478B CN 2011102724447 A CN2011102724447 A CN 2011102724447A CN 201110272444 A CN201110272444 A CN 201110272444A CN 102339478 B CN102339478 B CN 102339478B
- 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.)
- Active
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
Technical field
The application relates to area of geographic information, particularly relates to a kind of method and apparatus from the contour map generating digital elevation model.
Background technology
Digital elevation model (Digital Elevation Model, abbreviation DEM) be the planimetric coordinates (X of regular grid point in certain limit, Y) and the data set of elevation (Z), it is mainly the space distribution of describing the regional landforms form, be to carry out data acquisition (comprise sampling and measure) by level line or similar stereoscopic model, then carry out interpolation of data and form.DEM is the virtual representation to landform shape, can derive the 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 to 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 Geographic Information System, can carry out 3-D display and carry out on this basis terrain analysis.Contour map, as a kind of important and cheap Data Source that generates DEM, has very important research and using value.
In prior art, the common practice that generates DEM from contour map is to adopt the interpolation algorithm of DEM, the total thought of the interpolation 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 again a variety of in actual applications; Concrete, generate DEM from contour map, for the grid (being pixel) of height value to be calculated, according to prior art, need to from a series of known level lines, select the reference point set S1 of impact, S2, ..., Sn, and then carry out interpolation with these reference point and ask height value.Different algorithms, select 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 large 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 is to provide a kind of method and apparatus from the contour map generating digital elevation model, utilize on the contour map basis, to calculating the grid of height value, do not need to select with reference to the point set process, directly generate DEM, in the situation that guarantee the good quality of data, computing velocity efficiency is significantly improved.
In order to address the above problem, the application discloses a kind of method of contour map generating digital elevation model, comprising:
Grid generates step, and the vertical and horizontal resolution of digital elevation model as required, be divided into grid with the quadrilateral bounding box of contour map;
The first height value obtaining step, for each sampled point on each level line in described contour map,, according to the elevation rate of change of described sampled point and the height value of described sampled point, obtain the height value of grid and a plurality of grids around this grid at described sampled point place; Described elevation rate of change is according to isocontour spacing and the difference of elevation acquisition nearest with described sampled point;
Describedly obtain the process of a plurality of grid height value around described sampled point place grid and this grid and specifically comprise:
Obtain the distance of each the grid element center point in described sampled point and the described a plurality of grids that are close on every side; And judge that whether described distance is greater than threshold value, if less than threshold value, with the height value of the height value of this sampled point as this grid, if greater than threshold value, according to the elevation rate of change of described distance and corresponding sampled point, obtain the elevation variable quantity of corresponding described grid, and judge that described elevation variation increases progressively or successively decreases; If increase progressively, the height value of described corresponding sampled point is added result that the elevation variable quantity of corresponding described grid the obtains height value as described grid,, if successively decrease, the height value of described corresponding sampled point is deducted result that the elevation variable quantity of corresponding described grid the obtains height value as described grid;
Described judgement elevation variation increases progressively or the process of successively decreasing comprises: the central spot of determining described grid is inner or outside in this sampled point place level line; Increase progressively or successively decrease according to the described central point of described contour map distribution trend judgement; When contour map raise from innermost layer to outermost height, if described central point is inner at described high line, the elevation variation was successively decreased; When contour map reduced from innermost layer to outermost height, if described central point is inner at described high line, the elevation variation increased progressively;
The second height value obtaining step, travel through whole grid, to being in the grid of the state of being untreated, search the height value of the grid of having processed of the multiple directions the most contiguous with described grid, get the average of value of described a plurality of grids of having processed as the height value of the described grid that is untreated;
When a grid is processed, carry out following steps:
Judge whether this grid is processed: when this grid is untreated, with the height value of this grid of obtaining directly as the height value of this grid; When this grid has been processed, height value and the existing height value of this grid of this grid of acquisition are got the height value of average as this grid;
Digital elevation model output step, according to the whole digital elevation model of height value output of the whole grid that obtains.
Preferably, also comprise in step 123: for sampled point place grid, directly with the height value of the height value of sampled point as this grid.
Corresponding the application also discloses a kind of device from the contour map generating digital elevation model, it is characterized in that, comprising:
The grid generation module, digital elevation model vertical and horizontal resolution as required, become grid with the minimum quadrilateral partition of contour map;
The first height value acquisition module, be used for for each sampled point on each level line of described contour map,, according to the elevation rate of change of described sampled point and the height value of described sampled point, obtain the height value of grid and a plurality of grids around this grid at described sampled point place;
Described elevation rate of change is according to isocontour spacing and the difference of elevation acquisition nearest with described sampled point;
Described the first height value acquisition module specifically comprises:
Grid obtains submodule, is used for the grid according to described sampled point place, determines described grid a plurality of grids of vicinity on every side;
Distance is obtained judge module, is used for obtaining described sampled point and the described distance of each grid element center point of contiguous a plurality of grids on every side; And whether judge described distance greater than threshold value, and, if less than threshold value, change height value over to and obtain submodule one,, if greater than threshold value, change the elevation variable quantity over to and obtain the judgement submodule;
Height value is obtained submodule one, with the height value of the height value of this sampled point as this grid;
The elevation variable quantity obtains the judgement submodule, be used for the elevation rate of change according to described distance and corresponding sampled point, obtain the elevation variable quantity of corresponding described grid, and judge that described elevation variation increases progressively or successively decreases, if increase progressively, change height value over to and obtain submodule two,, if successively decrease, change height value over to and obtain submodule three;
Described elevation variable quantity obtains the judgement submodule and specifically comprises:
Judgement submodule one, be used for determining that the central spot of described grid is inner or outside in this sampled point place level line;
Judgement submodule two, be used for increasing progressively or successively decreasing according to the described central point of described contour map distribution trend judgement;
Height value is obtained submodule two, is used for the height value of described corresponding sampled point is added result that the elevation variable quantity of corresponding described grid the obtains height value as described grid;
Height value is obtained submodule three, is used for the height value of described corresponding sampled point is deducted result that the elevation variable quantity of corresponding described grid the obtains height value as described grid;
The second height value acquisition module, be used for traveling through whole grid, to being in the grid of the state of being untreated, search the height value of the grid of having processed of the multiple directions the most contiguous with described grid, get the average of value of described a plurality of grids of having processed as the height value of the described grid that is untreated;
The digital elevation model output module, be used for the whole digital elevation model of height value output according to whole grid; Described the second height value acquisition module also comprises:
The first submodule, determine that the central spot of described grid is inner or outside in this sampled point place level line;
The second submodule, increase progressively or successively decrease according to the described central point of described contour map distribution trend judgement; Specifically comprise: when contour map raise from innermost layer to outermost height, if described central point is inner at described high line, the elevation variation was successively decreased;
When contour map reduced from innermost layer to outermost height, if described central point is inner at described high line, the elevation variation increased progressively;
Overlapping processing module, be used for when this grid is untreated, with the height value of this grid of acquisition directly as the height value of this grid; When this grid has been processed, height value and the existing height value of this grid of this grid of acquisition are got the height value of average as this grid.
Compared with prior art, the application comprises following advantage:
the application passes through at first to the elevation rate of change of each known sampled point and the height value of this sampled point, obtain described 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 the height value of the grid of having processed of the multiple directions the most contiguous with described grid, get the average of value of described a plurality of grids of having processed as the height value of the described grid that is untreated, the process that obtains so whole grid height value does not need to select with reference to the process of point set directly to generate DEM, and can be in the situation that guarantee the good quality of data, computing velocity efficiency is significantly improved.
Description of drawings
Fig. 1 is the schematic flow sheet of a kind of method from the contour map generating digital elevation model of the application;
Fig. 2 is the actual level line legend of of the application;
Fig. 3 is the schematic flow sheet of preferred the first height value obtaining step of a kind of method from the contour map generating digital elevation model of the application;
Fig. 4 is the preferred schematic flow sheet of further processing a grid of a kind of method from the contour map generating digital elevation model of the application;
Fig. 5 is the digital elevation model figure that the application is generated by actual level line legend;
Fig. 6 is the structural representation of a kind of device from the contour map generating digital elevation model of the application;
Fig. 7 is the structural representation of preferred the first height value acquisition module of a kind of device from the contour map generating digital elevation model of the application.
Embodiment
, for above-mentioned purpose, the feature and advantage that make the application can become apparent more, below in conjunction with the drawings and specific embodiments, the application is described in further detail.
the application passes through at first to the elevation rate of change of each known sampled point and the height value of this sampled point, obtain described 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 the height value of the grid of having processed of the multiple directions the most contiguous with described grid, get the average of value of described a plurality of grids of having processed as the height value of the described grid that is untreated, the process that obtains so whole grid height value does not need to select with reference to the process of point set directly to generate DEM, and can be in the situation that guarantee the good quality of data, computing velocity efficiency is significantly improved.
With reference to Fig. 1, show the schematic flow sheet of a kind of method from the contour map generating digital elevation model of the application, its step comprises:
Grid generates step 110, and the vertical and horizontal resolution of digital elevation model as required, be divided into grid with the quadrilateral bounding box of contour map.
After obtaining contour map, general at first by calculating the bounding box of contour map, obtain comprising the minimum quadrilateral of this contour line chart, and then the resolution grid division number vertical, horizontal according to DEM, wherein when grid division take minimum quadrilateral as basis, progressively adjust the quadrilateral of the encirclement contour map meet an integer grid.After the grid division number, can these trellis state be set to " being untreated " in this step.
With reference to Fig. 2, it is the application's an actual level line legend, and this contour line chart is the contour map of a hillside type, and its isocontour value is successively decreased to outermost layer gradually by innermost layer, wherein the level line difference of elevation is 100 meters, and the isocontour height value of innermost layer is 1000 meters.
comprise all isocontour minimum quadrilaterals by calculating, it is also the frame rectangle in Fig. 2, obtain simultaneously long a and the wide b of this rectangle, the demand of DEM resolution according to actual needs, the resolution of DEM is defined as 50 meters, the sizing grid that namely obtains DEM is a/50, b/50, wherein, the result of this calculating is just not necessarily integer, words such as a=110, but number of grid is integer certainly, so, get the smallest positive integral 2 again+1 of a/50=2.2, be that we get 3, readjust so the long a value of rectangle, a=3*50=150 rice now, carry out aforesaid operations for b/50 with same principle and adjust the b value, finally make a/50 and b/50 be integer.Here the original state of each grid is " being untreated ".
The first height value obtaining step 120, for each sampled point on each level line in described contour map,, according to the elevation rate of change of described sampled point and the height value of described sampled point, obtain the height value of grid and a plurality of grids around this grid at described sampled point place.
Have known a plurality of sampled points in reality on each level line of contour map.To level line CL all in contour map
1, CL
2..., CL
nCarry out circular treatment, for level line CL
i(n 〉=i 〉=1), to all the sampled point P on it
1, P
2..., P
mCircular treatment.Wherein, P no matter
jWhere (m 〉=i 〉=1) is, it must be in a certain grid of DEM or on limit, centered by this grid, it has a plurality of grids on every side, preferably, get centered by a grid adjacent 8 grids (totally 8 of upper and lower, left and right, upper left, lower-left, upper right, bottom rights) around it in the application, can obtain the macrolattice A of 3 * 3, then according to the height value of the elevation rate of change of this sampled point and described sampled point, calculate the height value of 9 grids in macrolattice A the inside.
Further,, with reference to Fig. 3, be the schematic flow sheet of preferred the first height value obtaining step of a kind of method from the contour map generating digital elevation model of the application, specifically comprise:
For the sampled point (P on certain level line
x, P
y), find that level line nearest with this sampled point on two level lines adjacent with this sampled point place level line, and calculate described minimum distance L, in conjunction with the difference of elevation d between level line, obtain the elevation rate of change in this shop, generally by following formula, obtain elevation rate of change: T
j=d/L, wherein T
jThe elevation rate of change that represents this sampled point.
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 can be on a grid inside or limit, when a grid was inner, directly grid was chosen f grid on every side centered by this grid when it; When it on the limit of a grid time, choose at random one of them grid of the grid of common limit or concurrent and choose grid around it as central gridding; The application preferably gets adjacent 8 grids on every side of central gridding, can obtain the macrolattice A of 3 * 3.Such as in Fig. 2, have 35 sampled points in the level line of innermost layer,, to each sampled point P, by its locus P (x, y), calculated at the position of DEM grid (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)。
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 described grid to this sampled point apart from Q, whether judge it greater than threshold value,, if less than threshold value, change step 124 over to, if greater than threshold value 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) be Q to the distance of sampled point P (x, y), but this Q value is less than the length of a grid, the time, change so step 124 over to; But (P
x, P
y+1) to sampled point apart from the length of Q value greater than a grid, change step 125 over to.
Wherein, generally, the central gridding at sampled point place is directly used the height value of sampled point as the height value of this central gridding.
Because this grid is to the grid of sampled point distance less than threshold value, so directly with the height value of the height value of this sampled point as this grid.Such as the height value of a sampled point P (x, y) on aforementioned innermost layer level line is 1000, and (P
x+1, P
y+1) be length less than a grid to the distance of sampled point P (x, y), with 1000 directly as grid (P
x+1, P
y+1) height value.
, for to the grid of sampled point distance greater than threshold value,, according to the elevation rate of change of described distance and corresponding sampled point, obtain the elevation variable quantity of corresponding described grid, and judge that described elevation variable quantity increases progressively or successively decreases, and, if increase progressively, changes step 126 over to,, if successively decrease, change step 127 over to.Wherein the elevation variable quantity can calculate in the following way: W=Q * T
j, wherein W represents elevation variable quantity that should grid, T
jTo elevation rate of change that should sampled point, Q is the distance that aforesaid this grid element center is put this sampled point for aforesaid expression.Such as aforementioned (P
x, P
y+1) to sampled point apart from the length of Q value greater than a grid, calculate elevation variable quantity W that should grid, and judge that described elevation variation increases progressively or successively decreases.
Further, can judge that the elevation variable quantity increases progressively or the process of successively decreasing by following steps:
Steps A 1, determine that the central spot of described grid is inner or outside in this sampled point place level line;
Steps A 2, increase progressively or successively decrease according to the described central point of described contour map distribution trend judgement.
Further, described steps A 2 specifically comprises:
When contour map raise from innermost layer to outermost height, if described central point is inner at described high line, the elevation variation was successively decreased;
When contour map reduced from innermost layer to outermost height, if described central point is inner at described high line, the elevation variation increased progressively.
This process can be to increase progressively limit or successively decrease limit at this isocontour elevation according to this grid; Isocontour elevation increases progressively limit definition: if level line internal point elevation greater than on level line, the level line internal point belongs to elevation and increases progressively limit; If level line internal point elevation is less than on level line, the level line external point belongs to elevation and increases progressively limit.
Such as the elevation variable quantity of a grid is W, the height value of this sampled point is S, and the elevation that this grid is corresponding changes and increases progressively, the height value Z=S+W of this grid.
Can obtain like this height value of the central gridding at a sampled point place and grid on every side thereof.
What so circulate 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 processed, carry out following steps:
Judge whether this grid is processed: when this grid is untreated, with the height value of this grid of obtaining directly as the height value of this grid; When this grid has been processed, height value and the existing height value of this grid of this grid of acquisition 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 processes for the first time is g, and sampled point P
3Also comprise grid d in corresponding grid, then by abovementioned steps, obtained P
3Obtain after should another height value h of grid, with the mean value of g and h namely (g+h)/2 as the new height value of grid d.
Third high journey value obtaining step 130, travel through whole grid, to being in the grid of the state of being untreated, search the height value of the grid of having processed of the multiple directions the most contiguous with described grid, get the average of value of described a plurality of grids of having processed as the height value of the described grid that is untreated.
In reality, after a plurality of grids corresponding to sampled point are processed, 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 processing of multiple 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 that this grid closes on most, and (maximum 8 are respectively upper and lower, left and right, upper left, upper right, lower-left, bottom right., if on border, get in esse several direction) the grid of having processed, ask the mean value of height value of these grids as the height value of this grid.
Digital elevation model output step 140, according to the whole digital elevation model of height value output of the whole grid that obtains.
, with reference to Fig. 5, be the final DEM that is obtained through aforementioned process by contour map in Fig. 2.Represent that as the figure innermost layer height value is larger, belong to the summit, hillside, outermost layer represents the subregion that elevation is lower.
The application's method has been removed and has been selected directly to generate with reference to the process of point set DEM, and can be in the situation that guarantee the good quality of data, and computing velocity efficiency is significantly improved.
Disclosed herein as well is accordingly a kind of device from the contour map generating digital elevation model.
With reference to Fig. 6, be the structural representation of a kind of device from the contour map generating digital elevation model of the application, described device comprises:
The grid generation module, digital elevation model vertical and horizontal resolution as required, be divided into grid with the quadrilateral bounding box of contour map.
The first height value acquisition module, be used for for each sampled point on each level line of described contour map,, according to the elevation rate of change of described sampled point and the height value of described sampled point, obtain the height value of grid and a plurality of grids around this grid at described sampled point place.
The second height value acquisition module, be used for traveling through whole grid, to being in the grid of the state of being untreated, search the height value of the grid of having processed of the multiple directions the most contiguous with described grid, get the average of value of described a plurality of grids of having processed as the height value of the described grid that is untreated.
The digital elevation model output module, be used for the whole digital elevation model of height value output according to whole grid.
Further,, with reference to Fig. 7, be the structural representation of preferred the first height value acquisition module of a kind of device from the contour map generating digital elevation model of the application.Described the first height value acquisition module specifically comprises:
The elevation rate of change obtains submodule, is used for for the sampled point on a level line of described contour map,, according to isocontour spacing and the difference of elevation nearest with described sampled point, obtains the elevation rate of change of described sampled point.
Grid obtains submodule, is used for the grid according to described sampled point place, determines described grid a plurality of grids of vicinity on every side.
Distance is obtained judge module, is used for obtaining described sampled point and the described distance of each grid element center point of contiguous a plurality of grids on every side; And whether judge described distance greater than threshold value, and, if less than threshold value, change height value over to and obtain submodule one,, if greater than threshold value, change the elevation variable quantity over to and obtain the judgement submodule.
Height value is obtained submodule one, with the height value of the height value of this sampled point as this grid;
The elevation variable quantity obtains the judgement submodule, be used for the elevation rate of change according to described distance and corresponding sampled point, obtain the elevation variable quantity of corresponding described grid, and judge that described elevation variation increases progressively or successively decreases, if increase progressively, change height value over to and obtain submodule two,, if successively decrease, change height value over to and obtain submodule three.
Further, the elevation variable quantity obtains and judges that submodule comprises:
Judgement submodule one, be used for determining that the central spot of described grid is inner or outside in this sampled point place level line.
Judgement submodule two, be used for increasing progressively or successively decreasing according to the described central point of described contour map distribution trend judgement.
Height value is obtained submodule two, is used for the height value of described corresponding sampled point is added result that the elevation variable quantity of corresponding described grid the obtains height value as described grid.
Height value is obtained submodule three, is used for the height value of described corresponding sampled point is deducted result that the elevation variable quantity of corresponding described grid the obtains height value as described grid.
Further, also comprise overlapping processing module, be used for when this grid is untreated, with the height value of this grid of acquisition directly as the height value of this grid; When this grid has been processed, height value and the existing height value of this grid of this grid of acquisition are got the height value of average as this grid.
Further, described the second height value acquisition module also comprises:
The first submodule, determine that the central spot of described grid is inner or outside in this sampled point place level line.
The second submodule, increase progressively or successively decrease according to the described central point of described contour map distribution trend judgement.
For system embodiment, because it is substantially similar to embodiment of the method, so description is fairly simple, relevant part gets final product referring to the part explanation of embodiment of the method.
Each embodiment in this instructions all adopts the mode of going forward one by one to describe, and what each embodiment stressed is and the difference of other embodiment that between each embodiment, identical similar part is mutually referring to getting final product.
Above a kind of method and apparatus from the contour map generating digital elevation model that the application is provided, be described in detail, applied specific case herein the application's principle and embodiment are set forth, the explanation of above embodiment just is used for helping to understand the application's method and core concept thereof; Simultaneously, for one of ordinary skill in the art, the thought according to the application, all will change in specific embodiments and applications, and in sum, this description should not be construed as the restriction to the application.
Claims (3)
1. the method from the contour map generating digital elevation model, is characterized in that, comprising:
Grid generates step, and the vertical and horizontal resolution of digital elevation model as required, be divided into grid with the quadrilateral bounding box of contour map;
The first height value obtaining step, for each sampled point on each level line in described contour map,, according to the elevation rate of change of described sampled point and the height value of described sampled point, obtain the height value of grid and a plurality of grids around this grid at described sampled point place;
The second height value obtaining step, travel through whole grid, to being in the grid of the state of being untreated, search the height value of the grid of having processed of the multiple directions the most contiguous with described grid, get the average of value of described a plurality of grids of having processed as the height value of the described grid that is untreated;
Digital elevation model output step, according to the whole digital elevation model of height value output of the whole grid that obtains;
Wherein, described the first height value obtaining step further comprises:
Step 221,, for the sampled point on the level line in described contour map,, according to isocontour spacing and the difference of elevation nearest with described sampled point, obtain the elevation rate of change of described sampled point;
Step 222,, according to the grid at described sampled point place, determine described grid a plurality of grids of vicinity on every side;
Step 223, obtain described sampled point and described around the distance of each grid element center point in contiguous a plurality of grids; And whether judge described distance greater than threshold value,, if less than threshold value, change step 224 over to, if, greater than threshold value, change step 225 over to;
Step 224, with the height value of the height value of this sampled point as this grid;
Step 225,, according to the elevation rate of change of described distance and corresponding sampled point, obtain the elevation variable quantity of corresponding described grid, and judge that described elevation variation increases progressively or successively decreases, and, if increase progressively, changes step 226 over to,, if successively decrease, changes step 227 over to;
Step 226, add result that the elevation variable quantity of corresponding described grid the obtains height value as described grid with the height value of described corresponding sampled point;
Step 227, deduct result that the elevation variable quantity of corresponding described grid the obtains height value as described grid with the height value of described corresponding sampled point;
Wherein, described the second height value obtaining step further comprises:
Judge whether this grid is processed: when this grid is untreated, with the height value of this grid of obtaining directly as the height value of this grid; When this grid has been processed, height value and the existing height value of this grid of this grid of acquisition are got the height value of average as this grid;
Wherein, judge in described step 225 that the elevation variation increases progressively or the process of successively decreasing further comprises:
Steps A 1, determine that the central spot of described grid is inner or outside in this sampled point place level line;
Steps A 2, increase progressively or successively decrease according to the described central point of described contour map distribution trend judgement;
Wherein, described steps A 2 criterions are:
When contour map raise from innermost layer to outermost height, if described central point is inner at described level line, the elevation variation was successively decreased;
When contour map reduced from innermost layer to outermost height, if described central point is inner at described level line, the elevation variation increased progressively.
2. method according to claim 1 is characterized in that:
Step 123 further comprises: for sampled point place grid, directly with the height value of the height value of sampled point as this grid.
3. the device from the contour map generating digital elevation model, is characterized in that, comprising:
The grid generation module, be used for digital elevation model vertical and horizontal resolution as required, and the minimum quadrilateral partition of contour map is become grid;
The first height value acquisition module, be used for for each sampled point on each level line of described contour map,, according to the elevation rate of change of described sampled point and the height value of described sampled point, obtain the height value of grid and a plurality of grids around this grid at described sampled point place;
The second height value acquisition module, be used for traveling through whole grid, to being in the grid of the state of being untreated, search the height value of the grid of having processed of the multiple directions the most contiguous with described grid, get the average of value of described a plurality of grids of having processed as the height value of the described grid that is untreated;
The digital elevation model output module, be used for the whole digital elevation model of height value output according to whole grid;
Wherein, described the first height value acquisition module further comprises:
Module 221, be used for for the sampled point on a level line of described contour map,, according to isocontour spacing and the difference of elevation nearest with described sampled point, obtains the elevation rate of change of described sampled point;
Module 222, be used for the grid according to described sampled point place, determines described grid a plurality of grids of vicinity on every side;
Module 223, be used for obtaining described sampled point and the described distance of each grid element center point of contiguous a plurality of grids on every side; And whether judge described distance greater than threshold value, and, if less than threshold value, change module 224 over to,, if greater than threshold value, change module 225 over to;
Module 224, be used for the height value of the height value of this sampled point as this grid;
Module 225, be used for the elevation rate of change according to described distance and corresponding sampled point, obtains the elevation variable quantity of corresponding described grid, and judge that described elevation variation increases progressively or successively decreases, and, if increase progressively, changes module 226 over to,, if successively decrease, changes module 227 over to;
Module 226, be used for the height value of described corresponding sampled point is added result that the elevation variable quantity of corresponding described grid the obtains height value as described grid;
Module 227, be used for the height value of described corresponding sampled point is deducted result that the elevation variable quantity of corresponding described grid the obtains height value as described grid;
Wherein, described the second height value acquisition module further comprises:
Be used for judge the module whether this grid is processed, when this grid is untreated, with the height value of this grid of acquisition directly as the height value of this grid; When this grid has been processed, height value and the existing height value of this grid of this grid of acquisition are got the height value of average as this grid;
Wherein, described module 225 is carried out the described elevations of judgement and is changed and increase progressively or further comprise while successively decreasing:
Modules A 1, be used for determining that the central spot of described grid is inner or outside in this sampled point place level line;
Modules A 2, be used for increasing progressively or successively decreasing according to the described central point of described contour map distribution trend judgement;
Wherein, described modules A 2 criterions are:
When contour map raise from innermost layer to outermost height, if described central point is inner at described level line, the elevation variation was successively decreased;
When contour map reduced from innermost layer to outermost height, if described central point is inner at described level line, the elevation variation increased progressively.
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 CN102339478A (en) | 2012-02-01 |
| CN102339478B true 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) |
Families Citing this family (8)
| 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 |
| CN104331928B (en) * | 2014-10-28 | 2017-10-03 | 中国航空工业集团公司洛阳电光设备研究所 | A kind of contour automatic labeling method based on the triangulation network |
| CN105046267B (en) * | 2015-05-19 | 2018-08-28 | 广州杰赛科技股份有限公司 | Generation method, generating means, polymerization and the polyplant of dynamic grid |
| CN107393002B (en) | 2017-07-28 | 2019-02-01 | 北京金风科创风电设备有限公司 | Apparatus and method for extracting terrain boundary |
| CN108280880B (en) * | 2018-01-24 | 2018-11-16 | 长春工程学院 | A method of improving the digital elevation data resolution of massif using remote sensing image |
| CN108898681B (en) * | 2018-06-26 | 2022-07-05 | 中煤航测遥感集团有限公司 | Digital elevation model processing method and device |
| CN109191579B (en) * | 2018-07-02 | 2019-06-28 | 广东省地质测绘院 | Method, apparatus, computer equipment and the storage medium of coordinate conversion |
| CN112419495B (en) * | 2020-10-26 | 2022-11-15 | 天津大学 | Elevation point automatic extraction method based on multi-scale DEM space model |
Citations (1)
| 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 |
-
2011
- 2011-09-14 CN CN2011102724447A patent/CN102339478B/en active Active
Patent Citations (1)
| 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 (2)
| Title |
|---|
| GIS数字地形模型的建立及应用研究;傅炜;《陕西师范大学学报(自然科学版)》;20010630;第29卷(第02期);第92-97页 * |
| 傅炜.GIS数字地形模型的建立及应用研究.《陕西师范大学学报(自然科学版)》.2001,第29卷(第02期),第92-97页. |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102339478A (en) | 2012-02-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102339478B (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 | |
| CN104821013A (en) | Method and system for specific surface area extraction based on geodetic coordinate system digital elevation model | |
| CN104751733B (en) | The region method for drafting and device of map, path distance sorting technique and system | |
| CN108241369B (en) | Method and device for avoiding static obstacle for robot | |
| CN104077134B (en) | Obtain the method and system of data center machine room face temperature cloud atlas | |
| CN108629231A (en) | Obstacle detection method, device, equipment and storage medium | |
| WO2010106018A1 (en) | 3d building generalization for digital map applications | |
| CN102496187B (en) | Method for tracking contour line to boundary and fault based on triangular mesh | |
| CN107220372B (en) | A kind of automatic laying method of three-dimensional map line feature annotation | |
| CN108801273A (en) | A kind of generation method and device of road reference line | |
| CN102270236A (en) | Rasterized geographic information system (GIS)-based spatial relationship judging method and system | |
| CN113286982A (en) | System and method for generating, updating and enhancing large-scale high-precision 3D road map and multi-level road map | |
| US20110050689A1 (en) | Method and Apparatus for Creating of 3D Direction Displaying | |
| CN107564098A (en) | A kind of Fast rendering method of big Local Area Network three-dimensional noise map | |
| CN112233205B (en) | Electronic map making method and system for partitioning and cutting based on discrete data | |
| CN105844224A (en) | Point cloud fast ordering method for on-vehicle LiDAR road points | |
| CN111047682A (en) | Three-dimensional lane model generation method and system | |
| CN106383828A (en) | Drawing method and electronic equipment | |
| CN106327576B (en) | A kind of City scenarios method for reconstructing and system | |
| CN106408531A (en) | GPU acceleration-based hierarchical adaptive three-dimensional reconstruction method | |
| CN106156281A (en) | Arest neighbors point set method for quickly retrieving based on Hash Cube spatial level partition structure | |
| CN110096993A (en) | The object detection apparatus and method of binocular stereo vision | |
| CN105893590A (en) | Automatic processing method for real-situation cases of DTA (Digital Terrain Analysis) modelling knowledge | |
| CN105023288B (en) | Collimation error removing method of the two-dimensional vector solid line in three-dimensional scenic |
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 |