[go: up one dir, main page]

CN101887595B - Three-dimensional digital earth-space data organizing and rendering method based on quad-tree index - Google Patents

Three-dimensional digital earth-space data organizing and rendering method based on quad-tree index Download PDF

Info

Publication number
CN101887595B
CN101887595B CN200910062052.0A CN200910062052A CN101887595B CN 101887595 B CN101887595 B CN 101887595B CN 200910062052 A CN200910062052 A CN 200910062052A CN 101887595 B CN101887595 B CN 101887595B
Authority
CN
China
Prior art keywords
data
dimensional
map
quad
digital earth
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.)
Expired - Fee Related
Application number
CN200910062052.0A
Other languages
Chinese (zh)
Other versions
CN101887595A (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.)
WUHAN IMMERSION ENVIRONMENT CO Ltd
Original Assignee
WUHAN IMMERSION ENVIRONMENT 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 WUHAN IMMERSION ENVIRONMENT CO Ltd filed Critical WUHAN IMMERSION ENVIRONMENT CO Ltd
Priority to CN200910062052.0A priority Critical patent/CN101887595B/en
Publication of CN101887595A publication Critical patent/CN101887595A/en
Application granted granted Critical
Publication of CN101887595B publication Critical patent/CN101887595B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

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

Abstract

The invention provides a three-dimensional digital earth-space data organizing and rendering method based on a quad-tree index, which belongs to the technical fields of cartography, geography information systems and virtual reality. The method comprises the following steps of: unifying common spatial data formats under multi-scale and multi-projection conversion, multispectral, multi-temporal and high-resolution remote sensing satellite images and aerial images as well as digital thematic maps with different scales into same coordinate system, carrying out operation on the attribute of each element and the parameter regulation of quad-tree tiles, outputting in the form of the quad-tree tiles, carrying out quad-tree cutting on three-dimensional landscape map data, leading spatial data into a relevant database, and carrying out unified management. By using the method, common vector data, raster data, altitude data and three-dimensional map data are organically integrated and issued into a three-dimensional digital earth prototype, thereby shortening the time of data preprocessing, improving the execution efficiency and providing a new integration method for three-dimensional digital earth fundamental geographic data dissemination.

Description

Three-dimensional digital earth space data organization rendering method based on quadtree index
Technical Field
The invention belongs to the technical field of geographic information systems, cartography and virtual reality, and particularly relates to a three-dimensional digital geospatial data organization rendering method based on a quadtree index.
Background
The digital earth is based on computer technology, multimedia technology and large-scale storage technology, and carries out multi-resolution, multi-scale, multi-space-time and multi-variety three-dimensional description on the earth by using massive earth information by taking a broadband network as a link.
The digital earth is mainly composed of spatial data, text data, an operation platform and an application model. The space data mainly comprises global space data with large, medium and small scales, various multispectral, multi-temporal and high-resolution remote sensing satellite images, aerial images and various digital thematic maps with different scales of the earth. A basic idea of the digital earth is to cover a high-resolution remote sensing image on the whole globe, and establish an earth landscape model by establishing a digital elevation model and vector data covering the whole globe to form a virtual earth.
The spatial data is mainly realized by combining various different types of spatial data. The surface area of the earth is close to (1.5 multiplied by 10)12) The data volume of the three-dimensional digital earth is extremely large and the data type is complex due to square kilometers, and the data volume increases in a geometric progression along with the improvement of resolution and often reaches hundreds of G or even tens of TB bytes. In the current network environment or on a PC, the huge data is difficult to meet the requirements of real-time transmission and display, so that the transmission of massive geographic data and the real-time drawing of three-dimensional scenes become important contents in digital earth research. The updating of computer storage equipment and the reduction of price provide support for large-capacity data storage, and for a network-based system, the method reasonably organizes spatial data, improves the access speed of mass data, reduces unnecessary data transmission on the network as far as possible, and has important significance for improving the display rate of global three-dimensional scenes.
The global terrain data volume is very large, for example, 1: 25 ten thousand terrain data in the world has a total data volume of more than 20G, and the data volume is conceivable together with images. In order to reduce the consumption of computer memory, accelerate computer processing time, save network transmission time and access multi-resolution ultra-large geographic data at the fastest speed, a reasonable database storage structure must be designed for flexibly managing data from coarse to fine, and such an adaptive database structure must be implemented by using spatial indexes. Data is typically organized hierarchically, with the data at each level being divided into small chunks, and this method of organizing data is commonly referred to as a quadtree structure or a pyramid data structure. FIG. 1 illustrates the structural hierarchy of a quadtree.
Reference to the literature
[1] The basic institution of national spatial information deals with the digital earth [ M ]. beijing: qinghua university Press, 1999.
[2] Goji keyya, duchesnea, li qingquan, etc. contemporary geographic information technology [ M ]. beijing: scientific press, 2004: 87 to 88.
[3] The fourth generation of GIS software research, china graphics press, 2001b, 6A (9): 817 to 823.
[4] Songkefu, bells, research and development of modular geographic information systems, graphic presentations of images, 1998, 4: 314.31.
[5] the method comprises the steps of child dawn thoroughfare, cardia admission, Zhang Yong, construction and rapid display of a global multi-resolution data model, surveying and mapping science 2006.1.
[6] White-built military, Zhao national wins, Chen military based on the global terrain visualization of ellipsoid triangular network, Wuhan university newspaper, 2005.05.
[7] Li de ren, information highway, space data infrastructure and digital earth, proceedings of surveying and mapping, 1999.
[8] Li de ren, development of remote sensing and development of century, university of wuhan, 2003.02.
[9] Wu Yanlan, DEM and several applications, Wuhan university of surveying and mapping science and technology academic papers, 1998.
[10] Penhui, establishment of a digital elevation model based on MicroStation and application research thereof, academic paper of Wuhan university of surveying and mapping technology, 1998.
[11] Xue Yong, Wangjiaqin, Guo Huadong digital Earth grid computing was young [ J ]. remote sensing academic, 2004.8(I).
[12] Yangyong, Guodazhi, digital map calculation model [ J ] based on geographic coordinates, university of southwest traffic, 6 months 2005, Vol 40, No3.
Disclosure of Invention
The invention aims to provide a three-dimensional digital geospatial data organization rendering method based on a quadtree index. The method can unify common spatial data formats under global multi-scale and multi-projection changes, various multispectral, multi-temporal and high-resolution remote sensing satellite images, aerial images and various digital thematic maps of different scales of the earth into the same coordinate system, can perform operations such as attribute adjustment and quadtree tile parameter adjustment on each element, finally outputs visual data in a quadtree pyramid tile form, performs quadtree cutting on three-dimensional landscape map data, and imports geometric data and attribute data of the spatial data into a relational database for unified management. The invention provides a set of detection framework for the final data, and ensures that the data can be stably and correctly provided for the digital earth platform.
The technical scheme of the invention is a three-dimensional digital earth space data organization rendering method based on the quadtree index, which comprises the following design ideas (or the following basic steps):
1) hierarchical blocking mechanism using digital earth data
2) Generating and rendering tiles
3) Fusing multi-source data
4) Detection tool adopting three-dimensional digital earth prototype data
The invention has the beneficial effects that: by using the method, common vector data, raster data, elevation data, three-dimensional map data and the like can be organically fused and issued to the three-dimensional digital earth prototype in a split storage mode of a quadtree pyramid, so that the time for preprocessing the data is obviously shortened, the execution efficiency is improved, and a new integration method is provided for issuing the three-dimensional digital earth basic geographic data.
Drawings
FIG. 1 structural hierarchy of quadtrees
FIG. 2 Cartesian coordinate System
Width and height of the tiles of FIG. 3
Graph 436 degree division
FIG. 5 is a schematic diagram of different levels of division
FIG. 6 level quadtree tile chunking representation
FIG. 7 segmentation rules for layer n quad tree to layer n +1 quad tree
FIG. 8 Tile Pre-load and transition display processing
FIG. 9 map data Source
FIG. 10 three-dimensional landscape map rendering process in 3Ds Max
FIG. 11 single map cutting process
FIG. 12 multiple-map-related cutting process for super large map
FIG. 13 is a hierarchical relationship between linear vector layers and DEM and image
FIG. 14 supports geodetic to Lambert projection coordinates
FIG. 15 supports Jpeg2000 raster data
FIG. 16 supports DEM data
FIG. 17 simple element data model for OpenGIS
FIG. 18 Tile rendering method tool prototype screenshot
FIG. 19 Tile rendering method tool prototypes export data to three-dimensional digital Earth prototypes
FIG. 20 Multi-layer grid fusion with vector data
Detailed Description
The design idea of the invention is described below with reference to the drawings of the specification.
1 digital earth data layered block division mechanism
1.1 surface Tile partitioning mechanism
Cartesian coordinates are taken with the origin (X-0, Y-0.) at the lower left of the projection coordinates, i.e., south poles (-90, -180) (in longitude and latitude), see fig. 2 and 3.
As shown in FIG. 3, the Size of each Tile width and height is determined using "Level Zero Tile Size" (all tiles are squares. the standard Level Zero Tile Size has not been established, but it must be satisfied that it is divisible by 180.
The first layer is divided by 36 degrees, see fig. 4, with 360 degrees longitude and 180 degrees latitude of the earth as the criteria.
The total number of layers is (360/36) × (180/36) ═ 50 layers, the second layer is 18 degrees, the third layer is 9 degrees, and the third layer is similarly arranged, as shown in fig. 5 and 6.
1.2 Tile positioning mechanism
Defining:
A. the map level is from 1, and the number of rows and columns is from 0;
B. extreme values of the map range are expressed by XMin, XMax, YMin and YMax;
C. x Δ represents an abscissa difference, and y Δ represents an ordinate difference;
D. naming of the quad-tree map tiles by adopting a mode of containing map levels and row and column numbers, such as 'zoomx _ row _ column', wherein x represents the map level, and row and column analysis represents the row number and the column number of the current tile;
1.2.1 map quantity calculation
If the first level map is m rows and n columns, the number of rows and columns of the map at the second level is as follows:
line number: mx 2level-1,level≥1
The number of columns: n x 2level-1,level≥1
1.2.2 conversion of geographic coordinates to Tile position
1.2.2.1 calculating geographic range from Tile location
Inputting: rank pair (i, j)
And (3) outputting: extent Range (currentXMin, currentXMax, currentYMin, currentYMax)
The coordinate range of the ith row and j column position slice at the level of the level is calculated by the following formula:
<math> <mrow> <mi>currentXMin</mi> <mo>=</mo> <mi>XMin</mi> <mo>+</mo> <mi>j</mi> <mo>&times;</mo> <mfrac> <mi>&Delta;x</mi> <mrow> <mi>m</mi> <mo>&times;</mo> <msup> <mn>2</mn> <mrow> <mi>level</mi> <mo>-</mo> <mn>1</mn> </mrow> </msup> </mrow> </mfrac> </mrow> </math>
<math> <mrow> <mi>currentXMax</mi> <mo>=</mo> <mi>XMin</mi> <mo>+</mo> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&times;</mo> <mfrac> <mi>&Delta;x</mi> <mrow> <mi>m</mi> <mo>&times;</mo> <msup> <mn>2</mn> <mrow> <mi>level</mi> <mo>-</mo> <mn>1</mn> </mrow> </msup> </mrow> </mfrac> </mrow> </math>
<math> <mrow> <mi>currentYMin</mi> <mo>=</mo> <mi>YMax</mi> <mo>-</mo> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&times;</mo> <mfrac> <mi>&Delta;y</mi> <mrow> <mi>n</mi> <mo>&times;</mo> <msup> <mn>2</mn> <mrow> <mi>level</mi> <mo>-</mo> <mn>1</mn> </mrow> </msup> </mrow> </mfrac> </mrow> </math>
<math> <mrow> <mi>currentYMax</mi> <mo>=</mo> <mi>YMax</mi> <mo>-</mo> <mi>i</mi> <mo>&times;</mo> <mfrac> <mi>&Delta;y</mi> <mrow> <mi>n</mi> <mo>&times;</mo> <msup> <mn>2</mn> <mrow> <mi>level</mi> <mo>-</mo> <mn>1</mn> </mrow> </msup> </mrow> </mfrac> </mrow> </math>
1.2.2.2 calculating Tile position from coordinates
Inputting: (x, y) coordinates
And (3) outputting: i, j (number of rows and columns of picture)
Wherein,
<math> <mrow> <mi>j</mi> <mo>=</mo> <msqrt> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <msub> <mi>X</mi> <mi>Min</mi> </msub> <mo>)</mo> </mrow> <mo>/</mo> <mi>&Delta;x</mi> <mo>&times;</mo> <mi>m</mi> <mo>&times;</mo> <msup> <mn>2</mn> <mrow> <mi>level</mi> <mo>-</mo> <mn>1</mn> </mrow> </msup> </msqrt> </mrow> </math>
<math> <mrow> <mi>i</mi> <mo>=</mo> <msqrt> <mrow> <mo>(</mo> <msub> <mi>Y</mi> <mi>Max</mi> </msub> <mo>-</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>/</mo> <mi>&Delta;y</mi> <mo>&times;</mo> <mi>n</mi> <mo>&times;</mo> <msup> <mn>2</mn> <mrow> <mi>level</mi> <mo>-</mo> <mn>1</mn> </mrow> </msup> </msqrt> </mrow> </math>
with the data structure of the quadtree, referring to fig. 4, Level 0 divides the earth into 50 regions at 36 degrees, each block at the 0 th layer is divided into four subtrees, and the Level 1 layer is calculated as shown in fig. 7.
Each block has four subtrees, including its own information, and the data is defined as:
double West; v/four directional planar Cartesian coordinate values
double East;
double North;
double South;
Angle CenterLatitude; // tile center latitude and longitude
Angle CenterLongitude;
double LatitudeSpan;
double LongitudeSpan;
int Level; // Level
int Row; // row, column
int Col; // subtree
QuadTile northWestChild;
QuadTile southWestChild;
QuadTile northEastChild;
QuadTile southEastChild;
Each block corresponds to a 512 x 512 texture of the current Level, each edge of each block is divided into 40 equal parts, 41 x 41 or 43 x 43 vertexes are generated, and then the real-time drawing is performed by combining the height information.
Due to the layered and partitioned structure, when terrain data is loaded and rendered in real time, texture pictures and elevation information of a required layer and an area only need to be downloaded from a server, and some pictures can be preprocessed to form direct pictures of three-dimensional fruits without elevation information.
Transition treatment:
the Level (n +1) Level does not download the unprocessed pictures, and still displays the effect of the Level n Level, as shown in fig. 8.
2 generating tiles
2.1 three-dimensional Tile data
The map pre-generation is to cut a map in a specified range into square pictures in a plurality of rows and columns according to a specified size and a specified format (such as TIFF, JPEG, PNG, etc.), a map slice obtained by cutting the map is called a Tile (Tile), and the Tile acquisition flow is shown in fig. 9. Rendering the tiles to obtain only raster graphs corresponding to two-dimensional vector data, image data or digital elevation data, performing three-dimensional modeling on the cut graphs of the three-dimensional landscape map by using tools such as AutoCAD (auto computer aided design), 3DMax and the like, adjusting the cut graphs by using a VRay plug-in, and rendering the whole three-dimensional map.
The three-dimension of the three-dimensional electronic map is 2.5-dimension with rich observation surface, and is a static rendering picture with fixed view angle generated by a three-dimensional model. Before manufacturing, the angle of a camera is determined according to the orientation characteristics of a regional house, the camera and the lighting effect are set in three-dimensional model manufacturing software, a base map is manufactured, and a standard scene is formed and is used by a manufacturer in a unified mode. Objects that are occluded, invisible or in secondary positions may be omitted or simplified. And point, line and surface city entity information layers in the three-dimensional scene are manufactured in a layering mode during model manufacturing, and the three layers of models are overlapped to form the city three-dimensional landscape. Fig. 10 shows an implementation process of rendering a three-dimensional model into a landscape picture in 3d scax.
2.2 Tile cutting
The Net class library has a limit to the image volume and limits the length and width of a loaded single map to be within 15000 × 15000, so that the map cutting is divided into two types of single map cutting and multiple related cutting of an oversized map:
1. single map cutting process:
step 1: correcting the original value of the cutting range to make the length and height values of the original value be multiples of the tile;
step 2: the fixed range is cut from the upper left corner (northwest), a quad tree data structure is adopted among subsequent maps of different levels, a Tile (Tile) on the level two is split into 4 tiles from the level +1, the structure is beneficial to cutting and displaying, but the obtained map has no fixed scale, and the scale changes along with the vertical coordinates of the geography, so that the geography calculation is directly calculated according to the geography coordinates instead of the scale. The map naming rule is that the Nth level map picture is named as zoomN.
The tile naming rule of the x-th level chunk row chunk column is zomx _ chunk.
Wherein the variables define:
x: a current zoom level;
chunkrow: the number of rows of the current tile in the map;
chunkcol: the number of columns of the current tile in the map;
2. a large map multiple correlation cutting process:
step 1: the method comprises the steps of dividing a super large map into a plurality of mutually related small maps, wherein the naming rule of the small maps is uperchunk _ row _ col \ zoomn.
Step 2: cutting the small map by adopting a single map cutting process
The naming rule of the tile after cutting is hyperchronk _ row _ col | zoomx _ chunkrow _ chunkcol \ x _ chunkrow _ chunkcol.
Wherein the variables define:
row: the current small map occupies the row position in all the cut small maps;
and (5) Col: the current small map occupies the column position in all the cut small maps;
x: a current zoom level;
chunkrow: the number of rows of the current tile in the map;
chunkcol: the number of columns of the current tile in the map;
single map cutting process/multiple associated cutting processes for super large maps, as shown in fig. 11 and 12.
3 Multi-source data fusion
3.1 fusing data layers
With the difference of the distance of the viewpoints, the definition of the three-dimensional terrain changes in real time, and at the moment, the system automatically calls DEMs and image data with different resolutions to construct a ground model. When the terrain transits from one detail level to another detail level, the system calls a plurality of terrain data displays at the same time, and the segmentation of the terrain patches changes. If the ground features are not fused with the terrain, the space vector entity can be caused to have the phenomena of 'sinking' and 'lifting'. In order to increase the display speed and the graphic quality, the data of the vector layer is hierarchically managed as the DEM and the image, as shown in fig. 13.
3.2 support projective transformation
By reading the. prj file corresponding to the.shp file and parsing in the WKT format standard, fig. 14.
The parameters include:
the coordinate system is described in terms of OpenGIS WKT (Well known Text). It includes the following information:
A. the name of a general coordinate system.
B. The name of a geographic coordinate system.
C. The earth frame of reference.
D. Reference ellipsoid, ellipsoid semi-major axis and. The derivative of the ellipsoidal ellipticity (i.e., a/(a-b)).
E. A 0 degree meridian name car and its degree of deviation from this initial meridian.
F. One projection type.
G. A list of projection parameters.
H. A unit name and a conversion factor to meters or radians.
I. The name and ordering of the axes.
3.3 multiple vector and raster data type support
The OpenGIS standard is used as a uniform interface, and different data can be supported by the same platform, such as
■ vector data: shape data, mapinfo data, gml data, geojson data
■ raster data: geotiff data, ecw data, jp2 data (see FIG. 15)
■ digital elevation data: DEM data (see FIG. 16)
■ GPS data: GPX format
The GDAL library is used for realizing the support of vector data such as Shape data, mapinfo data, gml data, geojson data and the like.
The vector data is based on a rectangular coordinate system, and a data model or a data structure is obtained by describing geographic elements by points, lines and polygons. Each geographic element has a series of sequential coordinate (x, y) descriptions, which are combined with attributes. We have designed a Geometry model that strictly follows the simple element data specification of the OGC OpenGIS, where basic Geometry in vector graphics is defined, and some simple spatial correlation analysis functions. The most important geometries are points, lines, facets, and thus derived multi-points, multi-lines, multi-facets, curves, lines, loops, polygons, etc., see fig. 17.
The relationship of grid position (pixel coordinates or line coordinates) and geographic reference coordinates is described in two ways. The first, also the most commonly applied one, is the affine transformation.
The affine transform consists of six coefficients, which map the pixel/line coordinates to the geographic reference space with the following relation:
Xgeo=GT(0)+Xpixel*GT(1)+Yline*GT(2)
Ygeo=GT(3)+Xpixel*GT(4)+Yline*GT(5)
for the downward-ordered images, GT (2) and GT (4) are both zero, GT (1) is pixel wide, GT (5) is pixel high, and GT (0), GT (3)) represents the upper left coordinates of the upper left pixels of the grid.
Note that: the upper pixel and line coordinates are from the upper left coordinate (0.0 ) of the upper left pixel to the lower right coordinate (pixel wide, pixel high) of the lower right pixel. The center position of the top left pixel/line is then (0.5 ).
Defining: wave band
The raster data for one band corresponds to the GDALRasterBand class in GDAL. It describes the band/channel/layer of a single band. It does not describe the entire image at once. For example, a 24-bit RGB image is generally described in a dataset as three bands, corresponding to the red/green/blue colors.
Each band is specifically interpreted as:
GCI _ defined: default, unknown information.
GCI _ gray index: is an independent gray-scale image.
GCI_PaletteIndex:this raster acts as an index into a color table。
GCI _ RedBand: the red band of an RGB or RGBA image.
GCI _ GreenBand: the green band of an RGB or RGBA image.
GCI _ BlueBand: the blue band of an RGB or RGBA image.
GCI _ AlphaBand: the alpha channel of the RGBA image.
GCI _ hue band: hue of HLS image.
GCI _ SaturationBand: saturation of HLS images.
GCI _ LightnessBand: light intensity of HLS image.
GCI _ CyanBand: CMY or CMYK images.
GCI _ magenta band: magenta bands for CMY or CMYK images.
GCI _ YellowBand: yellow bands for CMY or CMYK images.
GCI _ BlackBand: black band of CMYK images.
Support data attribute tuning: supporting data attribute tuning
3.4 quadtree tile parameter adjustment
Adjusting the tile according to the quadtree parameters, then exporting the tile, and exporting a parameter configuration file, wherein the format is as follows:
<MapParam>
<tilesize>128</tilesize>
<ColumnNum>4</ColumnNum>
<RowNum>2</RowNum>
<centerX>-243383.913445131</centerX>
<centerY>4376213.31793891</centerY>
<ZoomParam>2</ZoomParam>
<MaxLevel>4</MaxLevel>
<Xmin>-4663419.25215648</Xmin>
<Xmax>4176651.42526622</Xmax>
<Ymin>2166195.64858323</Ymin>
<Ymax>6586230.98729458</Ymax>
</MapParam>
3.5 building of spatial database
A. Database connection
Comprises the following steps: checking in a Sql mode; and ② Windows identity authentication.
B. A new database is built, including database files (. MDF) and log files (. LDF).
C. Spatial database deployment: fields "blue. general _ COLUMNS" and "blue. spatial _ REF _ SYS" are generated in the specified database and used for storing the geometry layer and the spatial WKT information, respectively.
D. Vector data import
And importing the shape data into a database. The shp file generates fixed fields oid, the _ geom, zID, the _ geom _ Envelope _ MinX, the _ geom _ Envelope _ MinY, the _ geom _ Envelope _ MaxX, and the _ geom _ Envelope _ MaxY, which respectively store a stream number, geometric feature data, a map layer number, a bounding box abscissa minimum value, a bounding box ordinate minimum value, a bounding box abscissa maximum value, and a bounding box ordinate maximum value. From the dbf file, attribute information fields are generated, such as: "city" and "population".
3.6 geographic data bulk search
And searching the corresponding longitude and latitude by inputting the geographic position based on the longitude and latitude search of the Geocoding. Transmission format is "http://ditu.google.com/maps/geooutput=xml&qAnd analyzing the returned xml format file according to the request of the geographic position, wherein point. By using datasets and data adapters, a batch of data search is conducted based on sqlserver 2005.
4 the map is edited by the method and the data is exported to the digital earth for demonstration, see fig. 18, fig. 19 and fig. 20.

Claims (2)

1. A three-dimensional digital earth space data organization rendering method based on a quadtree index comprises the following steps:
(1) performing mutual operation on geographic coordinates and tiles by adopting a digital earth data hierarchical partitioning and tile positioning mechanism based on a quadtree pyramid, and calculating a geographic range according to tile positions;
(2) rendering a three-dimensional map based on a quadtree pyramid, layering and superposing three-dimensional models to form urban three-dimensional landscape, and cutting the rendered three-dimensional map into a single map and a plurality of map tiles of an oversized map;
(3) establishing a two-dimensional Cartesian coordinate system to support multi-source data fusion, abstracting vector data into multi-point, multi-line, multi-surface, curve, straight line, circular line and polygon for storage fusion based on a Geometry model following simple element data specification of an OGC OpenGIS, performing affine transformation on raster data according to the coordinate system, extracting data fusion of red, green and blue three wave bands, reading a universal wave band of elevation data, converting various projection transformations into a unified coordinate system of geodetic coordinates, and performing fusion processing on the multi-source data.
2. The method of claim 1, wherein in step (1), the earth is divided into 50 regions at level 0 by 36 degrees, each region corresponds to a 512 × 512 texture of the current Leve1, each edge of each block is divided into 40 equal parts, and 41 × 41 vertices are generated;
in the step (2), the city entity information layers of points, lines and surfaces in the three-dimensional scene are overlapped to form a city three-dimensional landscape during model making, and the rendered three-dimensional map is cut to ensure that the length range of the loaded single map is within 15000 x 15000;
in the step (3), a two-dimensional Cartesian coordinate system is established to support multi-source data fusion; and finally, storing bounding box information and attribute information of elements in the fusion data into a relational database, and adopting Geocoding for indexing.
CN200910062052.0A 2009-05-14 2009-05-14 Three-dimensional digital earth-space data organizing and rendering method based on quad-tree index Expired - Fee Related CN101887595B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200910062052.0A CN101887595B (en) 2009-05-14 2009-05-14 Three-dimensional digital earth-space data organizing and rendering method based on quad-tree index

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200910062052.0A CN101887595B (en) 2009-05-14 2009-05-14 Three-dimensional digital earth-space data organizing and rendering method based on quad-tree index

Publications (2)

Publication Number Publication Date
CN101887595A CN101887595A (en) 2010-11-17
CN101887595B true CN101887595B (en) 2014-05-28

Family

ID=43073504

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200910062052.0A Expired - Fee Related CN101887595B (en) 2009-05-14 2009-05-14 Three-dimensional digital earth-space data organizing and rendering method based on quad-tree index

Country Status (1)

Country Link
CN (1) CN101887595B (en)

Families Citing this family (87)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9256981B2 (en) 2011-05-13 2016-02-09 Beijing Electric Power Economic Research Institute Method and device for processing geological information
CN103164417A (en) * 2011-12-12 2013-06-19 国家电网公司 Geological information processing method and device
CN102750737B (en) * 2011-08-31 2017-09-19 新奥特(北京)视频技术有限公司 A method and system for creating a three-dimensional earth model based on a quadtree
CN102346921A (en) * 2011-09-19 2012-02-08 广州市凡拓数码科技有限公司 Renderer-baking light mapping method of three-dimensional software
CN102509022B (en) * 2011-11-18 2015-02-18 武汉大学 Method for quickly constructing raster database facing to Virtual Earth
US9177412B2 (en) 2011-12-14 2015-11-03 Intel Corporation Techniques for multiple pass rendering
US9196081B2 (en) * 2011-12-15 2015-11-24 Intel Corporation Techniques for enhancing multiple view performance in a three dimensional pipeline
CN103186549A (en) * 2011-12-28 2013-07-03 泰瑞数创科技(北京)有限公司 Management method for mass tile raster data file
CN103208225B (en) * 2012-01-12 2015-10-28 中国科学院遥感应用研究所 A kind of tile map method for making and system
CN102663119B (en) * 2012-04-19 2014-06-11 北京天下图数据技术有限公司 Method for rendering massive vector data
CN102722549A (en) * 2012-05-28 2012-10-10 中国科学院遥感应用研究所 Cluster-based real-time rendering service of remote sensing data set
CN102999592A (en) * 2012-11-19 2013-03-27 北京中海新图科技有限公司 B/S architecture based global Argo multi-source marine data management and visualization system and method
CN103034999B (en) * 2012-12-05 2016-02-10 上海创图网络科技发展有限公司 A kind of pyramid dicing method based on dem data
CN103927933B (en) * 2013-01-14 2016-07-13 北京中交兴路信息科技有限公司 A kind of magnanimity moves method and the device that target renders
CN103390369B (en) * 2013-01-24 2016-04-13 中国科学院遥感应用研究所 A kind of map rendering intent and system
US9449410B2 (en) * 2013-04-11 2016-09-20 Qualcomm Incorporated Intra-frame timestamps for tile-based rendering
CN103345490B (en) * 2013-06-25 2016-08-31 中国资源卫星应用中心 The method showing first order satellite image in three-dimensional digital earth
CN103366006B (en) * 2013-07-19 2016-06-22 武汉睿数信息技术有限公司 A kind of method setting up geographical space index based on three-dimensional laser point cloud
CN103412962B (en) * 2013-09-04 2016-09-07 国家测绘地理信息局卫星测绘应用中心 A kind of storage method of magnanimity tile data and read method
CN103455624B (en) * 2013-09-16 2017-07-25 湖北文理学院 A Lightweight Global Multidimensional Remote Sensing Image Network Map Service Implementation Method
CN103546690A (en) * 2013-10-30 2014-01-29 天彩电子(深圳)有限公司 Method for obtaining and displaying motion data of action camera
CN103678631B (en) * 2013-12-19 2016-10-05 华为技术有限公司 page rendering method and device
CN103744998A (en) * 2014-01-23 2014-04-23 中国人民解放军国防科学技术大学 Hierarchical-division based storage and organization method for large-scale space vector data
CN103986763B (en) * 2014-05-13 2017-08-25 中国科学院地理科学与资源研究所 A kind of visual vector data streaming method of web oriented and system
CN104063900B (en) * 2014-05-30 2017-06-06 中国科学院遥感与数字地球研究所 The implementation of the three-dimensional earth based on webGL
CN104134187B (en) * 2014-07-30 2017-05-24 成都信息工程大学 Point-to-point tile map dynamic projection method and device
CN104239431B (en) * 2014-08-27 2017-11-14 广东威创视讯科技股份有限公司 Three-dimension GIS model display methods and device
CN104267940A (en) * 2014-09-17 2015-01-07 武汉狮图空间信息技术有限公司 Quick map tile generation method based on CPU+GPU
CN105589861B (en) * 2014-10-22 2019-02-05 杭州海康威视系统技术有限公司 A method of it is quickly positioned in map datum
CN104376050B (en) * 2014-10-31 2017-11-24 浙江科澜信息技术有限公司 It is a kind of in three-dimensional scenic dynamic dispatching and the method that shows massive vector data
CN104574489B (en) * 2014-12-16 2017-11-03 中国人民解放军理工大学 Landform and motion vector integrated approach based on lamination quaternary tree atlas
CN104408202B (en) * 2014-12-17 2017-08-08 武汉大学 The adaptive image pyramid dicing method of facing relation database extension
CN104615739B (en) * 2015-02-11 2018-10-09 华中科技大学 Data rapid filing method suitable for three-dimensional brain tissue high-resolution mass atlas
CN104766360B (en) * 2015-03-09 2017-12-22 腾讯科技(深圳)有限公司 The rendering intent and device of a kind of three-dimensional space data
CN104821007B (en) * 2015-04-15 2018-08-10 中国资源卫星应用中心 A kind of system for carrying out regarding display soon directly on the three-dimensional earth
CN104809217B (en) * 2015-04-28 2018-01-12 武大吉奥信息技术有限公司 A kind of GIS raster datas cloud storage method
CN104899282B (en) * 2015-06-02 2018-11-27 北京博阳世通信息技术有限公司 A kind of processing method and search method of multi-source multi-temporal satellite remote sensing tile data
CN105354206B (en) * 2015-06-15 2019-02-01 中国地质大学(武汉) Vector data tile based on big figure cutting cuts out figure accelerated method
CN105427380B (en) * 2015-12-23 2018-05-08 昆明能讯科技有限责任公司 One kind is based on detail three-dimensional map data processing method
CN105469434B (en) * 2015-12-31 2018-10-30 中南林业科技大学 A method of browsing large-scale remote sensing image in android system
US10319071B2 (en) * 2016-03-23 2019-06-11 Qualcomm Incorporated Truncated square pyramid geometry and frame packing structure for representing virtual reality video content
CN105931284B (en) * 2016-04-13 2019-12-31 中测新图(北京)遥感技术有限责任公司 Fusion method and device of three-dimensional texture TIN data and large scene data
CN105976426B (en) * 2016-04-27 2019-03-05 中国电子科技集团公司第二十八研究所 A kind of quick three-dimensional atural object model building method
CN106021480B (en) * 2016-05-19 2019-09-17 福建农林大学 A kind of parallel spatial division methods and its system based on grid dividing
CN106327577B (en) * 2016-08-17 2018-02-02 中国地质大学(武汉) Dimensional topography optimal method based on local curvature's entropy and quad-tree structure
CN106372262A (en) * 2016-10-25 2017-02-01 四川农业大学 System and method for urban outdoor public space urban home furnishing management
CN106780416B (en) * 2016-11-10 2020-04-07 国家电网公司 Method and device for planning path by using vector topographic map
CN106611056B (en) * 2016-12-27 2020-06-26 广州都市圈网络科技有限公司 Method and device for cutting three-dimensional simulation map
CN106844519B (en) * 2016-12-29 2019-07-26 中国科学院电子学研究所苏州研究院 Visualization method of remote sensing special products based on 3D digital earth
CN106844537B (en) * 2016-12-30 2020-11-17 重庆大学 Organization management method for space-time map tile data
CN106920273A (en) * 2017-02-21 2017-07-04 浙江科澜信息技术有限公司 A kind of global landform mist effect rendering intent of on-fixed concentration
CN106897425A (en) * 2017-02-24 2017-06-27 中国科学院电子学研究所 A kind of three-dimensional visualization method of terrestrial gravitation field data
CN106898047B (en) * 2017-02-24 2020-06-05 朱庆 Self-adaptive network visualization method for dynamic fusion of tilt model and multivariate model
CN107479898B (en) * 2017-08-28 2020-10-27 荆门程远电子科技有限公司 Three-dimensional virtual earth system architecture based on mobile platform
CN107862048B (en) * 2017-11-08 2021-05-18 四川易利数字城市科技有限公司 Map cache slice distributed cluster management method based on quad-tree hierarchical grid
CN108228723A (en) * 2017-12-11 2018-06-29 厦门亿力吉奥信息科技有限公司 WEB three-dimensional modeling datas construction method and computer readable storage medium
CN109408657B (en) * 2018-11-13 2021-02-02 国家基础地理信息中心 Method and system for rapidly drawing super-large-scale spatial data
CN109741436B (en) * 2018-12-04 2023-06-06 北京世纪高通科技有限公司 Method and device for rendering vector data
CN109829021A (en) * 2019-01-04 2019-05-31 广州市城市规划勘测设计研究院 A kind of map methods of exhibiting and device
CN109981346B (en) 2019-02-21 2020-07-10 清华大学 A method and device for creating a cyberspace coordinate system based on an autonomous system
CN110738735B (en) * 2019-10-23 2023-11-07 黄河勘测规划设计研究院有限公司 Method for improving display effect of three-dimensional digital earth platform
CN111210516B (en) * 2019-12-30 2023-04-18 成都赫尔墨斯科技股份有限公司 Software platform for integrated display control of avionics equipment
CN111221933A (en) * 2019-12-31 2020-06-02 武汉市珞珈俊德地信科技有限公司 Three-dimensional tile construction method for fusion of massive map data and building information model
CN111310089B (en) * 2020-02-17 2023-04-28 自然资源部第三地理信息制图院 Vector river network data online rapid loading and rendering method suitable for scale
CN111325841B (en) * 2020-02-18 2021-03-30 国网河北省电力有限公司 Smart map data system modeling method and application thereof in power audit operation
CN111581323B (en) * 2020-05-18 2024-07-12 中国科学院自动化研究所 Two/three-dimensional geographic space display method and system
CN111915711B (en) * 2020-08-04 2025-01-28 北京吉威空间信息股份有限公司 Method and device for obtaining stereoscopic images of three-level land classification spots supporting virtual VR
CN112017286A (en) * 2020-08-28 2020-12-01 北京国遥新天地信息技术有限公司 Seamless splicing display simulation method for digital earth skirt-free terrain tiles
CN112330791A (en) * 2020-09-18 2021-02-05 广州中望龙腾软件股份有限公司 Three-dimensional model creation method, intelligent terminal and storage device
CN111968212B (en) * 2020-09-24 2024-05-03 中国测绘科学研究院 Viewpoint-based three-dimensional urban scene data dynamic scheduling method
CN112598793B (en) * 2020-12-31 2023-02-17 华东计算技术研究所(中国电子科技集团公司第三十二研究所) Three-dimensional map tile vector grid integrated system and method
CN112902968B (en) * 2021-02-01 2024-10-22 成都紫瑞青云航空宇航技术有限公司 High-efficiency three-dimensional digital navigation earth generation method and system
CN112966041B (en) * 2021-02-02 2024-04-26 苍穹数码技术股份有限公司 Data processing method, device, equipment and storage medium
CN112700547B (en) * 2021-03-23 2021-06-29 北京简巨科技有限公司 Map making method and related equipment
CN112905830B (en) * 2021-03-31 2024-08-16 自然资源部国土卫星遥感应用中心 Satellite remote sensing image acquisition method and related device
CN113096246B (en) * 2021-05-12 2022-04-12 泰瑞数创科技(北京)有限公司 Three-dimensional map data updating method
CN114092575B (en) * 2021-11-24 2022-04-12 北京清晨动力科技有限公司 Digital earth real-time coloring method and device
CN114490899B (en) * 2021-12-27 2022-09-23 北京吉威数源信息技术有限公司 Multi-layer space big data superposition and combination method, device, equipment and storage medium
CN114399405B (en) * 2022-01-20 2025-06-27 黄河勘测规划设计研究院有限公司 Large-scale water conservancy project management method and system based on GIS and BIM
CN114463505B (en) * 2022-02-15 2023-01-31 中国人民解放军战略支援部队航天工程大学士官学校 Outer space environment element model construction method and system
CN114882181B (en) * 2022-04-14 2025-06-13 江苏徐工工程机械研究院有限公司 A method and system for generating unmanned three-dimensional maps of open-pit mines
CN115170779B (en) * 2022-07-18 2023-03-28 生态环境部环境规划院 Remote sensing image rendering and network publishing method
CN116310093B (en) * 2023-02-23 2023-09-15 南京麦堤微林信息科技有限公司 Virtual three-dimensional urban geographic scene sand table model construction system and method thereof
CN116630566B (en) * 2023-07-24 2024-01-23 中国电子科技集团公司第十五研究所 Geometric modeling and rendering method of digital-earth-oriented ellipsoidal BLH slice
CN118674888A (en) * 2024-08-22 2024-09-20 慧航(江西)数字科技有限公司 Three-dimensional visual management method and system for underground space
CN118691487B (en) * 2024-08-27 2024-12-13 国网江苏省电力有限公司南京供电分公司 High-low precision Digital Elevation Model (DEM) fusion method and system for three-dimensional digital power grid
CN119648927B (en) * 2025-02-19 2025-05-02 杭州市勘测设计研究院有限公司 High-point analysis method and system based on double spatial indexes

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1959670A (en) * 2005-11-02 2007-05-09 北京灵图软件技术有限公司 Method for storing topographic image data in pyramid architecture
CN101021868A (en) * 2007-03-06 2007-08-22 武汉大学 Terrain data storing method based on object storage
CN101082934A (en) * 2007-06-19 2007-12-05 北京师范大学 Three-dimensional visual method of very-large-scale space data in network surroundings

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1959670A (en) * 2005-11-02 2007-05-09 北京灵图软件技术有限公司 Method for storing topographic image data in pyramid architecture
CN101021868A (en) * 2007-03-06 2007-08-22 武汉大学 Terrain data storing method based on object storage
CN101082934A (en) * 2007-06-19 2007-12-05 北京师范大学 Three-dimensional visual method of very-large-scale space data in network surroundings

Also Published As

Publication number Publication date
CN101887595A (en) 2010-11-17

Similar Documents

Publication Publication Date Title
CN101887595B (en) Three-dimensional digital earth-space data organizing and rendering method based on quad-tree index
CN108919944B (en) Virtual roaming method for realizing data lossless interaction at display terminal based on digital city model
CN111784833B (en) Flood evolution situation three-dimensional dynamic visualization display method based on WebGL
US8665266B2 (en) Global visualization process terrain database builder
CN103455624B (en) A Lightweight Global Multidimensional Remote Sensing Image Network Map Service Implementation Method
CN104778744B (en) Extensive three-dimensional forest Visual Scene method for building up based on Lidar data
CN105976426B (en) A kind of quick three-dimensional atural object model building method
Shan et al. Auxiliary use and detail optimization of computer VR technology in landscape design
Liang et al. Embedding user-generated content into oblique airborne photogrammetry-based 3D city model
CN115187647A (en) Vector-based road three-dimensional live-action structured modeling method
CN108921943A (en) A kind of road threedimensional model modeling method based on lane grade high-precision map
CN112530009A (en) Three-dimensional topographic map drawing method and system
CN103971414A (en) Method and system for making visualized true three-dimensional map
CN115859414A (en) Cross-coordinate system use method for global scale geographic information base map
Wu et al. [retracted] intelligent city 3d modeling model based on multisource data point cloud algorithm
Zhu et al. CyberCity GIS (CCGIS): integration of DEMs, images, and 3D models
CN114490907A (en) Method and device for constructing famous city management database and storage medium
CN106683183A (en) Online three-dimensional visualization method for construction project model
Kumar et al. A CityGML extension for handling very large TINs
Bi et al. Research on CIM basic platform construction
CN115409962B (en) Method for constructing coordinate system in illusion engine, electronic device and storage medium
CN113870409B (en) Three-dimensional graph lightweight method based on industrial simulation scene
CN116630830A (en) Method and system for presenting full-space information flow of power pipeline based on oblique photography
CN115409961A (en) Thermodynamic diagram construction method based on illusion engine, electronic equipment and storage medium
Song et al. Research on 3D campus integrated management based on ArcGIS Pro and CityEngine

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140528

Termination date: 20160514

CF01 Termination of patent right due to non-payment of annual fee