CN109086286B - Method for producing and publishing color topographic map - Google Patents
Method for producing and publishing color topographic map Download PDFInfo
- Publication number
- CN109086286B CN109086286B CN201711283552.8A CN201711283552A CN109086286B CN 109086286 B CN109086286 B CN 109086286B CN 201711283552 A CN201711283552 A CN 201711283552A CN 109086286 B CN109086286 B CN 109086286B
- Authority
- CN
- China
- Prior art keywords
- data
- geospatial
- elevation
- wms
- color
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/50—Lighting effects
- G06T15/60—Shadow generation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2200/00—Indexing scheme for image data processing or generation, in general
- G06T2200/32—Indexing scheme for image data processing or generation, in general involving image mosaicing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Remote Sensing (AREA)
- Processing Or Creating Images (AREA)
- Instructional Devices (AREA)
Abstract
The invention discloses a method for producing and publishing a color topographic map, which is based on the integration and support of related geographic information technology, designs and provides sub-processes of geographic space network resource storage management, digital elevation data processing, batch geographic space data rendering, standard geographic space network resource service publishing, standard geographic space network resource caching and the like, realizes the hierarchical block rapid production of topographic map data through related processing processes, and provides a standard method for publishing and accessing the online topographic map. The invention is mainly used for integrating geographic space data service resources, designing related data processing flows by combining data characteristics of digital orthographic images, digital elevations, vector maps and the like, generating clear and accurate topographic map data and services and improving the production and access efficiency of topographic data.
Description
Technical Field
The invention belongs to the technical field of geographic information, and mainly relates to a method for producing and publishing a color topographic map.
Background
In the current society of the rapid development of the internet industry and the information technology, the related industries of the geographic information technology develop rapidly, wherein various network map services gradually become an important support for the geographic information technology and the large data platform technology. The network map service based on the internet resources is an important component of a public basic geographic information platform, and the borne basic geographic information is rich and diverse, so that people can more intuitively know the geographic space. The topographic map data is used as an important component in the geographic information network map service, and is a multidimensional geographic information network service which is issued after vector geographic space data, grid geographic space data, digital elevation and the like are organically combined, so that the topographic distribution condition in a certain area can be visually known through height data models such as colors, mountain shadows and the like, and meanwhile, the accurate analysis and measurement of the geographic space data information in the area can be clearly carried out according to the element marking condition of the vector map. With the arrival of the big data era of the geographic information industry, the topographic map data and services of various natural elements and social elements have increasingly obvious effects on city planning, disaster assessment, national defense construction and the like. As the topographic map data has higher requirement on the integration level of various geographic spatial elements, the early topographic map is difficult to draw and produce, with the continuous development of remote sensing technology and the perfection of computer mapping software functions, remote sensing provides a data source for the production of the topographic map, a geographic information system is used as a comprehensive technical system for collecting, analyzing and processing spatial information based on computer mapping, and provides a means for acquiring information and a method for directly and quickly drawing the map, and the production and production technology of the topographic map is greatly stepped forward by the combination of the remote sensing data and the geographic information system. However, many problems still face today in terms of topographic map data source, map accuracy, data quality, service efficiency, etc.: the network topographic map service based on geographic information transmission and space entity expression purposes has the problems of different forms in different fields of different industries, serious loss of space and time semantics, incomplete expression objects, inaccurate scale, various bearing file forms, huge and complex production systems and the like, and meanwhile, the high-precision network topographic map service production is closely related to a remote sensing technology, a digital image processing technology and a big data parallel processing technology, and also needs support and cooperation of different industries such as surveying and mapping, urban construction, aerospace and the like, so that the early single topographic map data production mode cannot meet the existing comprehensive requirements, and the efficient, accurate and convenient topographic map service mode is very important. According to the problems, the existing geographic information related technology and network geographic information resources are integrated, a process-based network topographic map production and service release process is designed, the collection and visualization of multi-source geographic information resources are realized, the data production efficiency and flexibility are improved, the resource sharing and pushing capacity is enhanced by combining part of network technologies, and the verification is applied to related projects.
Disclosure of Invention
The invention aims to solve the technical problem of a standard topographic map rapid production and organization method, which completes the production and rapid access of a color topographic map facing to a geographic space network resource and realizes the comprehensive utilization of different types of geographic space data such as remote sensing images, digital elevations, vector maps and the like.
The invention adopts the following technical scheme for solving the technical problems:
a method for producing and publishing a color topographic map specifically comprises the following steps:
a, geographic space data storage management;
b, processing geospatial elevation data;
c, rendering the geospatial data in batches;
d, standard geospatial network resource service release;
e-standard geospatial network resource service caching.
As a further preferred scheme of the color topographic map production and release method, the geographic space data storage management method A specifically comprises the following steps:
step a1, geospatial data import:
step A11, the geospatial vector data and the grid data acquired from the network are sent to a geospatial database, and data projection information is defined according to different data sources when the geospatial vector data and the grid data are imported;
step A12, for data of different coordinate systems, projection conversion is carried out during importing, and the data are converted into a WGS84 coordinate system in a unified manner, so that subsequent processing is facilitated;
step A2, correcting expired or error data in the database;
a21, selecting a data item which needs to meet a certain condition according to the structured attribute information of the geospatial data and the spatial query operation of the Geometry object;
step A22, executing data updating operation, setting the attribute information of the inquired data item or the Geometry object representing the geospatial position of the data for correction;
step A3, supplementing the missing data information by means of insertion;
step A31, reading data information of the island of the south China sea island, such as name, terrain type and geospatial position information represented by WKT format, from the obtained other data sources;
step A32, performing database insertion operation, and inserting data items of the island of the south China sea island into a database in batches;
step A33, packaging the geographic space position information into a Geometry object in the insertion process, and performing coordinate system conversion if necessary;
step A4, performing complex and simple conversion on place names of data in the Taiwan region of China and realizing the ten thousand Chinese conversion of the global place names;
step A41, obtaining the name attribute of the data item in the Taiwan range through space query operation, and exporting the name attribute to csv file; step A42, establishing a place name data dictionary for complex and simple conversion according to the exported csv file, and importing the place name data dictionary into a database as a new data table;
step A43, the place name information of the data item in step A41 is replaced according to the imported place name data dictionary through fuzzy inquiry, and a field is added in the original data table for storing the original place name information;
as a further preferable scheme of the method for producing and publishing the color topographic map, the method for processing the geospatial elevation data comprises the following specific steps:
step B1, splicing elevation data;
step B11, writing the paths of the elevation data files divided according to the longitude and latitude into a spliced file list;
step B12, appointing a series of splicing parameters such as the maximum and minimum resolution ratio of the splicing operation processing, the resolution ratio in the x and y directions, the maximum and minimum xy ranges, the file list to be spliced, the name and the path of the target raster file, and the nodata value in the appointed splicing process;
step B13, executing splicing operation to generate virtual raster data or other data of data type that can be supported;
step B2, resampling the image:
step B21, determining the resolution value of the elevation data required under different magnification levels according to the corresponding relation between the zoom level and the resolution of the map service;
step B22, the formula can also be expressed by latitude range and latitude direction parameter value, the calculation result is consistent;
step B23, specifying the resolution and the resampling mode of the target elevation data file to be generated, wherein the resampling mode can select nearest neighbor sampling or bilinear interpolation sampling and the like;
step B24, resampling the elevation data to obtain elevation data with specified resolution;
and step B3, cutting the elevation data:
step B31, opening the vector file through a visualization tool, selecting an area needing to be reserved, and storing the area as a new layer; or drawing the area needing to be reserved by a drawing tool and storing the area as a layer;
step B32, if a plurality of discontinuous areas exist, storing the discontinuous areas as a plurality of layers respectively, then combining the layers, and outputting the combined layers as a shp file;
step B33, using a grid processing tool to perform cutting, inputting grid data as data to be cut, outputting a vector layer with a range for cutting, and setting a value of null data in cutting parameters;
and step B34, exporting the layers obtained after cutting into a raster data file.
Step B4, generating a raster image of the color topographic map:
step B41, checking basic information of the elevation data, calculating the range of the elevation value of the data and acquiring information such as a null data value;
step B42, determining color scheme of the color topographic map according to the range of the elevation value and the required display effect, and generating a color configuration file;
step B43, generating color terrain raster image data according to the designated color configuration file;
step B5, generating a mountain shadow map:
step B51, a series of parameter values such as vertical exaggeration coefficient and unit, light azimuth angle, light elevation angle and the like in the processing process are designated according to requirements;
step B52, calculating the mountain shadow value, the calculation formula is: hillshade × ((cos (Zneith) × cos (slope)) + (sin (Zneith) × sin (slope) × cos (Azimuth-Aspect))), wherein Zneith is the solar apex angle, Azimuth is the solar elevation angle, Slo is the solar elevation anglepe is the slope and Aspect is the slope. All angles take radian as a calculation unit;
in step B53, the gradient is calculated by the formula: slope ═ arctan ((h)2-h1)/(x2-x1) Wherein h) is2-h1Representing the difference in elevation, x, between two adjacent points2-x1Represents the horizontal distance between two adjacent points;
step B54, generating mountain shadow raster image data according to the calculation result;
step B6, generating a slope shadow map:
step B61, establishing a gradient grid image according to the elevation data, and designating a vertical unit in the establishing process;
step B62, establishing a color profile of the slope shadow image;
and step B63, generating a rendered gradient shadow raster image according to the generated gradient raster image and the color configuration file.
Step B7, generating contour data:
step B71, determining equal height distance according to the zooming level and the required display effect;
step B72, the wave band serial number used for generating the contour line is appointed, and parameters such as the field name of the elevation value, the contour distance information, the output contour line vector file type, the layer name and the like are saved;
step B73, a contour vector data file is generated.
As a further preferable scheme of the method for producing and publishing the color topographic map, the step C of rendering the geospatial data in batches comprises the following steps:
step C1, configuring the spatial reference frame of the Mapnik rendered data.
And C2, configuring the rendering style of the vector data and the raster data, and specifying the style application level and style parameter setting of point, line, surface and other labels.
And step C3, configuring the style name and the data source referred by the vector data layer.
And step C4, adjusting the position relation of the vector data layer and the raster data layer, and combining into a graph layer group.
Step C5, a rendering process is performed.
As a further preferable scheme of the color topographic map production and release method, the specific implementation steps of the standard D geospatial network resource service release are as follows:
step D1, preparing geospatial data and Mapnik rendering style files;
step D2, starting the WMS server, calling the Mapnik rendering engine at the bottom layer, and issuing WMS service;
and D3, monitoring the WMS request from the client and returning the corresponding network map service data to the client.
As a further preferred scheme of the color topographic map production and release method, the E-standard geographic space network resource service caching specifically comprises the following steps:
and step E1, configuring a cache path of the WMS cache server.
Step E2, configuring WMS layer parameters of the WMS cache server, such as a coordinate system, a range, a tile format, a request WMS layer group name and a WMS service request address;
step E3, starting a WMS cache server, and pre-caching the WMS layer configured by the cache server;
step E4, initiating a WMS request to the WMS server according to the pre-caching setting, and caching the requested network map service data to the caching server;
and E5, monitoring the WMS request from the client and returning the cached network map service data to the client.
Compared with the prior art, the invention adopting the technical scheme has the following technical effects:
1) the integration and visualization of multi-source geographic information resources are realized through a network topographic map production and service release mode based on a process;
2) comprehensive utilization of various different types of geospatial data maps such as remote sensing images, digital elevations, vector maps and the like is realized;
3) compared with the prior art, the method has the advantages of convenient data acquisition, wide application scene, high production data efficiency, flexible style and the like;
4) the whole data production process can be realized on different platforms, and the cross-platform performance is good. 5) The method has higher performance for the application in the high concurrency scene, and can provide good user experience.
Drawings
FIG. 1 is a block diagram of a method for producing and publishing a color terrain map based on geospatial network resources;
FIG. 2 is a flow diagram of an implementation of a geospatial network resource-based color terrain map production and distribution;
FIG. 3(a) is an overall effect diagram of color terrain map production and distribution based on geospatial network resources;
FIG. 3(b) is a partial effect of a geospatial network resource based color topography production and release with a resolution of 1200 meters;
FIG. 3(c) is a partial effect of a geospatial network resource based color topography production and release at a resolution of 160 meters;
FIG. 3(d) is a partial effect of a geospatial network resource based color topography production and release at a resolution of 80 meters;
FIG. 4(a) is a comparison of the effect of a geospatial network resource based color topographic map and a Google topographic map at a resolution of 1200 meters;
FIG. 4(b) is a comparison of the effect of a geospatial network resource based color topographic map and a Google topographic map at a resolution of 600 meters;
FIG. 4(c) is a comparison of the effect of a geospatial network resource based color topographic map and a Google topographic map at a resolution of 160 meters;
FIG. 4(d) is a comparison of the effect of a geospatial network resource-based colormap versus Google at a resolution of 80 meters
Fig. 4(e) is a comparison of the effect of a geospatial network resource-based color topographic map and google topographic map at a resolution of 40 meters.
Detailed Description
The technical scheme of the invention is further explained in detail by combining the attached drawings:
according to the invention, a set of complete production and release schemes of the color topographic map based on the geospatial network resources are realized by integrating a geospatial data storage technology, an elevation data processing technology, a batch geospatial data rendering technology, a standard geospatial network resource service release and cache technology and the like, as shown in FIG. 1. The following describes in detail an embodiment of the present invention with reference to fig. 2.
The implementation steps of the geospatial data storage management sub-process are as follows:
geospatial data import
And (3) the geospatial vector data and the raster data acquired from the network are put into a geospatial database, and the information is projected according to different definition data of data sources when the data is imported.
The data of different coordinate systems are subjected to projection conversion in the importing process and are converted into a WGS84 coordinate system in a unified mode, and the subsequent processing is facilitated.
Corrections are made for expired or erroneous data in the database.
And selecting data items which need to meet a certain condition according to the structured attribute information (such as administrative level information) of the geospatial data and by utilizing the spatial query operation on the Geometry object.
And executing data updating operation, and setting attribute information (such as name) of the inquired data item or a Geometry object representing the geographic spatial position of the data for correction.
The missing data information is supplemented by means of insertion. The following description will be given in detail by taking the example of data insertion in the island of south China sea.
And reading data information of the island of the south China sea island, such as name, terrain type and geospatial position information expressed in WKT format, from the obtained other data sources.
And performing database insertion operation, and inserting the data items of the island of south China sea into the database in batches.
And encapsulating the geographic space position information into a Geometry object in the insertion process, and performing coordinate system conversion if necessary.
The place name of data in the Taiwan area of China is subjected to traditional and simple conversion and the ten-thousand-language Chinese realization of the global place name (the desert name is used as an example). The detailed procedure of this step is described below.
And acquiring the name attribute of the data item in the range of Taiwan by space query operation, and exporting the name attribute as a csv file.
And establishing a place name data dictionary for the traditional and simple conversion according to the exported csv file, and importing the place name data dictionary into a database as a new data table.
Replacing the place name information of the data items in the step 1) according to the imported place name data dictionary through fuzzy query, and adding a field for storing original place name information in the original data table.
The process of Handization of the names of the ten thousand languages is similar to the above steps.
The geospatial elevation data processing sub-process comprises the following implementation steps:
and splicing the elevation data.
And writing the paths of the elevation data files divided according to the longitude and latitude into a spliced file list.
And a series of splicing parameters such as the maximum and minimum resolution, the resolution in the x and y directions, the maximum and minimum x and y ranges, a list of files to be spliced, the name and path of a target raster file, and a nondata value in the splicing process are specified.
And executing splicing operation to generate virtual raster data or other data of data types which can be supported.
Resampling an image
And determining the value of the resolution of the elevation data required under different magnification levels according to the corresponding relation between the zoom level and the resolution of the map service.
The resolution is calculated as follows: res ═ (Lon)max-Lonmin)/(CountTile×CountPixel) Wherein LonmaxAnd LonminRespectively representing the maximum and minimum values of the longitude range, CountTileRepresents the total number of tiles in the longitudinal direction of the magnification level, CountPixelThe total number of pixels representing the horizontal or vertical direction of a single map tile may be selected 256, 512 or other values according to the settings of the map service to be published.
The formula can also be expressed by latitude range and latitude direction parameter values, and the calculation results are consistent.
The resolution and the resampling mode of the target elevation data file to be generated are specified, and the resampling mode can select the nearest neighbor sample or bilinear interpolation sample, etc.
And resampling the elevation data to obtain the elevation data with the specified resolution.
And cutting the elevation data.
Opening the vector file through a visualization tool, selecting an area needing to be reserved, and storing the area as a new layer; or drawing the area needing to be reserved through a drawing tool and saving the area as the layer.
If a plurality of discontinuous areas exist, the discontinuous areas can be respectively stored as a plurality of layers, then layer combination is carried out, and the combined layers are output as an shp file.
And (4) using a grid processing tool to perform cutting, wherein input grid data are data to be cut, and an output range is a vector layer for cutting. The value of the null data is set in the clipping parameter.
And exporting the layers obtained after cutting into a raster data file.
Generating raster images of a color topographic map
And checking basic information of the elevation data, calculating an elevation value range of the data, and acquiring information such as a null data value.
And determining a color scheme of the color topographic map according to the range of the elevation values and the required display effect, and generating a color configuration file (comprising RGB values and a transparency parameter alpha).
And generating color terrain raster image data according to the specified color configuration file.
Generating mountain shadow maps
The vertical exaggeration coefficient and unit, the light azimuth angle, the light elevation angle, and a series of parameter values in the processing process are specified according to requirements.
Calculating the mountain shadow value by the following formula: hillshade × ((cos (Zneith) × cos (Slope)) + (sin (Zneith) × sin (Slope) × cos (Azimuth-Aspect)), where Zneith is the solar apex angle, Azimuth is the solar elevation angle, Slope is the Slope, and Aspect is the Slope direction.
The calculation formula of the gradient is as follows: slope ═ arctan ((h)2-h1)/(x2-x1) Wherein h) is2-h1Representing the difference in elevation, x, between two adjacent points2-x1Representing the horizontal distance between two adjacent points.
And generating mountain shadow raster image data according to the calculation result.
Generating grade shadow maps
And establishing a gradient raster image according to the elevation data, and designating a vertical unit in the establishing process.
The gradient calculation formula is as shown above.
A color profile of the gradient shadow image is established.
And generating a rendered gradient shadow raster image according to the generated gradient raster image and the color configuration file.
Generating contour data
Determining equal-height distances according to the zooming levels and the required display effect;
and the serial number of the wave band for generating the contour line is specified, and parameters such as the field name for storing the elevation value, the contour distance information, the output contour line vector file type, the layer name and the like are stored.
And generating a contour vector data file.
The implementation steps of the batch geospatial data rendering sub-process are as follows:
and configuring a spatial reference frame of the Mapnik rendering data.
And configuring rendering styles of vector data and raster data, and specifying style application levels and style parameter settings of marks such as points, lines and surfaces.
And configuring style names and data sources referenced by the vector data layers.
And adjusting the position relation of the vector data layer and the raster data layer to combine a graph layer group.
A rendering process is performed.
The standard geospatial network resource service publishing sub-process comprises the following specific implementation steps:
geospatial data and Mapnik rendering style files are prepared.
And starting the WMS server, calling the Mapnik rendering engine at the bottom layer, and issuing the WMS service.
And monitoring a WMS request from the client, and returning corresponding network map service data to the client.
The standard geospatial network resource service caching subprocess comprises the following implementation steps:
and configuring a cache path of the WMS cache server.
And configuring WMS layer parameters of the WMS cache server, such as a coordinate system, a range, a tile format, a request WMS layer group name and a WMS service request address.
And starting the WMS cache server, and pre-caching the WMS layer configured by the cache server.
And initiating a WMS request to the WMS server according to the pre-caching setting, and caching the requested network map service data to a cache server.
And monitoring a WMS request from the client, and returning the cached network map service data to the client.
The color terrain preview effect after release is shown in fig. 3. FIG. 3(a) is an overall effect diagram of color terrain map production and distribution based on geospatial network resources; FIG. 3(b) is a partial effect of a geospatial network resource based color topography production and release with a resolution of 1200 meters; FIG. 3(c) is a partial effect of a geospatial network resource based color topography production and release at a resolution of 160 meters; FIG. 3(d) is a partial effect of the resolution of 80 meters for color terrain production and distribution based on geospatial network resources.
As shown in fig. 4, for the same area, by comparing the google topographic map (left) with the color topographic map (right) issued by the present invention, the color topographic map has the characteristics of richer geographic elements, more obvious ground object contrast, clearer visualization degree, etc. Wherein, fig. 4(a) is the comparison of the effect of the color topographic map based on the geospatial network resources and the google topographic map when the resolution is 1200 meters; FIG. 4(b) is a comparison of the effect of a geospatial network resource based color topographic map and a Google topographic map at a resolution of 600 meters; FIG. 4(c) is a comparison of the effect of a geospatial network resource based color topographic map and a Google topographic map at a resolution of 160 meters; FIG. 4(d) is a comparison of the effect of a geospatial network resource based color topographic map and a Google topographic map at a resolution of 80 meters; 4(e) comparing the effect of the color topographic map based on the geospatial network resources with the google topographic map at a resolution of 40 meters.
Claims (5)
1. A method for producing and publishing a color topographic map is characterized by comprising the following steps: the method specifically comprises the following steps:
a, geographic space data storage management;
b, processing geospatial elevation data;
c, rendering the geospatial data in batches;
d, standard geospatial network resource service release;
e, standard geographic space network resource service caching;
b, processing the geospatial elevation data specifically comprises the following steps:
step B1, splicing elevation data;
step B11, writing the paths of the elevation data files divided according to the longitude and latitude into a spliced file list;
step B12, appointing the maximum and minimum resolution ratio of the splicing operation, the resolution ratio in the x and y directions, the maximum and minimum x and y ranges, the list of the files to be spliced, the name and the path of the target raster file, and the nondata value in the splicing process;
step B13, executing splicing operation to generate virtual raster data or other data of data type that can be supported;
step B2, resampling the image:
step B21, determining the resolution value of the elevation data required under different magnification levels according to the corresponding relation between the zoom level and the resolution of the map service;
step B22, or expressing the values of the resolution of the elevation data required under different amplification levels by using the parameter values in the latitude range and the latitude direction, and the calculation results are consistent;
step B23, specifying the resolution and the resampling mode of the target elevation data file to be generated, wherein the resampling mode selects the nearest sampling or bilinear interpolation sampling;
step B24, resampling the elevation data to obtain elevation data with specified resolution;
and step B3, cutting the elevation data:
step B31, opening the vector file through a visualization tool, selecting an area needing to be reserved, and storing the area as a new layer; or drawing the area needing to be reserved by a drawing tool and storing the area as a layer;
step B32, if a plurality of discontinuous areas exist, storing the areas as a plurality of layers respectively, then combining the layers, and outputting the combined layers as an shp file;
step B33, using a grid processing tool to perform cutting, inputting grid data as data to be cut, outputting a vector layer with a range for cutting, and setting a value of null data in cutting parameters;
step B34, exporting the layers obtained after cutting into a raster data file;
step B4, generating a raster image of the color topographic map:
step B41, checking the basic information of the elevation data, calculating the range of the elevation value of the data and acquiring the null data value;
step B42, determining color scheme of the color topographic map according to the range of the elevation value and the required display effect, and generating a color configuration file;
step B43, generating color terrain raster image data according to the designated color configuration file;
step B5, generating a mountain shadow map:
step B51, specifying the vertical exaggeration coefficient and unit, the light azimuth angle and the light elevation angle in the processing process according to requirements;
step B52, calculating the mountain shadow value, the calculation formula is: hillshade × (cos (Zneith) × cos (Slope) (sin) (Zneith) × sin (Slope) × cos (Azimuth-Aspect))), wherein Zneith is the solar apex angle, Azimuth is the solar elevation angle, Slope is the Slope, and Aspect is the Slope direction, and all angles take the radian as a calculation unit;
in step B53, the gradient is calculated by the formula: slope ═ arctan ((h)2-h1)/(x2-x1) Wherein h) is2-h1Representing the difference in elevation, x, between two adjacent points2-x1Represents the horizontal distance between two adjacent points;
step B54, generating mountain shadow raster image data according to the calculation result;
step B6, generating a slope shadow map:
step B61, establishing a gradient grid image according to the elevation data, and designating a vertical unit in the establishing process;
step B62, establishing a color profile of the slope shadow image;
step B63, generating a rendered gradient shadow raster image according to the generated gradient raster image and the color configuration file;
step B7, generating contour data:
step B71, determining equal height distance according to the zooming level and the required display effect;
step B72, appointing the wave band serial number for generating contour line, and storing the field name of elevation value, the contour distance information, the outputted contour line vector file type and the layer name;
step B73, a contour vector data file is generated.
2. A method for producing and distributing a color map according to claim 1, characterized in that: the method comprises the following specific steps of A, geospatial data storage management:
step a1, geospatial data import:
step A11, the geospatial vector data and the grid data acquired from the network are sent to a geospatial database, and data projection information is defined according to different data sources when the geospatial vector data and the grid data are imported;
step A12, for data of different coordinate systems, projection conversion is carried out during importing, and the data are converted into a WGS84 coordinate system in a unified manner, so that subsequent processing is facilitated;
step A2, correcting expired or error data in the database;
a21, selecting a data item which needs to meet a certain condition according to the structured attribute information of the geospatial data and the spatial query operation of the Geometry object;
step A22, executing data updating operation, setting the attribute information of the inquired data item or the Geometry object representing the geospatial position of the data for correction;
step A3, supplementing the missing data information by means of insertion;
step A31, reading data information of the island of the south China sea island from other acquired data sources, wherein the data information comprises names, terrain types and geospatial position information represented in a WKT format;
step A32, performing database insertion operation, and inserting data items of the island of the south China sea island into a database in batches;
step A33, packaging the geographic space position information into a Geometry object in the insertion process, and performing coordinate system conversion;
step A4, performing complex and simple conversion on place names of data in the Taiwan region of China and realizing the ten thousand Chinese conversion of the global place names;
step A41, obtaining the name attribute of the data item in the Taiwan range through space query operation, and exporting the name attribute to csv file; step A42, establishing a place name data dictionary for complex and simple conversion according to the exported csv file, and importing the place name data dictionary into a database as a new data table;
and step A43, replacing the place name information of the data item in the step A41 according to the imported place name data dictionary through fuzzy query, and adding a field for storing the original place name information in the original data table.
3. A method for producing and distributing a color map according to claim 1, characterized in that: c, rendering the geospatial data in batch, which comprises the following specific steps:
step C1, configuring a spatial reference system of Mapnik rendering data;
step C2, configuring the rendering style of the vector data and the raster data, and specifying the style application level and the style parameter setting of points, lines and surfaces;
step C3, configuring style names and data sources referenced by the vector data layers;
step C4, adjusting the position relation between the vector data layer and the raster data layer, and combining into a graph layer group;
step C5, a rendering process is performed.
4. A method for producing and distributing a color map according to claim 1, characterized in that: the specific implementation steps of the standard geospatial network resource service release are as follows:
step D1, preparing geospatial data and Mapnik rendering style files;
step D2, starting the WMS server, calling the Mapnik rendering engine at the bottom layer, and issuing WMS service;
and D3, monitoring the WMS request from the client and returning the corresponding network map service data to the client.
5. A method for producing and distributing a color map according to claim 1, characterized in that: the E-standard geographic space network resource service caching method specifically comprises the following steps:
step E1, configuring a cache path of the WMS cache server;
step E2, configuring WMS layer parameters of the WMS cache server, wherein the WMS layer parameters comprise a coordinate system, a range, a tile format, a request WMS layer group name and a WMS service request address;
step E3, starting a WMS cache server, and pre-caching the WMS layer configured by the cache server;
step E4, initiating a WMS request to the WMS server according to the pre-caching setting, and caching the requested network map service data to the caching server;
and E5, monitoring the WMS request from the client and returning the cached network map service data to the client.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711283552.8A CN109086286B (en) | 2017-12-07 | 2017-12-07 | Method for producing and publishing color topographic map |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711283552.8A CN109086286B (en) | 2017-12-07 | 2017-12-07 | Method for producing and publishing color topographic map |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109086286A CN109086286A (en) | 2018-12-25 |
| CN109086286B true CN109086286B (en) | 2021-06-29 |
Family
ID=64839552
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711283552.8A Active CN109086286B (en) | 2017-12-07 | 2017-12-07 | Method for producing and publishing color topographic map |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109086286B (en) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109947882A (en) * | 2019-02-26 | 2019-06-28 | 北京悦图遥感科技发展有限公司 | Map Services dissemination method and device |
| CN109949416A (en) * | 2019-02-28 | 2019-06-28 | 华南机械制造有限公司 | Network Three-dimensional GIS-Geographic Information System based on WebGL technology |
| CN109977510B (en) * | 2019-03-15 | 2022-12-06 | 长江水利委员会长江科学院 | Hydrological model network publishing method based on GIS technology |
| CN112052864A (en) * | 2020-10-20 | 2020-12-08 | 网易（杭州）网络有限公司 | Image drawing method and device, electronic equipment and readable storage medium |
| CN112700547B (en) * | 2021-03-23 | 2021-06-29 | 北京简巨科技有限公司 | Map making method and related equipment |
| CN114018243A (en) * | 2021-11-05 | 2022-02-08 | 苍穹数码技术股份有限公司 | Map data processing method, device, equipment and storage medium |
| CN114820967B (en) * | 2022-04-28 | 2022-09-27 | 四川见山科技有限责任公司 | Layered image data real-time drawing method and system |
| CN115809244A (en) * | 2022-11-20 | 2023-03-17 | 北京泰策科技有限公司 | Massive spatial data access, processing and release method |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101593209A (en) * | 2008-12-30 | 2009-12-02 | 中国科学院地理科学与资源研究所 | A kind of alternating type online statistical map generation and distributing device and method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE112009004442B4 (en) * | 2009-03-04 | 2014-11-06 | Mitsubishi Electric Corporation | Card distribution server and card information terminal |
-
2017
- 2017-12-07 CN CN201711283552.8A patent/CN109086286B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101593209A (en) * | 2008-12-30 | 2009-12-02 | 中国科学院地理科学与资源研究所 | A kind of alternating type online statistical map generation and distributing device and method |
Non-Patent Citations (1)
| Title |
|---|
| GeoMedia Web Map地图生成与发布技术;王晓军等;《第六届全国计算机应用联合学术会议论文集》;20021001;7-274至7-280 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109086286A (en) | 2018-12-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109086286B (en) | Method for producing and publishing color topographic map | |
| US11222465B2 (en) | Embedded urban design scene emulation method and system | |
| CN112115198B (en) | Urban remote sensing intelligent service platform | |
| CN112559667B (en) | Map editor, map engine, map editing system, map editing method, and storage medium | |
| CN110321443B (en) | Three-dimensional live-action model database construction method and device and data service system | |
| US20230245389A1 (en) | Generating three-dimensional geo-registered maps from image data | |
| CN103606188A (en) | Geographical information on-demand acquisition method based on image point cloud | |
| CN110660125B (en) | Three-dimensional modeling device for power distribution network system | |
| CN112113544B (en) | Remote sensing data processing method and system based on unmanned aerial vehicle image | |
| US10067950B1 (en) | Systems and methods for efficiently organizing map styling information | |
| CN112100297A (en) | Three-dimensional visualization method for data management in urban updating | |
| CN115471619A (en) | City three-dimensional model construction method based on stereo imaging high-resolution satellite image | |
| JP2010506337A (en) | System and method for visualizing and measuring real-world 3-D spatial data | |
| CN111986320B (en) | Smart city application-oriented DEM and oblique photography model space fitting optimization method | |
| CN115357675B (en) | Method and system for establishing image control point database through standardized processing of image control points | |
| CN116721226A (en) | Real estate unit dynamic modeling and ownership information association method based on live-action three-dimension | |
| WO2023123583A1 (en) | Fusion method, apparatus and system suitable for geological data and geographic information data | |
| CN106846388B (en) | Image map registration method and system | |
| CN114490907A (en) | Method and device for constructing famous city management database and storage medium | |
| CN112445951A (en) | Ocean data visualization method compatible with various maps and capable of being switched quickly and flexibly | |
| CN110019632A (en) | A kind of geographic element shape analysis method based on shape tree | |
| CN116302579B (en) | Space-time big data efficient loading rendering method and system for Web end | |
| CN117313217B (en) | Building indoor model modeling method based on three-dimensional scene data | |
| CN116206073A (en) | Data processing method and device, electronic equipment and storage medium | |
| CN115840927A (en) | Data fusion method and device of BIM model and GIS model and electronic equipment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |