CN105427380B - One kind is based on detail three-dimensional map data processing method - Google Patents

One kind is based on detail three-dimensional map data processing method Download PDF

Info

Publication number
CN105427380B
CN105427380B CN201510979549.4A CN201510979549A CN105427380B CN 105427380 B CN105427380 B CN 105427380B CN 201510979549 A CN201510979549 A CN 201510979549A CN 105427380 B CN105427380 B CN 105427380B
Authority
CN
China
Prior art keywords
dem
math
data
model
tile
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
Application number
CN201510979549.4A
Other languages
Chinese (zh)
Other versions
CN105427380A (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.)
Kunming Enersun Technology Co Ltd
Original Assignee
Kunming Enersun Technology 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 Kunming Enersun Technology Co Ltd filed Critical Kunming Enersun Technology Co Ltd
Priority to CN201510979549.4A priority Critical patent/CN105427380B/en
Publication of CN105427380A publication Critical patent/CN105427380A/en
Application granted granted Critical
Publication of CN105427380B publication Critical patent/CN105427380B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05—Geographic models
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/29—Geographical information databases
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20112—Image segmentation details
    • G06T2207/20128—Atlas-based segmentation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • Remote Sensing (AREA)
  • Computer Graphics (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Processing Or Creating Images (AREA)
  • Instructional Devices (AREA)

Abstract

One kind is based on detail three-dimensional map data processing method, including 3 steps：Digital orthoimage cutting, the cutting of digital elevation data, altitude data is inlayed and edge fit.Have three-dimensional map precision at present and be up to 20 meters, when view distance surface texture is nearer, single texture pixel can cover multiple pixels under the action of texture amplifier, this just makes image become very fuzzy.One kind provided by the invention is based on detail three-dimensional map data processing method, improves clarity of the three-dimensional map in displaying, landforms situation of the relief detail closer to reality；The real time access performance of three-dimensional map is improved by the database storage techniques based on threedimensional model, database server pressure when accessing three-dimensional map is reduced and improves concurrent access performance at the same time.

Description

One kind is based on detail three-dimensional map data processing method
Technical field
The invention belongs to power grid three-dimensional map data processing method.
Background technology
With the continuous development of GIS technology, three-dimension GIS has been more and more widely used in whole power industry, base In the three-dimensional simulation of real scene data, shown in various aspects such as administration of power networks, breakdown repair, security monitorings outstanding Effect., will digitally face mould using technologies such as advanced GIS, RS and virtual realities by superior image, data qualification Type, power transmission and transforming equipment model and various power department professional attributes information combine, and establish electric power three-dimensional geographic information Platform, it can be achieved that the power specialty data being combined with Fundamental Geographic Information Data inquiry, renewal, the maintenance of mains-power circuit and Safe prosecution, realizes space performance and analysis, the management function of power grid in large scene.Platform passes through advanced three-dimensional visualization hand Section, includes computer management, specification power transmission and transformation operation flow, strengthens power department by whole power transmission and transformation business and management overall process Cooperation and management function, improve geographical department's operation of transmission lines and substations, managerial ability and level of decision-making.
Have three-dimensional map precision at present and be up to 20 meters, when view distance surface texture is nearer, single line Pixel is managed under the action of texture amplifier, multiple pixels can be covered, this just makes image become very fuzzy.
On the other hand, high resolution image includes more surface details really, but these details are all often unusual phases As.Details as offer will undoubtedly consume substantial amounts of texture memory, although this can give decision height and action to provide line Rope, avoids the fuzzy of surface, but causes texture memory inadequate, so that many actually useful information can not in time, smoothly Ground loads.
The content of the invention
The present invention is precisely in order to the one kind for overcoming drawbacks described above and designing is handled based on detail three-dimensional map data Method, the present invention are achieved by the following technical solution：
One kind is based on detail three-dimensional map data processing method, the invention is characterised in that, including 3 steps：Number The cutting of word orthography, the cutting of digital elevation data, altitude data is inlayed and edge fit；Wherein：
Digital orthoimage is cut：Digital orthoimage cutting is to cut into satellite image picture according to pyramid model The tile of 256*256 pixels, the tile cut carry out coding name according to the form of l_r_c.png, and l is pyramid model Level number, r are the line number of pyramid model, and c is the row number of pyramid model；
Digital elevation data are cut：The cutting of digital elevation data is that contour and elevational point are established irregular triangle network, Then the small cubes for cutting into 256*256*256 pixels according to pyramid model on the basis of irregular establish dem data；Cut The small cubes cut establish index according to level, row, column, block, and with the form storage of OBJ to database；
Altitude data edge fit and inlay：Based on the small cubes of digital elevation data cutting, choose consecutive phantom and generated Dem data, the elevation of same level coordinate grid site is checked, if there are the poor lattice for being more than 2 times of DEM mean square error of height of elevation Site, then being considered as transfinites, and is regarded as rough error point, and re-establish body Model；Edge fit is carried out to there is the dem data of rough error point Repair edge fit again after survey；It is sequentially completed by the above process between all single model dem datas in succession；If institute in fixed range There is the edge fit of single model dem data is poor all to meet regulation requirement, then can use the tile that digital orthoimage is cut to DEM Inlayed；The elevation of all same level coordinate net lattice points when inlaying to participating in edge fit takes its average value, as respective lattice The height value of site.
One kind is based on detail three-dimensional map data processing method, the invention is characterised in that, using following technology hand Section：
1) tile pyramid structure is set up
Tile map pyramid model is a kind of multiresolution hierarchical model, from the pyramidal bottom of tile to top layer, point Resolution is more and more lower, but the geographic range represented is constant, mainly includes two parts of tile and pyramid；Wherein：
1.1) tile
The tile be exactly according to a certain percentage size piece image is divided into some small square grids, each grid Lattice are exactly a tile；
Map tile after section is grating image, does not have location information, but related slices calculation has been used in section After method, the position specifically positioned can be calculated.
1.2) pyramid model
According to the demand of user, piece image is divided into ascending region.After segmentation, image forms ratio by small To big, the pyramid structure of quantity from less to more：Fig. 1 is pyramid planar structure, and Fig. 2 is intuitively pyramid knot Structure.
Tile is small square block, and length and width are equal, and length and width length is 256 pixels, 8~18 grades of totally 11 grades of engineer's scales；Ratio Ruler rank increases level-one；After having cut coding name is carried out according to the form of " l_r_c.png "；Wherein l is engineer's scale rank, R is line number, and c is row number, 0<=r<=2aa, 0<=c<=2aa；Coding rule is cut according to above-mentioned picture, obtains picture row Mutual conversion between row coordinate and actual geographic coordinate, these are the bases of map denotation.
Longitude and latitude is converted to line number and row number under some proportion ruler, and transfer process is：
N=Math.pow (2, lev)
X=((Longitude+180)/360) * n
SinLat=Math.Sin (Math.PI*Latitude/180)；
Y=0.5-Math.Log ((1+sinLat)/(1-sinLat))/(4*Math.PI)；
Wherein, Longitude is longitude, and Latitude is latitude value, and lev is level of zoom, and Math.pow is calculating 2 Lev power, Math.PI is pi, and Math.Sin is SIN function, and Math.Log is takes the logarithm.
Line number and row number under some proportion ruler are converted into latitude, and transfer process is：
N=Math.pow (2, lev)
Longitude=x/n*360-180；
Latitude=Math.Atan (Math.Sinh (Math.PI* (1-2*y/n)))；
Latitude=Latitude*180.0/Math.PI；
Wherein, Longitude is longitude, and Latitude is latitude value, and lev is level of zoom, and Math.pow is calculating 2 Lev power, Math.PI is pi, and Math.Atan is arctan function, and Math.Sinh is hyperbolic sine function.
2) DEM terrain datas are cut
Digital elevation model, vehicle economy M；It is a kind of entity that ground elevation is represented with one group of orderly array of values form Ground model, is a branch of digital terrain model, thus other various topographic index can derive from；It is generally believed that DTM It is various geomorphologic factors of the description including elevation, it is linear and non-thread including such as gradient, slope aspect, the change of slope factor Property combination spatial distribution, wherein DEM is the simple individual event digital land value model model of zeroth order, other as the gradient, slope aspect and the gradient become The landforms such as rate characteristic can derive from the basis of DEM.
Irregular triangle network, abbreviation TIN；The reconstruction of topographical surface is directly carried out using original sample point, by continuous phase The triangular facet composition being mutually coupled with, the shapes and sizes of triangular facet depend on density and the position of the observation station of irregular distribution.
The cutting of digital elevation data is that contour and elevational point are established irregular triangle network, then on the basis of irregular The small cubes that 256*256*256 pixels are cut into according to the layer of pyramid model, row, column establish dem data；Can be with different layers Secondary resolution ratio describes surface configuration；Can more accurately it be represented more with less room and time under a certain specified resolution Add the surface of complexity；Especially when landform includes big measure feature such as geosutures, construction line, these features can be preferably taken into account； The small cubes cut establish joint index according to level, row, column, block, and with the form storage of OBJ to database.
Platform directly parses obj file, obtains the information such as model vertices, textures and is stored by specification format to data Preserved in storehouse, the key element of parsing is as shown in Figure 3.
3) altitude data edge fit and inlay
3.1) DEM is edited
The DEM grid points that dem data editor refers to form interpolation are one by one into edlin；Identical terrain category DEM grid Point edge fit limit difference is twice of error in the terrain category grid points, and different terrain classification DEM edge fits limit difference is two kinds of landform classes The sum of other DEM grid points edge fit limit difference.
3.2) dem data edge fit
Choose the dem data that is generated of consecutive phantom, check the elevation of same level coordinate grid site, if occur elevation compared with Difference is more than the grid points of 2 times of DEM mean square error of height, then being considered as transfinites, and is regarded as rough error point, and re-establish body Model；It is right The dem data for rough error point occur carries out edge fit and repaiies edge fit again after survey；All single model DEM numbers are sequentially completed by the above process Between in succession.
3.3) dem data is inlayed and cut
If the edge fit of all single model dem datas is poor in fixed range all meets regulation requirement, numeral can be used The tile of orthography cutting inlays DEM；Height when inlaying to all same level coordinate net lattice points for participating in edge fit Journey takes its average value, as the height value of respective grid points, while forms the match accuracy report on each bar side.
After the completion of DEM inlays, carry out rectangle according to start-stop grid point coordinates as defined in relevant regulations or technical requirements and cut When, a row or multi-row DEM grid can be extended out according to specific requirement.
4) DEM quality controls
The quality control of DEM includes production process quality control and end result quality control two parts.
4.1) production process quality control
Whether the correctness and data that the key points of quality control inspection firsthand information in production process uses are reasonable.
Whether lake, reservoir, the layering in two-wire river are reasonable, and whether the elevation assignment in waterside line and coastline is rationally correct, DEM height values should be consistent in the range of static water body, flows the DEM height values in waters and gentle transition, relation should close from top to bottom Reason.
DEM elevations supposition area should be divided into by not reaching the region of regulation height accuracy requirement.
Grid positioned at white space should assign height value -9999, and the processing to clear area will be fully recorded in first number In.
4.2) end result quality control
The quality control of end result visually inspects whether contour has mutation feelings by the method for DEM interpolated contours Condition, or compared with topographic map, when landform shape, same place elevation differ greatly, then judge there are quality problems and repaiied Change.
Check the correctness of dem data start-stop point coordinates, check whether height value effective range is correct.
It should be checked after DEM splicings, determine whether overlapping and crack, whether splicing precision reaches requirement.
Usually used DEM files are gray-scale maps, its processing procedure is as shown in Figure 4.
The coordinate information of model is stored in the Double arrays that a structure is fixed, can be first during model is shown Coordinate information is read, to judge whether the model shows, when display goes the data of database reading model again, is reduced with this The interaction times and data volume of network.Specific coordinate information storage organization is as shown in Figure 5.
Significant data of the striograph as three-dimensional scenic, is stored in oracle database, it is put in storage process such as Fig. 6 institutes Show.Fig. 7 is put in storage process for model.
The present invention key point be：
1st, the three-dimensional map data treatment technology of detail.
2nd, the three-dimensional map data library storage technology based on ObjModel.
Compared with prior art, the present invention has the advantage that or beneficial effect is：
1st, clarity of the three-dimensional map in displaying, landforms situation of the relief detail closer to reality are improved；
2nd, the real time access performance of three-dimensional map is improved by the database storage techniques based on threedimensional model, reduces and visits Ask that database server pressure during three-dimensional map improves concurrent access performance at the same time.
Brief description of the drawings
Fig. 1 is pyramid planar structure；
Fig. 2 is intuitively pyramid structure；
Fig. 3 wants sketch map for what platform directly parsed obj file；
Fig. 4 is the processing procedure of DEM files；
Fig. 5 is coordinate information storage organization figure；
Fig. 6 is put in storage process for flat image figure；
Fig. 7 is put in storage process for model.
Embodiment
See Fig. 1-Fig. 7, one kind is based on detail three-dimensional map data processing method, the invention is characterised in that, including 3 A step：Digital orthoimage cutting, the cutting of digital elevation data, altitude data is inlayed and edge fit；Wherein：
Digital orthoimage is cut：Digital orthoimage cutting is to cut into satellite image picture according to pyramid model The tile of 256*256 pixels, the tile cut carry out coding name according to the form of l_r_c.png, and l is pyramid model Level number, r are the line number of pyramid model, and c is the row number of pyramid model；
Digital elevation data are cut：The cutting of digital elevation data is that contour and elevational point are established irregular triangle network, Then the small cubes for cutting into 256*256*256 pixels according to pyramid model on the basis of irregular establish dem data；Cut The small cubes cut establish index according to level, row, column, block, and with the form storage of OBJ to database；
Altitude data edge fit and inlay：Based on the small cubes of digital elevation data cutting, choose consecutive phantom and generated Dem data, the elevation of same level coordinate grid site is checked, if there are the poor lattice for being more than 2 times of DEM mean square error of height of elevation Site, then being considered as transfinites, and is regarded as rough error point, and re-establish body Model；Edge fit is carried out to there is the dem data of rough error point Repair edge fit again after survey；It is sequentially completed by the above process between all single model dem datas in succession；If institute in fixed range There is the edge fit of single model dem data is poor all to meet regulation requirement, then can use the tile that digital orthoimage is cut to DEM Inlayed；The elevation of all same level coordinate net lattice points when inlaying to participating in edge fit takes its average value, as respective lattice The height value of site.
One kind is based on detail three-dimensional map data processing method, the invention is characterised in that, using following technology hand Section：
1) tile pyramid structure is set up
Tile map pyramid model is a kind of multiresolution hierarchical model, from the pyramidal bottom of tile to top layer, point Resolution is more and more lower, but the geographic range represented is constant, mainly includes two parts of tile and pyramid；Wherein：
1.1) tile
The tile be exactly according to a certain percentage size piece image is divided into some small square grids, each grid Lattice are exactly a tile；
Map tile after section is grating image, does not have location information, but related slices calculation has been used in section After method, the position specifically positioned can be calculated；
1.2) pyramid model
Tile is small square block, and length and width are equal, and length and width length is 256 pixels, 8~18 grades of totally 11 grades of engineer's scales；Ratio Ruler rank increases level-one；After having cut coding name is carried out according to the form of " l_r_c.png "；Wherein l is engineer's scale rank, R is line number, and c is row number, 0<=r<=2aa, 0<=c<=2aa；
Coding rule is cut according to above-mentioned picture, obtains the mutually exchange between picture ranks coordinate and actual geographic coordinate Calculate, these are the bases of map denotation；
Longitude and latitude is converted to line number and row number under some proportion ruler, and transfer process is：
N=Math.pow (2, lev)
X=((Longitude+180)/360) * n
SinLat=Math.Sin (Math.PI*Latitude/180)；
Y=0.5-Math.Log ((1+sinLat)/(1-sinLat))/(4*Math.PI)；
Wherein, Longitude is longitude, and Latitude is latitude value, and lev is level of zoom, and Math.pow is calculating 2 Lev power, Math.PI is pi, and Math.Sin is SIN function, and Math.Log is takes the logarithm；
Line number and row number under some proportion ruler are converted into latitude, and transfer process is：
N=Math.pow (2, lev)
Longitude=x/n*360-180；
Latitude=Math.Atan (Math.Sinh (Math.PI* (1-2*y/n)))；
Latitude=Latitude*180.0/Math.PI；
Wherein, Longitude is longitude, and Latitude is latitude value, and lev is level of zoom, and Math.pow is calculating 2 Lev power, Math.PI is pi, and Math.Atan is arctan function, and Math.Sinh is hyperbolic sine function；
2) DEM terrain datas are cut
Digital elevation model, vehicle economy M；It is a kind of entity that ground elevation is represented with one group of orderly array of values form Ground model, is a branch of digital terrain model, thus other various topographic index can derive from；It is generally believed that DTM It is various geomorphologic factors of the description including elevation, it is linear and non-thread including such as gradient, slope aspect, the change of slope factor Property combination spatial distribution, wherein DEM is the simple individual event digital land value model model of zeroth order, other as the gradient, slope aspect and the gradient become The landforms such as rate characteristic can derive from the basis of DEM；
Irregular triangle network, abbreviation TIN；The reconstruction of topographical surface is directly carried out using original sample point, by continuous phase The triangular facet composition being mutually coupled with, the shapes and sizes of triangular facet depend on density and the position of the observation station of irregular distribution；
The cutting of digital elevation data is that contour and elevational point are established irregular triangle network, then on the basis of irregular The small cubes that 256*256*256 pixels are cut into according to the layer of pyramid model, row, column establish dem data；Can be with different layers Secondary resolution ratio describes surface configuration；Can more accurately it be represented more with less room and time under a certain specified resolution Add the surface of complexity；Especially when landform includes big measure feature such as geosutures, construction line, these features can be preferably taken into account； The small cubes cut establish joint index according to level, row, column, block, and with the form storage of OBJ to database；
3) altitude data edge fit and inlay
3.1) DEM is edited
The DEM grid points that dem data editor refers to form interpolation are one by one into edlin；Identical terrain category DEM grid Point edge fit limit difference is twice of error in the terrain category grid points, and different terrain classification DEM edge fits limit difference is two kinds of landform classes The sum of other DEM grid points edge fit limit difference；
3.2) dem data edge fit
Choose the dem data that is generated of consecutive phantom, check the elevation of same level coordinate grid site, if occur elevation compared with Difference is more than the grid points of 2 times of DEM mean square error of height, then being considered as transfinites, and is regarded as rough error point, and re-establish body Model；It is right The dem data for rough error point occur carries out edge fit and repaiies edge fit again after survey；All single model DEM numbers are sequentially completed by the above process Between in succession；
3.3) dem data is inlayed and cut
If the edge fit of all single model dem datas is poor in fixed range all meets regulation requirement, numeral can be used The tile of orthography cutting inlays DEM；Height when inlaying to all same level coordinate net lattice points for participating in edge fit Journey takes its average value, as the height value of respective grid points, while forms the match accuracy report on each bar side；
After the completion of DEM inlays, carry out rectangle according to start-stop grid point coordinates as defined in relevant regulations or technical requirements and cut When, a row or multi-row DEM grid can be extended out according to specific requirement；
4) DEM quality controls
The quality control of DEM includes production process quality control and end result quality control two parts；
4.1) production process quality control
Whether the correctness and data that the key points of quality control inspection firsthand information in production process uses are reasonable；
Whether lake, reservoir, the layering in two-wire river are reasonable, and whether the elevation assignment in waterside line and coastline is rationally correct, DEM height values should be consistent in the range of static water body, flows the DEM height values in waters and gentle transition, relation should close from top to bottom Reason；
DEM elevations supposition area should be divided into by not reaching the region of regulation height accuracy requirement；
Grid positioned at white space should assign height value -9999, and the processing to clear area will be fully recorded in first number In；
4.2) end result quality control
The quality control of end result visually inspects whether contour has mutation feelings by the method for DEM interpolated contours Condition, or compared with topographic map, when landform shape, same place elevation differ greatly, then judge there are quality problems and repaiied Change；
Check the correctness of dem data start-stop point coordinates, check whether height value effective range is correct；
It should be checked after DEM splicings, determine whether overlapping and crack, whether splicing precision reaches requirement.

Claims (1)

1. one kind is based on detail three-dimensional map data processing method, it is characterised in that including 3 steps：Numeral is just penetrated Image cutting, the cutting of digital elevation data, altitude data is inlayed and edge fit；Wherein：
Digital orthoimage is cut：Digital orthoimage cutting is that satellite image picture is cut into 256* according to pyramid model The tile of 256 pixels, the tile cut carry out coding name according to the form of l_r_c.png, and l is the layer of pyramid model Number, r is the line number of pyramid model, and c is the row number of pyramid model；
Digital elevation data are cut：The cutting of digital elevation data is that contour and elevational point are established irregular triangle network, then The small cubes that 256*256*256 pixels are cut into according to pyramid model establish dem data；The small cubes cut according to Level, row, column, block establish index, and with the form storage of OBJ to database；
Altitude data edge fit and inlay：Based on the small cubes of digital elevation data cutting, the DEM that consecutive phantom is generated is chosen Data, check the elevation of same level coordinate grid site, if there are the poor grid points for being more than 2 times of DEM mean square error of height of elevation, Then being considered as transfinites, and is regarded as rough error point, and re-establish body Model；Survey is repaiied to there is the dem data of rough error point progress edge fit Again edge fit afterwards；It is sequentially completed by the above process between all single model dem datas in succession；If all lists in fixed range The edge fit of model dem data is poor all to meet regulation requirement, then the tile cut using digital orthoimage inlays DEM It is embedding；The elevation of all same level coordinate net lattice points when inlaying to participating in edge fit takes its average value, as respective grid points Height value;
This method uses following technological means：
1) tile pyramid structure is set up
Tile map pyramid model is a kind of multiresolution hierarchical model, from the pyramidal bottom of tile to top layer, resolution ratio It is more and more lower, but the geographic range represented is constant, includes two parts of tile and pyramid；Wherein：
1.1) tile
The tile be exactly according to a certain percentage size piece image is divided into some small square grids, each grid is just It is a tile；
Map tile after section is grating image, does not have a location information, but section used related slices algorithm it Afterwards, the position specifically positioned can be calculated；
1.2) pyramid model
Tile is small square block, and length and width are equal, and length and width length is 256 pixels, 8~18 grades of totally 11 grades of engineer's scales；Engineer's scale level Level-one is not increased；After having cut coding name is carried out according to the form of " l_r_c.png "；Wherein l is engineer's scale rank, and r is Line number, c are row number, 0<=r<=28, 0<=c<=28；
Coding rule is cut according to above-mentioned picture, obtains the mutual conversion between picture ranks coordinate and actual geographic coordinate, These are the bases of map denotation；
Longitude and latitude is converted to line number and row number under some proportion ruler, and transfer process is：
N=Math.pow (2, lev)
X=((Longitude+180)/360) * n
SinLat=Math.Sin (Math.PI*Latitude/180)；
Y=0.5-Math.Log ((1+sinLat)/(1-sinLat))/(4*Math.PI)；
Wherein, Longitude is longitude, and Latitude is latitude value, and lev is level of zoom, and Math.pow is to calculate 2 Lev power, Math.PI are pi, and Math.Sin is SIN function, and Math.Log is to take the logarithm；
Line number and row number under some proportion ruler are converted into latitude, and transfer process is：
N=Math.pow (2, lev)
Longitude=x/n*360-180；
Latitude=Math.Atan (Math.Sinh (Math.PI* (1-2*y/n)))；
Latitude=Latitude*180.0/Math.PI；
Wherein, Longitude is longitude, and Latitude is latitude value, and lev is level of zoom, and Math.pow is to calculate 2 Lev power, Math.PI are pi, and Math.Atan is arctan function, and Math.Sinh is hyperbolic sine function；
2) DEM terrain datas are cut
Digital elevation model, vehicle economy M；It is a kind of entity ground that ground elevation is represented with one group of orderly array of values form Model, is a branch of digital terrain model, and DTM is various geomorphologic factors of the description including elevation, the gradient, slope aspect, Including the change of slope factor linear processes combination spatial distribution, wherein DEM be the simple individual event of zeroth order digitally Looks model, other such as gradients, slope aspect and change of slope landforms characteristic derive from the basis of DEM；
Irregular triangle network, abbreviation TIN；The reconstruction of topographical surface is directly carried out using original sample point, is mutually interconnected by continuous The triangular facet composition connect, the shapes and sizes of triangular facet depend on density and the position of the observation station of irregular distribution；
The cutting of digital elevation data is that contour and elevational point are established irregular triangle network, then according to pyramid model The small cubes that layer, row, column cut into 256*256*256 pixels establish dem data；It can be described with the resolution ratio of different levels Surface configuration；More complicated surface can be more accurately represented under a specified resolution with less room and time；It is local When shape includes big measure feature geosutures, construction line, these features can be preferably taken into account；The small cubes cut are according to layer Level, row, column, block establish joint index, and with the form storage of OBJ to database；
3) altitude data edge fit and inlay
3.1) DEM is edited
The DEM grid points that dem data editor refers to form interpolation are one by one into edlin；Identical terrain category DEM grid points connect Margin difference is twice of error in the terrain category grid points, and different terrain classification DEM edge fits limit difference is two kinds of terrain category DEM The sum of grid points edge fit limit difference；
3.2) dem data edge fit
The dem data that consecutive phantom is generated is chosen, checks the elevation of same level coordinate grid site, if it is poor big elevation occur In the grid points of 2 times of DEM mean square error of height, then being considered as transfinites, and is regarded as rough error point, and re-establish body Model；To occurring The dem data of rough error point carries out edge fit and repaiies after survey edge fit again；Be sequentially completed by the above process all single model dem datas it Between in succession；
3.3) dem data is inlayed and cut
If the edge fit of all single model dem datas is poor in fixed range all meets regulation requirement, digital orthoimage is used The tile of cutting inlays DEM；The elevation of all same level coordinate net lattice points when inlaying to participating in edge fit takes it flat Average, as the height value of respective grid points, while forms the match accuracy report on each bar side；
After the completion of DEM inlays, according to as defined in relevant regulations or technical requirements start-stop grid point coordinates carry out rectangle cut when, root According to specific requirement to extend out a row or multi-row DEM grid；
4) DEM quality controls
The quality control of DEM includes production process quality control and end result quality control two parts；
4.1) production process quality control
Whether the correctness and data that the quality control checking firsthand information in production process uses are reasonable；
Whether lake, reservoir, the layering in two-wire river are reasonable, and whether the elevation assignment in waterside line and coastline is rationally correct, static Water body in the range of DEM height values should be consistent, flow waters DEM height values should gentle transition, relation be reasonable from top to bottom；
DEM elevations supposition area should be divided into by not reaching the region of regulation height accuracy requirement；
Grid positioned at white space should assign height value -9999, and the processing to clear area will be fully recorded in metadata；
4.2) end result quality control
Whether the quality control of end result has catastrophe by the methods of DEM interpolated contours, visual inspection contour, or Compared with topographic map, when landform shape, same place depth displacement different time, then judge there are quality problems and modify；
Check the correctness of dem data start-stop point coordinates, check whether height value effective range is correct；
It should be checked after DEM splicings, determine whether overlapping and crack, whether splicing precision reaches requirement.
CN201510979549.4A 2015-12-23 2015-12-23 One kind is based on detail three-dimensional map data processing method Active CN105427380B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510979549.4A CN105427380B (en) 2015-12-23 2015-12-23 One kind is based on detail three-dimensional map data processing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510979549.4A CN105427380B (en) 2015-12-23 2015-12-23 One kind is based on detail three-dimensional map data processing method

Publications (2)

Publication Number Publication Date
CN105427380A CN105427380A (en) 2016-03-23
CN105427380B true CN105427380B (en) 2018-05-08

Family

ID=55505559

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510979549.4A Active CN105427380B (en) 2015-12-23 2015-12-23 One kind is based on detail three-dimensional map data processing method

Country Status (1)

Country Link
CN (1) CN105427380B (en)

Families Citing this family (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105790771A (en) * 2016-03-29 2016-07-20 黄河勘测规划设计有限公司 3D topographic data compression method of superhigh compression ratio
CN105956124B (en) * 2016-05-06 2019-07-23 中国神华能源股份有限公司 A kind of map methods of exhibiting and device
CN106204448B (en) * 2016-07-04 2019-09-13 北京空间飞行器总体设计部 The stepless zooming display methods and display system of the big figure of spacecraft stratification power supply and distribution
CN107705241B (en) * 2016-08-08 2021-08-03 国网新疆电力公司 Sand table construction method based on tile terrain modeling and projection correction
CN106651927B (en) * 2016-12-30 2020-03-24 北京航天泰坦科技股份有限公司 Homonymous point gross error eliminating method for matching embedded edges of orthophoto images
CN106991143B (en) * 2017-03-22 2019-07-19 苏州中科图新网络科技有限公司 Multi-layer image file, generation method and device, read method and device
CN108733660A (en) * 2017-04-13 2018-11-02 北大方正集团有限公司 Geographic information processing method, apparatus, server and system
CN107220499B (en) * 2017-05-26 2021-04-13 粉蓝医疗科技（杭州）有限公司 Slice file storage method and device
CN107393008B (en) * 2017-05-31 2020-09-15 陈树铭 System and method for modeling irregular triangular net pyramid of global geographic frame
CN109427063A (en) * 2017-08-30 2019-03-05 广州极飞科技有限公司 The storage method and device of image data
CN107992588B (en) * 2017-12-11 2020-08-25 国网技术学院 Terrain display system based on elevation tile data
CN109933559A (en) * 2017-12-15 2019-06-25 国网上海市电力公司 A kind of data displaying optimization method based on power business
CN108052642A (en) * 2017-12-22 2018-05-18 重庆邮电大学 Electronic Chart Display method based on tile technology
CN108108575B (en) * 2018-01-17 2021-06-01 河南省水利勘测设计研究有限公司 High-precision engineering three-dimensional terrain manufacturing method based on BIM (building information modeling) cross-platform technology
CN108053475A (en) * 2018-02-09 2018-05-18 城市生活（北京）资讯有限公司 A kind of three-dimensional terrain display method and device
CN108830929A (en) * 2018-05-21 2018-11-16 东南大学 Multi-resolution Terrain pyramid model generation method and system based on database
CN109118563B (en) * 2018-07-13 2023-05-02 深圳供电局有限公司 Method for extracting digital orthographic image from LOD paging surface texture model
CN109636718B (en) * 2018-10-31 2021-04-20 百度在线网络技术（北京）有限公司 Method, device and equipment for detecting point cloud splicing quality and storage medium
CN109657024A (en) * 2018-12-12 2019-04-19 武汉中科天宇科技有限公司 A kind of storage of three-dimensional map data, inquiry and display systems and corresponding method
CN109753551B (en) * 2019-01-28 2021-03-16 中国科学院遥感与数字地球研究所 Data organization and storage method and device for dynamic geographic phenomena
CN109920042A (en) * 2019-02-19 2019-06-21 网易（杭州）网络有限公司 Virtual landforms generation method, device and readable storage medium storing program for executing
CN110134752B (en) * 2019-05-16 2021-01-15 洛阳众智软件科技股份有限公司 Three-dimensional large-scene modeling data processing method and device
CN110310367A (en) * 2019-07-03 2019-10-08 武汉大势智慧科技有限公司 Based on large scene outdoor scene three-dimensional multi-angle 2.5D image lightweight browsing method
CN110347769B (en) * 2019-07-09 2021-03-30 东方网力科技股份有限公司 Processing method, device, equipment and storage medium for multi-level map tiles
CN110544306B (en) * 2019-08-29 2023-06-27 中国南方电网有限责任公司 Three-dimensional terrain data acquisition method, three-dimensional terrain data acquisition device, computer equipment and storage medium
CN111060075B (en) * 2019-12-10 2021-01-12 中国人民解放军军事科学院国防科技创新研究院 Local area terrain ortho-image rapid construction method and system based on unmanned aerial vehicle
CN113066178B (en) * 2020-01-02 2023-08-15 沈阳美行科技股份有限公司 Map data processing method, device, equipment and storage medium
CN112084283B (en) * 2020-09-11 2022-03-22 广州南方智能技术有限公司 Terrain image slicing method based on pyramid structure and hierarchical detail model
CN112530009A (en) * 2020-11-06 2021-03-19 四川见山科技有限责任公司 Three-dimensional topographic map drawing method and system
CN112529361B (en) * 2020-11-13 2024-03-22 许昌华杰公路勘察设计有限责任公司 Highway investigation route selection method based on smart phone and digital topography

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101713649A (en) * 2009-11-05 2010-05-26 中国测绘科学研究院 Disturbing gravity-based quasi-geoid land-sea seamless splicing method
CN101887595A (en) * 2009-05-14 2010-11-17 武汉如临其境科技创意有限公司 Three-dimensional digital earth-space data organizing and rendering method based on quad-tree index
CN102117498A (en) * 2009-12-30 2011-07-06 新奥特（北京）视频技术有限公司 Method and device for processing data block joint in three-dimensional modeling
CN102117500A (en) * 2009-12-30 2011-07-06 新奥特（北京）视频技术有限公司 Three-dimensional modeling method and system
CN102117492A (en) * 2009-12-30 2011-07-06 新奥特（北京）视频技术有限公司 Method and device for modeling partitioned data
CN102117493A (en) * 2009-12-30 2011-07-06 新奥特（北京）视频技术有限公司 Method and device for updating data in three-dimensional modeling
CN102117495A (en) * 2009-12-30 2011-07-06 新奥特（北京）视频技术有限公司 New data partitioning method and device in three-dimensional modeling
CN102750294A (en) * 2011-08-10 2012-10-24 新奥特（北京）视频技术有限公司 Method, device and system for generating topography image data
CN103425801A (en) * 2013-09-04 2013-12-04 中测新图(北京)遥感技术有限责任公司 Digital-earth-oriented data fusion method and device
CN103823981A (en) * 2014-02-28 2014-05-28 武汉大学 DEM (Digital Elevation Model)-assisted satellite image block adjustment method
CN105659820B (en) * 2011-09-21 2014-06-18 北京航天飞行控制中心 A kind of based on planetary vehicle vision and the multi-modal localization method of satellite remote sensing images landing point

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101887595A (en) * 2009-05-14 2010-11-17 武汉如临其境科技创意有限公司 Three-dimensional digital earth-space data organizing and rendering method based on quad-tree index
CN101713649A (en) * 2009-11-05 2010-05-26 中国测绘科学研究院 Disturbing gravity-based quasi-geoid land-sea seamless splicing method
CN102117498A (en) * 2009-12-30 2011-07-06 新奥特（北京）视频技术有限公司 Method and device for processing data block joint in three-dimensional modeling
CN102117500A (en) * 2009-12-30 2011-07-06 新奥特（北京）视频技术有限公司 Three-dimensional modeling method and system
CN102117492A (en) * 2009-12-30 2011-07-06 新奥特（北京）视频技术有限公司 Method and device for modeling partitioned data
CN102117493A (en) * 2009-12-30 2011-07-06 新奥特（北京）视频技术有限公司 Method and device for updating data in three-dimensional modeling
CN102117495A (en) * 2009-12-30 2011-07-06 新奥特（北京）视频技术有限公司 New data partitioning method and device in three-dimensional modeling
CN102750294A (en) * 2011-08-10 2012-10-24 新奥特（北京）视频技术有限公司 Method, device and system for generating topography image data
CN105659820B (en) * 2011-09-21 2014-06-18 北京航天飞行控制中心 A kind of based on planetary vehicle vision and the multi-modal localization method of satellite remote sensing images landing point
CN103425801A (en) * 2013-09-04 2013-12-04 中测新图(北京)遥感技术有限责任公司 Digital-earth-oriented data fusion method and device
CN103823981A (en) * 2014-02-28 2014-05-28 武汉大学 DEM (Digital Elevation Model)-assisted satellite image block adjustment method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《TIN与等高线的生成方法》;王野等;《科学发展与社会责任（A卷）——第五届沈阳科学学术年会文集》;20081013;268-273 *
《栅格型空间数据服务体系结构与算法研究》;邓雪清;《中国优秀博硕士学位论文全文数据库 (博士) 基础科学辑》;20031215;第2003年卷(第04期);A008-3 *
《电力信息***中三维GIS关键技术的应用研究》;魏祖宽等;《计算机与现代化》;20100515;第2010年卷(第5期);83-88 *

Also Published As

Publication number Publication date
CN105427380A (en) 2016-03-23

Similar Documents

Publication Publication Date Title
CN105427380B (en) One kind is based on detail three-dimensional map data processing method
Danielson et al. Topobathymetric elevation model development using a new methodology: Coastal national elevation database
Kumler An intensive comparison of triangulated irregular networks (TINs) and digital elevation models (DEMs)
Ahmad et al. Interactive Geoinformation Three-Dimensional Model of a Landscape Park Using Geoinformatics Tools
Van Westen et al. Multi-hazard risk assessment: Distance education course-Risk City Exercise book 2011
CN112184908A (en) 3D Tiles format model bounding box data generation method for realizing oblique photography model based on Cesum
CN108389255A (en) Landform geometric parameter extracting method based on layering elevation cloud atlas
CN104933899A (en) Ship seaworthiness early-warning system and method
CN103065361A (en) Three-dimensional (3d) island sandbox achieving method
CN105354882A (en) Method for constructing big data architecture based three-dimensional panoramic display platform for large-spatial-range electricity transmission
CN109063193A (en) A kind of thematic maps recommends the method and device of display
CN115687675A (en) Airborne vector map data processing method
KR101214996B1 (en) Method for serving and constructing a sea navigation information
Stupariu et al. Geometric approaches to computing 3D-landscape metrics
JP2009251250A (en) Numerical map data processing method, numerical map data processing program, and numerical map data processing apparatus
Pacina et al. Pre‐dam valley reconstruction based on archival spatial data sources: Methods, accuracy, and 3D printing possibilities
Hall et al. Assessment of natural stream sites for hydroelectric dams in the pacific northwest region
Beattie 3D visualization models as a tool for reconstructing the historical landscape of the Ballona Creek watershed
Jansson et al. Uncertainty in area determination
Yang et al. Predicting flood inundation and risk using geographic information system and hydrodynamic model
Parker et al. Blending bathymetry with topography: The Tampa Bay demonstration project
Devátý AUTOMATED ASSESSMENT OF PROFILES FOR BUILDING SMALL WATER RESERVOIRS
KR20230170430A (en) System, method and computer program for making three-dimensional geospatial model of river
Ismail et al. A REVIEW ON RUNOFF ESTIMATION OF VRISHABHAVATHI WATERSHED USING GIS
Zheng Continental-Scale High-Resolution River Geometry and Real-Time Inundation Mapping

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant