CN113689573B - DGGS interoperation quality evaluation method for geographic whole elements - Google Patents

Abstract

The invention discloses a DGGS interoperability quality evaluation method oriented to geographic full elements, which classifies and processes grid interoperability quality evaluation, analyzes the common features of geographic full elements expressed by grids for heterogeneous grids, and treats them as For area elements, study geometric similarity metrics, and design similarity vectors for geographic full elements. First, the geographical elements expressed in different grid systems are superimposed, and the interoperability quality is roughly evaluated by using the area ratio of the superposition area; the shape similarity is refined by establishing the shape descriptor between the regional boundary grids; finally, the Similarity measure vector to evaluate the quality of geographic full features in DGGS interoperability. For isomorphic grids, the geometric properties and dual relations of the grid itself are used for measurement. The area ratio and shape similarity of the overlapping area of the elements are selected to form the interoperability similarity measurement vector. Through the linear operation of the vector, the accuracy of the initial judgment of the grid interoperability quality is provided, which effectively improves the quality evaluation speed.

Description

DGGS Interoperability Quality Evaluation Method Oriented to All Geographical Elements

technical field

The invention relates to the technical field of global discrete grid interoperability quality evaluation, in particular to a DGGS interoperability quality evaluation method oriented to geographic full elements.

Background technique

Interoperability is the basis for data exchange and integration of different grid systems, and the basis for building a big data sharing platform. OGC's "Discrete Global Grid System Core Standard" (Discrete Global Grid System Core Standard) pointed out that interoperability and quality evaluation are important guarantees for expanding global grid applications, but the standard does not give specific evaluation principles and implementation methods . Therefore, how to perform quality evaluation has become one of the key issues that need to be solved urgently in the interoperability research of different grids.

At present, the research on the quality evaluation of geographic elements mainly focuses on the design and calculation of similarity metrics for vector data.

Point element quality evaluation. Depending on the similarity between points, it is generally represented by the distance between points. Common distance measurement methods mainly include Euclidean distance, Mahalanobis distance, Hausdorff distance, Manhattan distance, etc. The deviation between the distance and the threshold is the similarity.

Quality assessment of line features. For length, direction and curvature. The length is the sum of the distances between adjacent points in the line element; the direction is generally defined as the angle between the starting point of the line element and the horizontal direction; the curvature describes the curvature of the line element. The rate of change of length, direction, and curvature is used as the similarity. In addition, some researchers use the Fourier shape descriptor of the line element as its shape vector, and use the Euclidean distance between vectors as the similarity; the longest common subsequence measures the similarity of the line element.

The quality evaluation of area elements is mainly aimed at area and shape, where area is the basic feature of area elements, and commonly used area calculation indicators include MBR area, convex hull area, etc. A surface feature is composed of a closed line feature and its internal area, so the similarity measurement of the shape of the surface feature is the similarity measurement of the closed line feature. In addition to the above-mentioned common methods, it also includes using the boundary point to the center point The distance to construct the shape feature vector; extract the inner chord length, outer chord length and the average projection length of the arc to the chord to construct the chord feature matrix, and establish the shape similarity measurement model.

At present, there is no research on the quality evaluation of geographic elements in grid interoperability. In addition, the existing evaluation methods are independent of each other, evaluating point, line, and surface elements separately, lacking an overall comprehensive evaluation for all geographic elements. Geographical elements represented by grids have common characteristics. If the vector similarity measure is directly used for interoperability quality evaluation, the grid data needs to be vectorized, which loses the grid features of the grid and increases the computational complexity. At the same time, quality evaluation indicators are established for point, line, and surface elements, and there are repeated definitions and calculations in grid interoperability.

At present, there are few studies on the quality evaluation of geographic elements in grid interoperability, and the general idea is to directly use vector similarity measures for interoperability quality evaluation. This process requires grid data to be vectorized, and the grid characteristics of the grid are lost. , increasing the computational complexity.

Contents of the invention

Aiming at the defects of the prior art, the present invention provides a DGGS interoperability quality evaluation method oriented to all geographic elements.

In order to realize above object of the invention, the technical scheme that the present invention takes is as follows:

A quality evaluation method of DGGS interoperability oriented to geographic full elements, including: according to the difference of interoperable grid systems, grid interoperability is divided into heterogeneous grid interoperability and homogeneous grid interoperability. Among them, heterogeneous grid interoperability refers to the interoperability between grids with different structures, such as GeoSOT grid and spherical icosahedral triangular grid (Ico-TRI); isomorphic grid interoperability refers to grids with the same initial Interoperability between structured grids, such as Icosahedral hexagonal grid (Ico-HEX). The isomorphic grids have the same initial structure, so there is a simple geometric correspondence between the isomorphic grids, such as the dual relationship between Ico-TRI and Ico-HEX.

An interoperable quality evaluation method for heterogeneous grids, including: expressing geographic elements in different grid systems with known resolutions. The feature expression in the GeoSOT grid corresponds the geographic coordinates of the point, line and area features to the GeoSOT grid. The elements in the Ico-TRI grid are expressed, the spherical coordinates of the points are known, the point elements are projected onto the icosahedron datum using the Snyder projection, and the grid code is determined according to the grid resolution to realize the grid expression of the point elements. For line features, the coordinates of the points that make up the line feature are known, starting from the starting grid, searching for adjacent grids, finding the intersecting grid, using the currently determined grid as the starting grid, continuing to search until the end point, and realizing the line feature grid expression. For surface features, first determine the grid in the MRB of the surface feature, and then further determine the grid in the surface feature and the grid intersecting with the boundary of the surface feature to realize the grid expression of the surface feature.

For heterogeneous grids, the similarity measure vector is calculated according to the gridded results of the elements, and the gridded results of the elements are defined as the grid G _geosot and the grid G _tri , where G _geosot represents the grid set that expresses the geographical features by GeoSOT, G _tri represents the grid collection of geographical features represented by Ico-TRI;

For isomorphic grids, the geographical elements are first expressed in different grid systems with known resolutions. For the expression of elements in the Ico-HEX grid, refer to the Ico-TRI grid element expression method. Calculate the similarity measure vector according to the result of element gridding, and define the result of element gridding as grid G _tri and grid G _hex respectively, where G _hex represents the set of grids expressing geographical elements by Ico-HEX;

Further, the specific steps of calculating the similarity measure vector in the heterogeneous grid are as follows:

和/>

其中S表示面积，下标分别对应G_geosot和G_tri的相交面积、G_geosot格网面积和G_tri格网面积；S1: Calculation of the area ratio of the superposition area; according to the theory of spherical geometry, calculate the intersection area of the GeoSOT grid G _geosot and the Ico-TRI grid G _tri

and />

Where S represents the area, and the subscripts correspond to the intersection area of G _geosot and G _tri , the grid area of G _geosot and the grid area of G _tri , respectively;

S2: shape similarity;

Before the calculation, determine the boundary grids in the grids G _geosot and G _tri respectively, and obtain the boundary grid sets B _geosot and B _tri , where the point elements have no boundary grids, and the line elements are all boundary grids. For the boundary grid The shape similarity is calculated by the set, where B represents the boundary grid set, and the subscripts correspond to the GeoSOT and Ico-TRI grid systems respectively;

S3: construction of similarity measure vector;

The area ratio of the superimposed area and the shape similarity vector constitute the similarity measure vector (ε,σ _{Ggeosot→Gtri} , σ _{Gtri→Ggeosot} ), where ε represents the area ratio of the superimposed area, σ represents the distance variance, and the subscripts respectively represent GeoSOT→Ico-TRI and conversion of Ico-TRI→GeoSOT;

Further, the specific steps for calculating the similarity measure vector in the isomorphic grid are as follows:

S1: Define the dual relationship between triangular grid and hexagonal grid;

The duality relationship between triangle and hexagon is as follows: the apex of the triangle corresponds to the center of the hexagon, and the center of the triangle corresponds to the apex of the hexagon, and the three corners of the triangle are marked as one apex and two base angles. The point at the top corner corresponds to the center of the hexagon.

S2: similarity of geographic element area;

According to the dual relationship between triangular and hexagonal grids, the overlapping area of point elements can be directly obtained according to the dual relationship, the overlapping area of line elements depends on the direction of the G _tri grid, and the overlapping area of area elements depends on the boundary grid.

S3: geographic feature shape similarity;

Before calculation, determine the boundary grids in the grids G _tri and G _hex respectively, and obtain the boundary grid sets B _tri and B _hex , in which the point elements have no boundary grids, all the line elements are boundary grids, and some area elements exist Boundary grid, which calculates the shape similarity for the boundary grid set, where B _hex represents the boundary grid in the Ico-HEX expression geographic feature grid set G _hex ;

S4: construction of similarity measure vector;

The area ratio of the superimposed area and the shape similarity vector form a similarity measure vector (ε,σ _Gtri→Ghex ,σ _Gtri→Ghex ), where σ is the distance variance, and the subscripts represent the distance from Ico-TRI to Ico-HEX and Ico-HEX Conversion to Ico-TRI;

The area ratio of the superimposed area and the shape similarity vector constitute the similarity measure vector (ε,σ _Gtri→Ghex ,σ _Ghex→Gtri ).

Further, the specific sub-steps of the heterogeneous grid S1 are as follows:

和/>

S11: Since the GeoSOT grid and Ico-TRI are both equal-area grids, count the number of elements in G _geosot and G _tri respectively, and according to the resolution, the area of geographic elements expressed by the two grid systems can be obtained

and />

S12: The elements in GeoSOT grid and Ico-TRI are superimposed, which is equivalent to the intersection of spherical polygons. First calculate the intersecting intersection point, connect the intersection points sequentially to obtain the spherical polygon, use the theory of spherical geometry to directly calculate the area of the intersecting spherical polygon, get

S13: Calculate the area ratio of the superposition area, the calculation formula is

Further, the specific sub-steps of the heterogeneous grid S2 are as follows:

S21: Calculate the vertical distance of the boundary formed by the center point G _geosotc of each element corresponding to GeoSOT to the center point of G _tri in the Ico-TRI grid system, and obtain the vertical distance set d of the boundary grid from GeoSOT to Ico-TRI _{Ggeosot→Gtri} , where d represents the vertical distance of the boundary grid, and the subscript represents the conversion from GeoSOT to Ico-TRI;

S22: Using the same method, obtain the vertical distance set d _{Gtri→Ggeosot} of the boundary grid from Ico-TRI to GeoSOT, where d represents the vertical distance of the boundary grid, and the subscript represents the conversion from Ico-TRI to GeoSOT;

S23: Calculate the variance of d _{Ggeosot→Gtri} and d _{Gtri→Ggeosot} respectively to obtain σ _{Ggeosot→Gtri} and σ _{Gtri→Ggeosot} to form a shape similarity vector (σ _{Ggeosot→Gtri} , σ _{Gtri→Ggeosot} ).

Further, the specific sub-steps of the isomorphic grid S2 are as follows:

和/>

S21: Count the number of grids of geographic elements in the two grid systems, N and M respectively, and the area of each grid is

and />

其中s_Gtri和s_Ghex是单个三角形格网和六边形格网的面积，/>

表示G_tri和G_hex的叠加面积；S22: For point elements, the superimposed area in a single grid G _tri and G _hex

where s _Gtri and s _Ghex are the areas of individual triangular and hexagonal grids, />

Indicates the superposition area of G _tri and G _hex ;

其中n_i表示每个六边形出现的次数，n_i＞3；同样的，在Ico-HEX→Ico-TRI转换过程中，统计三角形出现的次数，每个三角形最多出现三次。S23: The area of the overlapping area of line elements. According to the conversion rules, one triangular grid corresponds to three hexagonal grids. During the Ico-TRI→Ico-HEX conversion process, record the number of hexagonal occurrences. When the hexagonal When the cumulative number exceeds 3, it proves that there are more than three triangles corresponding to the hexagon, record the hexagon and its occurrence times, and use it to calculate the intersection area. The covered area is

Where n _i represents the number of occurrences of each hexagon, and n _i >3; similarly, in the Ico-HEX→Ico-TRI conversion process, the number of occurrences of triangles is counted, and each triangle appears three times at most.

S24: The area of the overlapping area of the area element is determined by the area of the non-overlapping area, because there is no starting point and end point grid for the area element, the polyline is used as the line element, and the calculation method of the area of the overlapping area of the line element is used.

S25: Calculate the area ratio of the superposition area, the calculation formula is

Further, the specific sub-steps of the isomorphic grid S3 are as follows:

S31: According to the dual relationship between the triangular and hexagonal grids, calculate the center point distance of Ico-TRI and Ico-HEX, as the maximum tolerable interoperability offset of point elements, there is d _Gtri→Ghex =d _Ghex→Gtri , Where d represents the vertical distance of the boundary grid, and the subscript represents the conversion from Ico-TRI to Ico-HEX and from Ico-HEX to Ico-TRI;

S32: Calculate the offset of each grid in the line element grid, calculate the variance of d _Gtri→Ghex and d _Ghex→Gtri respectively, and obtain σ _Gtri→Ghex and σ _Ghex→Gtri to form a shape similarity vector (σ _{Gtri →Ghex} , σGhex _→Gtri ).

S33: Area features are composed of boundary grids and internal grids. The calculation of shape similarity depends on boundary grids, so the calculation method is the same as that of line features.

Compared with the prior art, the present invention has the advantages of:

1. The grid interoperability quality evaluation is classified and processed, the heterogeneous grid calculates the overlapping area to measure the position and geometry, and calculates the boundary grid distance to obtain the shape similarity vector; the homogeneous grid uses the geometric characteristics of the grid itself And measure the dual relationship.

2. Different from the traditional vector-based evaluation method, the homogeneous grid interoperability quality evaluation method can effectively reduce the amount of calculation and improve the calculation speed; the overlapping area of heterogeneous grids provides a preliminary accuracy judgment for the grid interoperability quality. Lay the foundation for subsequent high-precision calculations.

Description of drawings

Fig. 1 is the flow chart of the DGGS interoperability quality evaluation method for geographic full elements in an embodiment of the present invention;

Fig. 2 is the schematic diagram of the expression of geographical full elements in different grids in the embodiment of the present invention, wherein Fig. 2 (a) is the geographical elements expressed by GeoSOT, Fig. 2 (b) is the geographical elements expressed by Ico-TRI, and Fig. 2 (c ) is the geographical element expressed by Ico-HEX;

Figure 3 shows the dual relationship of the isomorphic grid;

Figure 4 is a schematic diagram of the overlay of geographical elements of the isomorphic grid, where Figure 4 (a) is the overlay process of Ico-HEX to Ico-TRI conversion, and Figure 4 (b) is the overlay process of Ico-TRI to Ico-HEX conversion;

5 is a schematic diagram of calculating shape similarity according to an embodiment of the present invention;

Figure 6 is a comparison of the efficiency of isomorphic grids and traditional methods.

Detailed ways

In order to make the object, technical solution and advantages of the present invention clearer, the present invention will be further described in detail below with reference to the accompanying drawings and examples.

According to different interoperable grid systems, grid interoperability can be divided into heterogeneous grid interoperability and homogeneous grid interoperability. Heterogeneous grid interoperability refers to the interoperability between grids with different structures, such as GeoSOT grid and spherical icosahedral triangular grid (Ico-TRI); Interoperability between structural grids, such as spherical icosahedral hexagonal grid (Icosahedral hexagonal grid, Ico-HEX).

For the interoperable quality evaluation method of heterogeneous grids, firstly, geographical elements are expressed according to the rules under the two grid systems. The feature expression in the GeoSOT grid corresponds the geographic coordinates of the point, line and area elements to the GeoSOT grid, as shown in Figure 2(a). The representation of elements in the Ico-TRI grid is shown in Figure 2(b). The spherical coordinates of the points are known, and the Snyder projection is used to project the point elements onto the icosahedron datum, and the grid code is determined according to the grid resolution to realize the grid expression of the point elements. For line features, the coordinates of the points that make up the line feature are known, starting from the starting grid, searching for adjacent grids, finding the intersecting grid, using the currently determined grid as the starting grid, continuing to search until the end point, and realizing the line feature grid expression. For surface features, first determine the grid in the MRB of the surface feature, and then further determine the grid in the surface feature and the grid intersecting with the boundary of the surface feature to realize the grid expression of the surface feature.

Calculate the similarity measure vector according to the gridded result of the element, define the gridded result of the element as grid G _geosot and grid G _tri respectively, as shown in Figure 1, the specific steps for calculating the similarity measure vector are as follows:

S1: Calculation of the area ratio of the superposition area;

和/>

具体步骤如下：According to spherical geometry theory, calculate the intersection area of GeoSOT grid G _geosot and Ico-TRI grid G _tri

and />

Specific steps are as follows:

和/>

S11: Since the GeoSOT grid and Ico-TRI are equal-area grids, count the number of elements in G _geosot and G _tri respectively, and according to the resolution, we can get

and />

S12: The elements in GeoSOT grid and Ico-TRI are superimposed, which is equivalent to the intersection of spherical polygons. First calculate the intersecting intersection point, connect the intersection points sequentially to obtain the spherical polygon, use the theory of spherical geometry to directly calculate the area of the intersecting spherical polygon, get

S13: Calculate the area ratio of the superposition area, the calculation formula is

S2: shape similarity;

Before the calculation, determine the boundary grids in the grids G _geosot and G _tri respectively, and obtain the boundary grid sets B _geosot and B _tri , where the point elements have no boundary grids, and the line elements are all boundary grids. For the boundary grid The set calculates the shape similarity, the specific steps are as follows:

S21: Calculate the vertical distance in Ico-TRI corresponding to the boundary of each element center point G _Trigeosotc in GeoSOT to the element center point in G _tri , and obtain the boundary grid vertical distance set d _Ggeosot→ from GeoSOT to Ico-TRI _Gtri , as shown in Figure 5;

S22: Using the same method, obtain the vertical distance set d _{Gtri→Ggeosot} of the boundary grid from Ico-TRI to GeoSOT;

S23: Calculate the variance of d _{Ggeosot→Gtri} and d _{Gtri→Ggeosot} respectively to obtain σ _{Ggeosot→Gtri} and σ _{Gtri→Ggeosot} to form a shape similarity vector (σ _{Ggeosot→Gtri} , σ _{Gtri→Ggeosot} ).

S3: construction of similarity measure vector;

The area ratio of the superposition area and the shape similarity vector constitute the similarity measure vector (ε,σ _{Ggeosot→Gtri} , σ _{Gtri→Ggeosot} ).

For the interoperable quality evaluation method of isomorphic grids, firstly, geographical elements are expressed in different grid systems with known resolutions. For the expression of elements in Ico-HEX grid, refer to the expression method of Ico-TRI grid elements. The results are shown in the figure 2(c).

Calculate the similarity measure vector according to the result of element gridding, and define the result of element gridding as grid G _tri and grid G _hex respectively. The specific steps to calculate the similarity measure vector are as follows:

S1: Define the dual relationship between triangular grid and hexagonal grid;

As shown in Figure 3, the duality relationship between the triangle and the hexagon is that the vertex of the triangle corresponds to the center of the hexagon, the center of the triangle corresponds to the vertex of the hexagon, and the three corners of the triangle are marked as a vertex and The two base corners and the point specifying the top corner correspond to the center of the hexagon.

S2: similarity of geographic element area;

According to the dual relationship between triangular and hexagonal grids, the overlapping area of point elements can be directly obtained according to the dual relationship, the overlapping area of line elements depends on the direction of the G _tri grid, and the overlapping area of area elements depends on the boundary grid.

和/>

S21: Count the number of grids of geographic elements in the two grid systems, N and M respectively, and the area of each grid is

and />

其中s_Gtri和s_Ghex是单个三角形格网和六边形格网的面积；S22: For point elements, the superimposed area in a single grid G _tri and G _hex

where s _Gtri and s _Ghex are the areas of individual triangular and hexagonal grids;

其中n_i表示每个六边形出现的次数，n_i＞3，具体叠加情况如图4(a)所示；同样的，在Ico-HEX→Ico-TRI转换过程中，统计三角形出现的次数，每个三角形最多出现三次，具体叠加情况如图4(b)所示。S23: The area of the overlapping area of line elements. According to the conversion rules, one triangular grid corresponds to three hexagonal grids. During the Ico-TRI→Ico-HEX conversion process, record the number of hexagonal occurrences. When the hexagonal When the cumulative number exceeds 3, it proves that there are more than three triangles corresponding to the hexagon, record the hexagon and its occurrence times, and use it to calculate the intersection area. The covered area is

Among them, n _i represents the number of occurrences of each hexagon, and n _i > 3, the specific superposition situation is shown in Figure 4(a); similarly, in the process of Ico-HEX→Ico-TRI conversion, count the number of occurrences of triangles , each triangle appears at most three times, and the specific superposition situation is shown in Figure 4(b).

S24: The area of the overlapping area of the area element is determined by the area of the non-overlapping area, because there is no starting point and end point grid for the area element, the polyline is used as the line element, and the calculation method of the area of the overlapping area of the line element is used.

S25: Calculate the area ratio of the superposition area, the calculation formula is

S3: geographic feature shape similarity;

Before calculation, determine the boundary grids in the grids G _tri and G _hex respectively, and obtain the boundary grid sets B _tri and B _hex , in which the point elements have no boundary grids, all the line elements are boundary grids, and some area elements exist Boundary grid, calculating the shape similarity for the boundary grid set;

S31: According to the dual relationship between the triangular and hexagonal grids, calculate the center point distance of Ico-TRI and Ico-HEX, as the maximum tolerable interoperability offset of point elements, there is _dGtri→Ghex = _dGhex→Gtri ;

S32: Calculate the offset of each grid in the line element grid, calculate the variance of d _Gtri→Ghex and d _Ghex→Gtri respectively, and obtain σ _Gtri→Ghex and σ _Ghex→Gtri to form a shape similarity vector (σ _{Gtri →Ghex} , σGhex _→Gtri ).

S33: Surface elements are composed of boundary grids and internal grids. The calculation of shape similarity depends on boundary grids, so the calculation method is the same as that of line elements.

S4: construction of similarity measure vector;

The area ratio of the superimposed area and the shape similarity vector constitute the similarity measure vector (ε,σ _Gtri→Ghex ,σ _Ghex→Gtri ).

1:10000 map of a certain area is selected, and after comprehensive cartographic processing, 5000 point elements, 5000 line elements and 5000 area elements are respectively selected to compare the calculation efficiency with the traditional vector-based method in the isomorphic grid system. The result is shown in Figure 6. It can be seen from the figure that the grid quality evaluation efficiency of the present invention is higher than that of the traditional method.

Those skilled in the art will appreciate that the embodiments described here are to help readers understand the implementation method of the present invention, and it should be understood that the protection scope of the present invention is not limited to such specific statements and embodiments. Those skilled in the art can make various other specific modifications and combinations based on the technical revelations disclosed in the present invention without departing from the essence of the present invention, and these modifications and combinations are still within the protection scope of the present invention.

Claims (5)

1. A DGGS interoperability quality evaluation method oriented to geographic full elements, characterized in that, the grid interoperability quality evaluation is classified and processed, the heterogeneous grid calculates the overlapping area to measure the position and geometry, and calculates the boundary grid The shape similarity vector is obtained from the distance; the homogeneous grid uses the geometric characteristics and duality of the grid itself to measure; for the heterogeneous grid, firstly, the geographical elements are expressed according to the rules under the two grid systems, and the GeoSOT grid Element expression, corresponding the geographic coordinates of point, line and area elements to the GeoSOT grid; element expression in the Ico-TRI grid, according to the code conversion of the regular icosahedral triangular grid and hexagonal grid based on the Sndyer projection Rules to express geographic elements separately; for point elements, the spherical coordinates of points are known, and the Snyder projection is used to project the point elements onto the icosahedron datum plane, and the grid code is determined according to the grid resolution to realize the grid expression of point elements ;For the line element, the point coordinates of the line element are known, start from the starting grid, search the adjacent grid, find the intersecting grid, use the currently determined grid as the starting grid, continue to search, and end at the end point to realize the line Element grid expression; for surface elements, first determine the grid in the area element MRB, and then further determine the grid inside the area element and the grid intersecting with the boundary of the area element to realize the grid expression of the area element; For heterogeneous grids, the similarity measure vector is calculated according to the gridded results of the elements, and the gridded results of the elements are defined as the grid G _geosot and the grid G _tri , where G _geosot represents the grid set that expresses the geographical features by GeoSOT, G _tri represents the grid collection of geographical features represented by Ico-TRI; For isomorphic grids, geographic elements are first expressed in different grid systems with known resolutions. For the expression of elements in the Ico-HEX grid, refer to the Ico-TRI grid element expression method; the similarity is calculated according to the gridded results of the elements Measure the vector, define the result of element gridding as grid G _tri and grid G _hex respectively, and G _hex represents the set of grids expressing geographical elements by Ico-HEX; The specific steps of calculating the similarity measure vector in the heterogeneous grid are as follows: S1: Calculation of the area ratio of the superposition area; according to the theory of spherical geometry, calculate the intersection area of the GeoSOT grid G _geosot and the Ico-TRI grid G _tri

and />

Where S represents the area, and the subscripts correspond to the intersection area of G _geosot and G _tri , the grid area of G _geosot and the grid area of G _tri , respectively; S2: shape similarity; Before the calculation, determine the boundary grids in the grids G _geosot and G _tri respectively, and obtain the boundary grid sets B _geosot and B _tri , where the point elements have no boundary grids, and the line elements are all boundary grids. For the boundary grid The shape similarity is calculated by the set, where B represents the boundary grid set, and the subscripts correspond to the GeoSOT and Ico-TRI grid systems respectively; S3: construction of similarity measure vector; The area ratio of the superimposed area and the shape similarity vector constitute the similarity measure vector (ε,σ _{Ggeosot→Gtri} , σ _{Gtri→Ggeosot} ), where ε represents the area ratio of the superimposed area, σ represents the distance variance, and the subscripts respectively represent GeoSOT→Ico-TRI and conversion of Ico-TRI→GeoSOT; The specific steps of calculating the similarity measure vector in the isomorphic grid are as follows: The specific steps to calculate the similarity measure vector in the isomorphic grid are as follows: S1: Define the dual relationship between triangular grid and hexagonal grid; The duality relationship between triangle and hexagon is as follows: the apex of the triangle corresponds to the center of the hexagon, and the center of the triangle corresponds to the apex of the hexagon, and the three corners of the triangle are marked as one apex and two base angles. The point of the corner corresponds to the center of the hexagon; S2: similarity of geographic element area; According to the duality relationship between triangular and hexagonal grids, the overlapping area of point elements can be directly obtained according to the duality relationship, the overlapping area of line elements depends on the direction of the G _tri grid, and the overlapping area of area elements depends on the boundary grid; S3: geographic feature shape similarity; Before calculation, determine the boundary grids in the grids G _tri and G _hex respectively, and obtain the boundary grid sets B _tri and B _hex , in which the point elements have no boundary grids, all the line elements are boundary grids, and some area elements exist Boundary grid, which calculates the shape similarity for the boundary grid set, where B _hex represents the boundary grid in the Ico-HEX expression geographic element grid set G _hex ; S4: construction of similarity measure vector; The area ratio of the superimposed area and the shape similarity vector form a similarity measure vector (ε,σ _Gtri→Ghex ,σ _Gtri→Ghex ), where σ is the distance variance, and the subscripts represent the distance from Ico-TRI to Ico-HEX and Ico-HEX Conversion to Ico-TRI; The area ratio of the superimposed area and the shape similarity vector constitute the similarity measure vector (ε,σ _Gtri→Ghex ,σ _Ghex→Gtri ). 2. a kind of DGGS interoperability quality evaluation method oriented to geographic full element according to claim 1, is characterized in that, the concrete sub-step of heterogeneous grid S1 superposition area ratio calculation is as follows: S11: Since the GeoSOT grid and Ico-TRI are both equal-area grids, count the number of elements in G _geosot and G _tri respectively, and according to the resolution, the area of geographic elements expressed by the two grid systems can be obtained

and />

S12: The elements in GeoSOT grid and Ico-TRI are superimposed, which is equivalent to intersecting spherical polygons; first, calculate the intersection point, connect the intersection points sequentially, and obtain the spherical polygon, and use the theory of spherical geometry to directly calculate the area of the intersecting spherical polygon.

S13: Calculate the area ratio of the superposition area, the calculation formula is

3. a kind of DGGS interoperability quality evaluation method oriented to geographic full element according to claim 1, is characterized in that, the concrete sub-step of heterogeneous grid S2 shape similarity calculation is as follows: S21: Calculate the vertical distance of the boundary formed by the center point G _geosotc of each element corresponding to GeoSOT to the center point of G _tri in the Ico-TRI grid system, and obtain the vertical distance set d of the boundary grid from GeoSOT to Ico-TRI _{Ggeosot→Gtri} , where d represents the vertical distance of the boundary grid, and the subscript represents the conversion from GeoSOT to Ico-TRI; S22: Using the same method, obtain the vertical distance set d _{Gtri→Ggeosot} of the boundary grid from Ico-TRI to GeoSOT, where d represents the vertical distance of the boundary grid, and the subscript represents the conversion from Ico-TRI to GeoSOT; S23: Calculate the variance of d _{Ggeosot→Gtri} and d _{Gtri→Ggeosot} respectively to obtain σ _{Ggeosot→Gtri} and σ _{Gtri→Ggeosot} to form a shape similarity vector (σ _{Ggeosot→Gtri} , σ _{Gtri→Ggeosot} ). 4. a kind of DGGS interoperability quality evaluation method oriented to geographical full element according to claim 1, is characterized in that, the concrete sub-step of isomorphic grid S2 geographic element area similarity calculation is as follows: S21: Count the number of grids of geographic elements in the two grid systems, N and M respectively, and the area of each grid is

and />

S22: For point elements, the superimposed area in a single grid G _tri and G _hex

where s _Gtri and s _Ghex are the areas of individual triangular and hexagonal grids, />

Indicates the superposition area of G _tri and G _hex ; S23: The area of the overlapping area of line elements. According to the conversion rules, one triangular grid corresponds to three hexagonal grids. During the Ico-TRI→Ico-HEX conversion process, record the number of hexagonal occurrences. When the hexagonal When the cumulative number exceeds 3, it proves that there are more than three triangles corresponding to the hexagon, record the hexagon and its occurrence times, and use it to calculate the intersection area. The covered area is

Among them, n _i represents the number of occurrences of each hexagon, and n _i >3; similarly, during the Ico-HEX→Ico-TRI conversion process, the number of occurrences of triangles is counted, and each triangle appears up to three times; S24: The area of the overlapping area of the area element is determined by the area of the non-overlapping area, because there is no starting point and end point grid for the area element, the polyline is used as the line element, and the calculation method of the area of the overlapping area of the line element is adopted; S25: Calculate the area ratio of the superposition area, the calculation formula is

5. a kind of DGGS interoperability quality evaluation method oriented to geographical full element according to claim 1, is characterized in that, the concrete sub-step of isomorphic grid S3 geographic element shape similarity calculation is as follows: S31: According to the dual relationship between the triangular and hexagonal grids, calculate the center point distance of Ico-TRI and Ico-HEX, as the maximum tolerable interoperability offset of point elements, there is d _Gtri→Ghex =d _Ghex→Gtri , Where d represents the vertical distance of the boundary grid, and the subscript represents the conversion from Ico-TRI to Ico-HEX and from Ico-HEX to Ico-TRI; S32: Calculate the offset of each grid in the line element grid, calculate the variance of d _Gtri→Ghex and d _Ghex→Gtri respectively, and obtain σ _Gtri→Ghex and σ _Ghex→Gtri to form a shape similarity vector (σ _{Gtri → Ghex} , σ _{Ghex → Gtri} ); S33: Area features are composed of boundary grids and internal grids. The calculation of shape similarity depends on boundary grids, so the calculation method is the same as that of line features.

Images