Vector tile-based spatial data mapping system and method

Abstract

The application discloses a spatial data mapping system and method based on vector tiles, comprising a spatial data processing subsystem, a spatial information analysis subsystem and a spatial mapping configuration subsystem; the spatial data processing subsystem is used for preliminary extraction, storage of data and real-time update when data changes; the spatial information analysis subsystem is used for receiving data information in the spatial data processing subsystem, the spatial information analysis subsystem analyzes mixed overlap degree, buffer zone and density of the received data, and transmits the analyzed data; the spatial mapping configuration subsystem is used for extraction and matching of the spatial information analysis subsystem according to system requirements and display after completion; on the basis of vector tiles, the application sets different levels to utilize the combined display method to reduce the storage of space and effectively utilize the space.

Description

TECHNICAL FIELD

[0001] The application relates to the technical field of spatial data mapping, in particular to a spatial data mapping system and method based on vector tiles. BACKGROUND

[0002] With the continuous development of the national economy, continuous improvement and enrichment of national basic geographic information data, a map as a carrier of data information visualization plays an increasingly important role in the decision of the upper level. Spatial map drawing mainly includes raster tile and vector tile technologies. When the raster tile technology is used, the style of geographic elements is not easy to update, and a high-resolution map is not easy to obtain, so the vector tile technology is usually used in practice.

[0003] However, in the actual map drawing process, when the same land has different properties due to setting planning, the phenomenon of repeated mapping and waste of storage space occurs, the abstract ground objects cannot be analyzed from multiple angles and quantified on the display level of the spatial data image, the influence relationship caused by the interlaced relationship between different ground objects due to the environment and natural climate cannot be visualized, and the actual survival relationship between the characters and the environment cannot be directly and simply observed. SUMMARY

[0004] The application aims to provide a spatial data mapping system and method based on vector tiles to solve the problems in the background technology.

[0005] In order to solve the above technical problems, the application provides the following technical scheme: a spatial data mapping system based on vector tiles, comprising a spatial data processing subsystem, a spatial information analysis subsystem and a spatial mapping configuration subsystem.

[0006] The spatial data processing subsystem is used for preliminary extraction, storage of data and real-time update when the data changes; the spatial information analysis subsystem is used for receiving data information in the spatial data processing subsystem, the spatial information analysis subsystem analyzes the mixed overlap degree, buffer zone and density of the received data, and transmits the analyzed data; and the spatial mapping configuration subsystem is used for extracting and matching the spatial information analysis subsystem according to system requirements and displaying the completed data.

[0007] Further, the spatial data processing subsystem comprises a data integration module, a data management module and a data update module.

[0008] The data integration module is used for integrating spatial data, the integration including comparison of spatial data and integration of spatial data, the spatial data including regional blocks, feature names and population density, the feature names being rivers, lakes, buildings and vegetation; the data integration module compares and judges the regional blocks, the feature names and the population density in grades and classifies them into corresponding sets for integration, and transmits the integrated data to the data management module;

[0009] The data management module includes a central database, which receives the data integrated by the data integration module, and transmits the spatial data to the data updating module in real time;

[0010] The data updating module is used for updating spatial data, receives the spatial data from the central database to judge whether the data need to be updated, transmits the data needing to be updated to the data integration module, and transmits the data re-integrated by the data integration module to the data management module.

[0011] Further, the spatial information analysis subsystem includes a spatial overlay analysis module, a spatial buffer analysis module and a spatial density analysis module;

[0012] The spatial overlay analysis module is used for analyzing the mixed overlap degree of data information in the spatial data processing subsystem, and re-plans the spatial data having the mixed overlap degree; and the spatial overlay analysis module stores the re-planned spatial data into the spatial data processing subsystem;

[0013] The spatial buffer analysis module is used for buffer analysis of the spatial data re-planned by the spatial overlay analysis module, analyzes the spatial data one by one based on the data of the spatial overlay analysis module, and performs corresponding buffer processing; the spatial buffer analysis module analyzes the spatial density in the analysis, and transmits a signal to the spatial density analysis module;

[0014] The spatial density analysis module is used for specific analysis of the density of the spatial data based on the spatial buffer analysis module.

[0015] Further, the spatial mapping configuration subsystem includes a demand extraction module and a matching mapping module;

[0016] The demand extraction module is used for extracting mapping demands of the system, and the mapping demands are extracted from the data in the spatial information analysis subsystem; the demand extraction module transmits a demand signal to the matching mapping module; and the matching mapping module is used for matching the demands in the demand extraction module, thereby producing corresponding images.

[0017] The spatial data mapping method based on vector tiles includes the following steps:

[0018] Step S100: the collected spatial data is processed, the processing includes comparison and integration of the spatial data; the integrated data is stored; when new spatial data needs to be updated, the new spatial data needs to be compared and integrated again, and the new integrated data is stored again;

[0019] Step S200: based on the spatial data stored in step S100, analysis is performed, which includes overlay analysis of the spatial data, and whether the spatial data has mixed overlap is determined;

[0020] Step S300: based on the spatial data in step S200, the analyzed data is mapped and matched, which includes demand extraction and image production.

[0021] Further, the comparison and integration of the spatial data in step S100 includes the following specific process:

[0022] Step S110: the collected spatial data is compared and layered, the spatial data includes regional blocks, feature names and population density; the data packet capacity of the regional blocks is denoted as a set {Q i}, i is a natural number, Q i represents the i-th regional block in the set {Q i}; the data packet capacity of the feature names is denoted as a set {D i}, D i represents the feature name corresponding to the i-th regional block in the set {D i}; the data packet capacity of the population density is denoted as a set {R i}, R i represents the population density corresponding to the i-th regional block in the set {R i};

[0023] Layered calculation is beneficial to the convenience and neatness of spatial data analysis, and the analysis of each level of data after layering is more comprehensive, and the extraction and data storage after layering are more space-saving.

[0024] Step S120: based on the set in step S110, the three sets are calculated respectively, i.e., {Q i}∩{D i}=a1, {Q i}∩{R i}=a2, {D i}∩{R i}=a3, and a1, a2, a3∈{{Q i}、{D i}、{R i}};

[0025] In the set intersection process, the data packet capacity of the set corresponding to the region block, the feature name and the population density can be determined. The region block can contain the feature name spatial data and the population density spatial data. When the feature name and the population density are both concrete spatial data, there is a possibility of mutual inclusion. Therefore, the intersection of the sets is used for determination.

[0026] Step S130: Determine the number of {Q i}, {D i}, and {R i}. When the number of {{Q i}, {D i}, and {R i}} is in {2, 1, 0}, the level of the set {Q i}, the set {D i}, and the set {R i} can be determined. When a1={D i}, a2={R i}, and a3={R i}, the set {Q i} is the first level vector tile, the set {D i} is the second level vector tile, and the set {R i} is the third level vector tile.

[0027] The above sets are sequentially intersected. According to the rules of the sets, the intersection of the largest set with the other two sets is one of the two sets, and the largest set will not appear after sequential intersection. The number of intersection set possibilities is used to determine the size of the set for layering, which ensures that the sets have different sizes, excludes the case of the same set, and maximizes the set.

[0028] Further, step S200 includes the following specific processes:

[0029] Step S210: Extract the spatial data of the first level vector tile. The region blocks in the region block set divided according to the soil type are denoted as set , and the region blocks in the region block set divided according to the city function type are denoted as set .

[0030] The soil type and the city function are used to divide all region blocks to achieve the diversity of spatial data and avoid single direct data acquisition. Therefore, two ways of division are adopted for region block acquisition. However, the two divisions also bring spatial pressure to data storage and complexity to data query, which causes a certain degree of overlap in region block division.

[0031] Step S220: Calculate the set The number of isolated region blocks in the set i1 , the set The number of isolated region blocks in the set i2 , using the formula x i = j i1 + j i2 - 1, evaluate the relationship between x i and The number of new isolated region blocks y i actually produced when overlapping, x i represents the i The number of new isolated region blocks when overlapping with is only increased by one, y i represents The number of new isolated region blocks when overlapping with ;

[0032] Step S230: Based on the spatial data of the step S220, when y i ≤ x i , record and There is no mixed overlap, when y i > x i , record and There is mixed overlap; overlap and simplify the storage of the region blocks that do not produce mixed overlap, that is, merge the new boundaries of the isolated region blocks in the region blocks; combine and store the new region blocks that produce mixed overlap;

[0033] It is set that the boundary when only one isolated region block is increased is the boundary that produces mixed overlap, because when two isolated region blocks of the region blocks are overlapped, there are two cases, one is that the complete overlap still retains the original two region blocks, and the other is that one region block divides one of the isolated region blocks of the other region block, and only divides into two new isolated region blocks, so this case does not produce mixed overlap, and when storing, in order to save space, the divided boundary can be simplified to still be two isolated region blocks; when more than one isolated region block is increased, it means that the isolated region blocks have crossed and mixed, then the original isolated region blocks are combined and stored, which saves the storage space and facilitates the diversity of the mapping system data;

[0034] Step S240: Based on the mixed overlap of step S230, analyze the buffer area of the spatial data of the new region blocks with mixed overlap.

[0035] Buffer analysis is performed on new regions with mixed overlap because the overlapping attributes of the new regions after mixing increase, containing multiple spatial data. Buffer analysis of these spatial data can effectively refine the different attribute features on different regions, thus enabling better drawing of detailed spatial data images.

[0036] Furthermore, step S240 includes the following specific processes:

[0037] Step S241: Define a new set of region blocks with mixed overlap as N, and define isolated region blocks in set N as N0. i N i Let N represent the i-th isolated region block in set N. i The number of isolated blocks in the n j In the orphan block n j The second-level vector tiles are analyzed, which include rivers, lakes, buildings, and vegetation;

[0038] Step S242: Based on the spatial data from step S241, the rivers and lakes in the second-level vector tiles are planned as spatial lines, the buildings in the second-level vector tiles are planned as spatial points, and the vegetation in the second-level vector tiles is planned as spatial surfaces; define the spatial line as l. j , l j Indicates the nth j The j-th spatial line in an isolated block; define the spatial plane as s. j s j Indicates the nth j The j-th spatial plane in an isolated block; define spatial point d. j d j Indicates the nth j The j-th spatial point in an isolated block;

[0039] The abstract land feature names in the second-level vector tiles are visualized, and the analysis is effectively combined with the points, lines and planes in space. Furthermore, the abstract elements in the spatial image are quantitatively analyzed to obtain a specific and clear spatial image.

[0040] Step S243: Based on the spatial data from step S242, determine the spatial line l j With space surface s j Do they intersect? If they intersect, let spatial line l be... j With space surface s j Let point O be the point where the two intersecting planes lie on the edge of the space plane. Draw a tangent line to the space plane from point O. The space plane and this tangent line will form two angles, α1 and α2, where α1 + α2 = 180°, and α1 < α2. Using the formula... Calculate the actual angle threshold difference; compare this actual angle threshold difference with the system preset angle difference value. When comparing, When, the system sets up a green buffer zone on one side of the spatial line with angle α1; when If no greenery buffer zone is set, proceed to step S244 if the plants do not intersect.

[0041] The actual ground relationship between rivers, lakes and vegetation is transformed into an angle relationship between lines and surfaces. The size of the angle determines the amount of vegetation distributed on both sides of the river and lake. A buffer zone is set up on the side with less vegetation to prevent water level rise caused by natural rain disasters from causing the embankment to overflow. Setting up a buffer zone on the side with less vegetation can effectively observe areas where danger may occur from the image.

[0042] Step S244: Based on the spatial data from step S243, analyze spatial point d. j With space surface s j The relationship between occupancy rate.

[0043] Furthermore, step S244 includes the following specific processes:

[0044] Step S2441: Extract the horizontal distance a from the center of the spatial surface and the distance b from its corresponding perpendicular intersection. Calculate the difference between |ab| using the formula and compare it with the system's preset difference range. When |ab| is greater than or equal to the system's preset difference range, compensate the spatial surface into a rectangle, in which case a>b; when |ab| is less than the system's preset difference range, compensate the spatial surface into a circle, in which case a=b.

[0045] When the horizontal and vertical distances differ by a certain range, compensating with a rectangle is an effective way to reduce errors. When the difference between the two distances is small, planning it as a circle can reduce errors. Furthermore, planning the spatial surface as a regular shape is for the effective analysis of the spatial surface, which is beneficial for the calculation of population density and the extraction of spatial points.

[0046] Step S2442: Based on the spatial data in step S2441, the circular or rectangular spatial surface is uniformly divided into four equal spatial surface sub-regions U according to the horizontal and vertical intersection lines. e e={1,2,3,4}, and then divide the horizontal and vertical intersection lines into spatial surface subgrids u at equal intervals. i Extract the number of spatial points on all spatial face subgrids into a set Y = {u1, u2, u3, ... u}. i}, calculate the total number of spatial points in the four spatial sub-regions respectively: Using the formula Calculate the maximum number of spatial points corresponding to the four spatial sub-regions;

[0047] The maximum value of the spatial point number corresponding to the spatial surface sub-area is extracted to determine which area block has the highest building density, and then infer that the population density contained in the area block with the highest building density is also the highest, which is beneficial to the analysis of the population density in the following;

[0048] Step S2443: Based on the spatial data in step S2442, extract the spatial surface sub-area U containing the maximum spatial point number e , and extract the population number in the third level vector tile of the spatial surface sub-area as W e , calculate the area of the spatial surface sub-area by using the formula or , and obtain the per capita green plant occupancy rate by using the formula or .

[0049] The per capita green plant occupancy rate is analyzed and calculated to highlight the greening environment of the area in the spatial data image, and the per capita green plant occupancy rate also effectively expresses the rationality and improvement of the urban planning, so that the display of the spatial data image not only accurately describes the spatial geographical position, but also effectively expresses the improvement of the humanistic feelings and people's livelihood environment.

[0050] Step S2444: Based on the spatial data in step S2443, compare the green plant occupancy rate of each person with the system preset green plant standard occupancy rate; when the green plant occupancy rate of each person is greater than or equal to the system preset green plant standard occupancy rate, the system does not display the buffer area; when the green plant occupancy rate of each person is less than the system preset green plant standard occupancy rate, the system displays the buffer area prompt in the corresponding spatial surface sub-area.

[0051] Further, step S300 includes the following specific processes:

[0052] Step S310: The demand extraction includes the display of the first level vector tile, the display of the second level vector tile and the display of the third level vector tile;

[0053] Step S320: Based on the demand in step S310, when the spatial image is made for the first level vector tile, the spatial image will display the different attributes before the level in detail, and different colors are displayed for different attributes; when the spatial image is made for the second level vector tile, the spatial image displays the corresponding area block containing the buffer area in a virtual display; when the spatial image is made for the third level vector tile, the spatial image displays the per capita green plant occupancy rate on the corresponding area block; when the demand is combined and extracted, the corresponding mapping can also be combined.

[0054] Compared with the prior art, the present application has the advantages that: on the basis of vector tiles, different levels are respectively set to use the combined display method to reduce the storage of space, and effectively use the space; secondly, the spatial data is quantized at multiple angles, abstract features are converted into concrete features, the mathematical relationship between features is calculated, and the result is visually displayed on the image; finally, the population density is combined with the environmental feature relationship to effectively display on the spatial data image, and the survival relationship and effective utilization relationship between man and environment are displayed. BRIEF DESCRIPTION OF DRAWINGS

[0055] The accompanying drawings are included to provide a further understanding of the present application, and constitute a part of the specification, illustrate the present application and explain the technical scheme of the present application, and do not constitute a limitation on the present application. In the drawings:

[0056] Figure 1 is a structural schematic diagram of the spatial data mapping system based on vector tiles of the present application;

[0057] Figure 2 is a mapping overall flowchart of the spatial data mapping method based on vector tiles of the present application;

[0058] Figure 3 is a comparison and integration step diagram of spatial data of the spatial data mapping method based on vector tiles of the present application;

[0059] Figure 4 is a superimposed analysis step diagram of spatial data of the spatial data mapping method based on vector tiles of the present application;

[0060] Figure 5 is a buffer analysis step diagram of spatial data of the spatial data mapping method based on vector tiles of the present application;

[0061] Figure 6 is an occupancy rate analysis step diagram of spatial data of the spatial data mapping method based on vector tiles of the present application;

[0062] Figure 7 is a mapping matching step diagram of the spatial data mapping method based on vector tiles of the present application. DETAILED DESCRIPTION

[0063] The technical scheme in the embodiments of the present application will be described clearly and completely below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor fall within the scope of protection of the present application.

[0064] Please refer to Figures 1-7The application provides a technical scheme: a spatial data mapping system based on vector tiles, comprising a spatial data processing subsystem, a spatial information analysis subsystem and a spatial mapping configuration subsystem;

[0065] The spatial data processing subsystem is used for preliminary extraction, storage and real-time update of data when data changes; the spatial information analysis subsystem is used for receiving data information in the spatial data processing subsystem, the spatial information analysis subsystem analyzes the received data in terms of mixed overlap, buffer zone and density, and transmits the analyzed data; and the spatial mapping configuration subsystem is used for extraction and matching of the spatial information analysis subsystem according to system requirements and then display.

[0066] The spatial data processing subsystem comprises a data integration module, a data management module and a data update module.

[0067] The data integration module is used for integrating spatial data, which includes comparison and integration of spatial data, and the spatial data contains regional blocks, feature names and population density, the feature names are rivers, lakes, buildings and vegetation; the data integration module compares and judges the regional blocks, feature names and population density in grades and classifies them into corresponding sets for integration, and transmits the integrated data to the data management module.

[0068] The data management module comprises a central database, which receives the integrated data of the data integration module and transmits the spatial data to the data update module in real time.

[0069] The data update module is used for updating spatial data, receives the spatial data of the central database to judge whether the data needs to be updated, transmits the data needing to be updated to the data integration module, and transmits the re-integrated data to the data management module again.

[0070] The spatial information analysis subsystem comprises a spatial overlay analysis module, a spatial buffer zone analysis module and a spatial density analysis module.

[0071] The spatial overlay analysis module is used for analyzing the mixed overlap of data information in the spatial data processing subsystem, and the spatial overlay analysis module re-plans the spatial data with mixed overlap; the spatial overlay analysis module stores the re-planned spatial data to the spatial data processing subsystem.

[0072] The spatial buffer analysis module is used for buffer analysis on the spatial data re-planned by the spatial overlay analysis module, the spatial buffer analysis module analyzes the spatial data one by one based on the data of the spatial overlay analysis module and performs corresponding buffer processing; the spatial buffer analysis module performs spatial density analysis in the analysis, and the spatial buffer analysis module transmits a signal to the spatial density analysis module;

[0073] The spatial density analysis module is used for specific density analysis on the spatial data based on the spatial buffer analysis module.

[0074] The spatial mapping configuration subsystem includes a demand extraction module and a matching mapping module;

[0075] The demand extraction module is used for extracting the mapping demand of the system, the mapping demand is used for extracting the data in the spatial information analysis subsystem, and the demand extraction module transmits a demand signal to the matching mapping module; the matching mapping module is used for corresponding matching of the demand in the demand extraction module, so as to draw a corresponding image.

[0076] The spatial data mapping method based on vector tiles includes the following steps:

[0077] Step S100: The collected spatial data is processed, the related processing includes comparison and integration of the spatial data; the integrated data is stored; after the storage is completed, when new spatial data needs to be updated, the new spatial data needs to be compared and integrated again, and the new integrated data is stored again;

[0078] The comparison and integration of the spatial data in step S100 include the following specific processes:

[0079] Step S110: The collected spatial data is compared and layered, the spatial data includes regional blocks, feature names and population density; the data packet capacity of the regional blocks is denoted as a set {Q i}, i is a natural number, Q i represents the i th regional block in the set {Q i}; the data packet capacity of the feature name is denoted as a set {D i}, D i represents the feature name corresponding to the i th regional block in the set {D i}; the data packet capacity of the population density is denoted as a set {R i}, R i represents the population density corresponding to the i th regional block in the set {R i};

[0080] The hierarchical calculation is beneficial to the convenience and neatness of the spatial data analysis. After the spatial data is divided into layers, the analysis of each layer data is more comprehensive, and the extraction and data storage after the division are more space-saving.

[0081] Step S120: Based on the set in step S110, the set intersection of the three is calculated respectively, that is, i}∩{D i}=a1, {Q i}∩{R i}=a2, {D i}∩{R i}=a3, and a1, a2, a3∈{{Q i}、{D i}、{R i}}.

[0082] In the set intersection process, the data packet capacity of the sets corresponding to the region block, the feature name and the population density can be judged. Since the region block can contain the feature name spatial data and the population density spatial data, when the feature name and the population density both belong to the concretized spatial data, there is a possibility of mutual inclusion, so the intersection of the set is used for judgment.

[0083] Step S130: The number of {Q i}、{D i}、{R i} is judged. When the number of {{Q i}、{D i}、{R i}} is ∈{2, 1, 0}, it is determined that the level of the set {Q i}, the set {D i} and the set {R i} is set. When a1={D i}, a2={R i} and a3={R i}, the set {Q i} is the first level vector tile, the set {D i} is the second level vector tile, and the set {R i} is the third level vector tile.

[0084] For example: In the intersection process, when {Q i}∩{D i}={Q i}, {Q i}∩{R i}={R i}, {D i}∩{R i}={Ri}, this time meet the existence of { {Q i}, {D i}, {R i}} The number of ∈ {2, 1, 0}, but not the system settings corresponding to; When {Q i} ∩ {D i}=0, {Q i} ∩ {R i}=0, {D i} ∩ {R i}={R i}, not meet the existence of { {Q i}, {D i}, {R i}} The number of ∈ {2, 1, 0} ;

[0085] The above set is intersected in turn, and according to the rule of the set, it can be obtained that the intersection of the largest set and the other two sets is one of the two sets, and the largest set will not appear after the intersection in turn; The number of intersection set possibilities is used to judge the size of the set to divide it into layers, which ensures that the set has size difference, excludes the case of the same set, and makes the set largest.

[0086] Step S200: based on the spatial data stored in step S100, analysis includes overlay analysis of spatial data, and judges whether the spatial data exists mixed overlap degree;

[0087] Step S200 includes the following specific process:

[0088] Step S210: extracting the spatial data of the first level vector tile, the region block set divided according to the soil type in the region block set is denoted as set , the region block set divided according to the city function type in the region block set is denoted as set ;

[0089] Dividing all region blocks according to soil type and city function is to realize the diversity of spatial data, avoid the single direct way of data acquisition, so two ways of division are adopted for region block acquisition, but at the same time, the two divisions bring space pressure to data storage, and complexity to data query, so that the division of region blocks appears a certain degree of overlap;

[0090] Step S220: calculating the number of isolated region blocks in set is j i1 , the number of isolated region blocks in set is j i2 , using the formula x i =j i1 +ji2 -1, evaluate x i with the number of new isolated region blocks actually generated when the overlap of i x i represents the i-th with the number of new isolated region blocks when the overlap of i represents with the number of new isolated region blocks actually generated when the overlap of

[0091] Step S230: based on the spatial data of step S220, when y i ≤x i , record with no mixed overlap degree, when y i >x i , record with mixed overlap degree; the region blocks that do not produce mixed overlap degree are overlapped and simplified to store, that is, the new boundaries of isolated region blocks in the region blocks are merged; the new region blocks that produce mixed overlap degree are combined and stored;

[0092] For example , , when y i =4, the new region block is= ;

[0093] The boundary is set when only one isolated region block is added to produce mixed overlap degree, because when two isolated region blocks are overlapped, there are two cases, one is complete overlap and still retains the original two region blocks, the other is that one region block divides one of the isolated region blocks of the other region block, and only divides it into two new isolated region blocks, so this case does not produce mixed overlap, and when storing, in order to save space, the boundary of the division can be simplified to still be two isolated region blocks; when more than one isolated region block is added, it means that the isolated region blocks have crossed and mixed, then the original isolated region blocks are combined and stored, which saves the storage space and facilitates the diversity of the mapping system data;

[0094] Step S240: based on the mixed overlap degree of step S230, the new region blocks with mixed overlap degree are analyzed for the buffer of the spatial data.

[0095] Buffer analysis is performed on new regions with mixed overlap because the overlapping attributes of the new regions after mixing increase, containing multiple spatial data. Buffer analysis of these spatial data can effectively refine the different attribute features on different regions, thus enabling better drawing of detailed spatial data images.

[0096] Step S240 includes the following specific processes:

[0097] Step S241: Define a new set of region blocks with mixed overlap as N, and define isolated region blocks in set N as N0. i N i Let N represent the i-th isolated region block in set N. i The number of isolated blocks in the n j In the orphan block n j The second-level vector tiles are analyzed, which include rivers, lakes, buildings, and vegetation;

[0098] Step S242: Based on the spatial data from step S241, the rivers and lakes in the second-level vector tiles are planned as spatial lines, the buildings in the second-level vector tiles are planned as spatial points, and the vegetation in the second-level vector tiles is planned as spatial surfaces; define the spatial line as l. j , l j Indicates the nth j The j-th spatial line in an isolated block; define the spatial plane as s. j s j Indicates the nth j The j-th spatial plane in an isolated block; define spatial point d. j d j Indicates the nth j The j-th spatial point in an isolated block;

[0099] The abstract land feature names in the second-level vector tiles are visualized, and the analysis is effectively combined with the points, lines and planes in space. Furthermore, the abstract elements in the spatial image are quantitatively analyzed to obtain a specific and clear spatial image.

[0100] Step S243: Based on the spatial data from step S242, determine the spatial line l j With space surface s j Do they intersect? If they intersect, let spatial line l be... j With space surface s j Let point O be the point where the two intersecting planes lie on the edge of the space plane. Draw a tangent line to the space plane from point O. The space plane and this tangent line will form two angles, α1 and α2, where α1 + α2 = 180°, and α1 < α2. Using the formula... Calculate the actual angle threshold difference; compare this actual angle threshold difference with the system preset angle difference value. When , the system sets a green buffer zone on one side of the spatial line with an angle of a1; when , no green buffer zone is set; if not intersecting, proceed to step S244;

[0101] For example: the spatial line l j intersects with the spatial surface s j , , the system presets an angle difference , at this time , so , the system does not need to set a buffer zone;

[0102] The actual ground relationship between rivers, lakes and vegetation is converted into the angle relationship between lines and surfaces, the amount of vegetation distributed on both sides of rivers and lakes is determined by the size of the angle, and a buffer zone is set on the side with less vegetation distribution to prevent water level rise caused by natural rain disasters, resulting in the occurrence of dike overflow. Setting a buffer zone on the side with less vegetation distribution can effectively observe the possible dangerous areas from the image;

[0103] Step S244: based on the spatial data in step S243, analyze the occupancy relationship between the spatial point d j and the spatial surface s j .

[0104] Step S244 includes the following specific processes:

[0105] Step S2441: extract the horizontal distance of the spatial surface center position as a, and the distance of the vertical intersection line corresponding thereto as b, and compare the size of the difference value calculated by the formula with the system preset difference value range; when |a-b| is greater than or equal to the system preset difference value range, the spatial surface is compensated as a rectangle, at this time a>b; when |a-b| is less than the system preset difference value range, the spatial surface is compensated as a circle, at this time a=b;

[0106] When the horizontal distance and the vertical distance differ within a certain range, compensating into a rectangle is an effective method to reduce errors, and planning into a circle when the distance difference between the two is small can reduce errors. Planning the spatial surface into a regular shape is for effective analysis of the spatial surface, which is beneficial to population density calculation and spatial point extraction;

[0107] Step S2442: based on the spatial data in step S2441, the circular spatial surface or rectangular spatial surface is uniformly divided into four equal parts of the spatial surface sub-area U e , e={1, 2, 3, 4}, and the horizontal intersection line and the vertical intersection line are respectively equally divided to make the spatial surface sub-grid u i; extract the number of spatial points on all spatial face subgrids as a set Y = {u1, u2, u3,... u i}, and calculate the total number of spatial points in the four spatial face subregions as follows: , and calculate the maximum value of the number of spatial points corresponding to the four spatial face subregions using the formula ;

[0108] For example, the horizontal intersection line and the vertical intersection line are equally divided into two parts, each spatial face subregion contains four spatial face subgrids, and Y = {5, 2, 6, 1, 2, 5, 3, 7, 1, 2, 5, 4, 3, 2, 6, 1}, then t1 = ∑u i = 5 + 2 + 6 + 1 = 14, t2 = ∑u i = 2 + 5 + 3 + 7 = 17, t3 = ∑u i = 1 + 2 + 5 + 4 = 12, and t4 = ∑u i = 3 + 2 + 6 + 1 = 12; Maxu = {t e = ∑u i} = 17;

[0109] The maximum value of the number of spatial points corresponding to the spatial face subregion is extracted in order to determine which regional block has the highest building density first, and then infer that the population density contained in the regional block with the highest building density is also the highest, which is beneficial to the analysis of population density in the following text;

[0110] Step S2443: Based on the spatial data in step S2442, extract the spatial face subregion U e with the maximum number of spatial points, and extract the number of population in the third level vector tile of the spatial face subregion as W e ; calculate the area of the spatial face subregion using the formula or ; and obtain the per capita green plant occupancy rate using the formula or ;

[0111] For example, when the horizontal distance is a = 100 m, the vertical distance b = 105 m, the system preset difference range is 15 m, and |a - b| = |100 - 105| = 5 < 15 at this time, the spatial face is circular, and the radius a = 100 m, then ; the number of population W e = 200, then ;

[0112] The purpose of analyzing and calculating the per capita green space occupancy rate is to highlight the green environment of a region in spatial data images. The per capita green space occupancy rate also effectively expresses the rationality and improvement of the city's planning. This makes the display of spatial data images not only a precise description of spatial location, but also an effective expression of humanistic sentiments and the direction of improvement of people's livelihood and environment.

[0113] Step S2444: Based on the spatial data in step S2443, compare each person's green plant occupancy rate with the system's preset green plant standard occupancy rate; when each person's green plant occupancy rate is greater than or equal to the system's preset green plant standard occupancy rate, the system does not display the buffer; when each person's green plant occupancy rate is less than the system's preset green plant standard occupancy rate, the system displays a buffer reminder in the corresponding spatial sub-area.

[0114] Step S300: Based on the spatial data in step S200, perform map matching on the analyzed data. Map matching includes extracting requirements and creating images.

[0115] Step S300 includes the following specific processes:

[0116] Step S310: The extraction of requirements includes the display of vector tiles at the first level, the display of vector tiles at the second level, and the display of vector tiles at the third level;

[0117] Step S320: Based on the requirements of step S310, when creating a spatial image for the first-level vector tiles, the spatial image will display the different attributes that constitute the first level in detail, and display different colors for different attributes; when creating a spatial image for the second-level vector tiles, the spatial image will display the corresponding area blocks containing the buffer in a blurred display; when creating a spatial image for the third-level vector tiles, the spatial image will display the per capita green plant occupancy rate on the corresponding area blocks; when the requirement is to combine and extract, the corresponding maps can also be combined.

[0118] It should be noted that, in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such process, method, article, or apparatus.

[0119] Finally, it should be noted that the above only describes the preferred embodiments of the present application and is not intended to limit the present application. Although the present application has been described in detail with reference to the foregoing embodiments, those skilled in the art will appreciate that the technical solutions described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalent ones. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A vector tile-based spatial data mapping system, characterized by, The system comprises a spatial data processing subsystem, a spatial information analysis subsystem and a spatial mapping configuration subsystem; The spatial data processing subsystem is used for preliminary extraction, storage and real-time update of data when data changes; The spatial information analysis subsystem is used for receiving data information in the spatial data processing subsystem, and performing mixed overlap analysis, buffer zone analysis and density analysis on the received data, and transmitting the analyzed data, and the spatial mapping configuration subsystem is used for extracting and matching the spatial information analysis subsystem according to system requirements and displaying the matched data; The spatial mapping configuration subsystem comprises a demand extraction module and a matching mapping module; The demand extraction module is used for extracting mapping requirements of the system, the mapping requirements are used for extracting data in the spatial information analysis subsystem, and the demand extraction module transmits a demand signal to the matching mapping module; and the matching mapping module is used for matching the demand in the demand extraction module, so as to generate a corresponding image.

2. The vector tile-based spatial data mapping system of claim 1, wherein: The spatial data processing subsystem comprises a data integration module, a data management module and a data update module; The data integration module is used for integrating spatial data, the integration comprises comparison and integration of spatial data, and the spatial data comprises regional blocks, feature names and population density, the feature names are rivers, lakes, buildings and vegetation; the data integration module compares and judges the levels of the regional blocks, feature names and population density, and classifies them into corresponding sets for integration, and the data integration module transmits the integrated data to the data management module; The data management module comprises a central database, the central database receives the integrated data of the data integration module, and the central database transmits spatial data to the data update module in real time; The data update module is used for updating spatial data, the data update module receives spatial data of the central database to judge whether the spatial data need to be updated, the data update module transmits the data needing to be updated to the data integration module, and the data integration module transmits the re-integrated data to the data management module again.

3. The vector-tile-based spatial data mapping system of claim 1, wherein: The spatial information analysis subsystem comprises a spatial overlay analysis module, a spatial buffer zone analysis module and a spatial density analysis module; The spatial overlay analysis module is used for analyzing the mixed overlap degree of data information in the spatial data processing subsystem, and the spatial overlay analysis module re-plans spatial data with the mixed overlap degree; The spatial overlay analysis module stores the re-planned spatial data into the spatial data processing subsystem; The space buffer analysis module is used for buffer analysis on the space data re-planned by the space overlay analysis module, and the space buffer analysis module sequentially analyzes the individual space data based on the data of the space overlay analysis module and performs corresponding buffer processing; the space buffer analysis module performs spatial density analysis in the analysis, and the space buffer analysis module transmits a signal to the space density analysis module; The space density analysis module is used for specific density analysis on the space data based on the space buffer analysis module.

4. A method for cartography of spatial data based on vector tiles, characterized in that, The method comprises the following steps: Step S100: performing relevant processing on the collected space data, the relevant processing comprising comparison and integration of the space data; storing the integrated data; when new space data needs to be updated after the storage is completed, the new space data needs to be compared and integrated again, and the new integrated data is stored again; Step S200: performing analysis based on the stored space data in step S100, the analysis comprising overlay analysis on the space data, and determining whether the space data has mixed overlap; Step S300: performing mapping matching on the analyzed data based on the space data in step S200, the mapping matching comprising demand extraction and image production.

5. The vector-tile-based spatial data mapping method according to claim 4, characterized in that: The comparison and integration of the space data in step S100 comprises the following specific process: Step S110: comparing and layering the collected spatial data, including regional blocks, feature names and population density; let the data packet capacity of regional blocks be a set {Q i}, i is a natural number, Q i represents the i-th regional block in the set {Q i}; let the data packet capacity of feature names be a set {D i}, D i represents the feature name corresponding to the i-th regional block in the set {D i}; let the data packet capacity of population density be a set {R i}, R i represents the population density corresponding to the i-th regional block in the set {R i}; Step S120: based on the set in the step S110, respectively, three set intersection calculation, namely {Q i}∩{D i}=a1, {Q i}∩{R i}=a2, {D i}∩{R i}=a3, and a1, a2, a3∈{{Q i}、{D i}、{R i}}; Step S130: judging the number of {Q i} i} i} i} i} i}∈{2, 1, 0}, the level of the set {Q i}, the set {D i}, and the set {R i} can be determined; the system is set when a1={D i}, a2={R i}, and a3={R i}, the set {Q i} is the first level of vector tiles, the set {D i} is the second level of vector tiles, and the set {R i} is the third level of vector tiles.

6. The vector-tile-based spatial data mapping method according to claim 5, characterized in that: Step S200 comprises the following specific process: Step S210: extracting the spatial data of the first-level vector tiles, and denoting the region blocks divided according to the soil types in the region block set as set denoting the region blocks divided according to the city function types in the region block set as set ; Step S220: calculating the set The number of isolated region blocks in the set is j i1 , the set The number of isolated region blocks in the set is j i2 , the relationship between x i = j i1 + j i2 - 1 is evaluated, x i and The number of new isolated region blocks y i actually generated when the overlap of i The number of new isolated region blocks when the overlap of the i and only increases one isolated region block, y i The number of new isolated region blocks actually generated when the overlap of and is overlapped; Step S230: Based on the spatial data from step S220, when y i ≤x i At that time, I recorded and There is no mixed overlap when y i >x i At that time, I recorded and There is mixed overlap; regions that do not produce mixed overlap are simplified and stored by merging the new boundaries of isolated regions; new regions that produce mixed overlap are stored in combination. Step S240: performing buffer analysis on the space data of the new area block with mixed overlap based on the mixed overlap in step S230.

7. The vector-tile-based spatial data mapping method according to claim 6, characterized in that: Step S240 comprises the following specific process: Step S241: define a new set of region blocks with mixed overlap degree as N, the isolated region blocks in set N are N i , N i represents the i-th isolated region block in set N, and the isolated region blocks in isolated region block N i are n j , the second level vector tiles are analyzed in isolated region block n j , and the second level vector tiles include rivers and lakes, buildings, and vegetation; Step S242: Based on the spatial data from step S241, the rivers and lakes in the second-level vector tiles are planned as spatial lines, the buildings in the second-level vector tiles are planned as spatial points, and the vegetation in the second-level vector tiles is planned as spatial surfaces; define the spatial line as l. j , l j Indicates the nth j The j-th spatial line in an isolated block; define the spatial plane as s. j s j Indicates the nth j The j-th spatial plane in an isolated block; define spatial point d. j d j Indicates the nth j The j-th spatial point in an isolated block; Step S243: judging whether the spatial line l j intersects with the spatial face s j , if yes, let the point where the spatial line l j intersects with the spatial face s j and is located on the edge of the spatial face be o, and draw a tangent of the spatial face through the point o, at this time, the spatial line l forms two angles ɑ1 and ɑ2 with the tangent, and ɑ1+ɑ2=180º, and ɑ1<ɑ2 are defined; the real angle threshold difference is calculated by the formula The real angle threshold difference is compared with a system preset angle difference value When , the system sets a green buffer zone on one side of the space line with an angle of a1; when , no green buffer zone is set; if not intersecting, proceed to step S244; Step S244: Based on the spatial data of step S243, analyze the occupancy rate of the spatial point d j in relation to the spatial surface s j .

8. The vector-tile-based spatial data mapping method according to claim 7, characterized in that: Step S244 comprises the following specific process: Step S2441: extracting the horizontal distance a of the space surface center position and the distance b of the vertical intersection line corresponding to the space surface center position, and comparing the size of the difference value |a-b| calculated by the formula with the system preset difference value range; when |a-b| is greater than or equal to the system preset difference value range, the space surface is compensated as a rectangle, and at this time a>b; when |a-b| is less than the system preset difference value range, the space surface is compensated as a circle, and at this time a=b; Step S2442: Based on the spatial data in step S2441, the circular or rectangular space surface is evenly divided into four equal parts along the horizontal and vertical intersecting lines to form four space surface sub-regions U e , e={1, 2, 3, 4}, and the space surface sub-grid u is obtained by equally dividing the horizontal and vertical intersecting lines i ; the number of space points on all space surface sub-grids is extracted as a set Y={u1, u2, u3,......u i}, and the total number of space points in the four space surface sub-regions is calculated as: , and the maximum value of the number of space points corresponding to the four space surface sub-regions is calculated using the formula ; Step S2443: Based on the spatial data in the step S2442, extract the spatial surface sub-area U containing the maximum number of spatial points e , and extract the population number W in the vector tile of the third level of the spatial surface sub-area e , calculate the area of the spatial surface sub-area by the formula or ​ Using the formula or to obtain the greenery occupancy per capita; Step S2444: comparing the green plant occupancy rate of each person with the system preset green plant standard occupancy rate based on the space data in step S2443; when the green plant occupancy rate of each person is greater than or equal to the system preset green plant standard occupancy rate, the system does not display the buffer zone; when the green plant occupancy rate of each person is less than the system preset green plant standard occupancy rate, the system displays the buffer zone reminder in the corresponding space surface sub-area.

9. The vector-tile-based spatial data mapping method according to claim 8, characterized in that: Step S300 comprises the following specific process: Step S310: the demand extraction comprises display of the first level of vector tiles, display of the second level of vector tiles and display of the third level of vector tiles; Step S310: the demand extraction comprises display of the first level of vector tiles, display of the second level of vector tiles and display of the third level of vector tiles; Step S320: Based on the requirement of step S310, when making the spatial image of the first level vector tile, the spatial image will show the different attributes before the level in detail, and different colors are displayed for different attributes; when making the spatial image of the second level vector tile, the spatial image shows the corresponding area block of the buffer zone in virtual display; when making the spatial image of the third level vector tile, the spatial image shows the per capita green plant occupancy on the corresponding area block; when the requirement is combined for extraction, the corresponding mapping can also be combined.

Images