A grid coding method and computer system

Abstract

The embodiment of the present application discloses a grid coding method, including: obtaining a target projection, the target projection is a projection of the earth's surface on a plane, and the area difference between the two projection areas corresponding to any two areas of the same area on the earth's surface on the target projection is less than a preset value; according to the target projection, determining a square plane, the area where the target projection is located is covered by the square plane; bisecting the square plane multiple times in the horizontal and vertical directions to obtain multiple square grids; encoding the multiple square grids to obtain multiple codes, each code corresponding to a square grid, and each code is used to indicate the area covered by the corresponding square grid. The present application performs grid division on the square plane constructed on the target projection, so that the division size of the earth's surface corresponding to the square grids obtained by the same level of division is basically the same, which is conducive to calculations sensitive to the size of grid data.

Description

Technical Field

[0001] The present application relates to the field of computers, and in particular to a grid coding method and a computer system. Background Art

[0002] The discrete global grid (DGG) is a kind of earth body fitting grid based on the sphere (or ellipsoid) that can be infinitely subdivided without changing its shape. When subdivided to a certain extent, it can achieve the purpose of simulating the earth's surface. DGG has the characteristics of hierarchy and global continuity, overcoming many constraints and uncertainties that limit the application of geographic information systems, so that any spatial data of any resolution (different accuracy) obtained at any location on the earth can be expressed and analyzed in a standardized manner, and multi-scale operations can be performed with a certain accuracy.

[0003] DGG includes GeoSOT grids, which belong to the quadtree grid system with equal longitude and latitude. First, the space on the earth's surface is expanded into a square plane space. The GeoSOT level 0 grid is defined as a plane space grid centered on the intersection of the equator and the prime meridian. The GeoSOT level 1 grid is evenly divided into four equal parts based on level 0. The subsequent subdivision levels are based on the above quadtree principle and so on.

[0004] However, since GeoSOT belongs to a quadtree grid system with equal longitude and latitude, the Euclidean distance of the grids is different at each level, and the size of the grids gradually decreases from the equator to the poles, making it difficult to generalize calculations that are sensitive to grid size. Summary of the invention

[0005] In a first aspect, the present application provides a grid coding method, the method comprising:

[0006] Obtain a target projection, wherein the target projection is a projection of the earth's surface on a plane, and the area difference between the two projection areas corresponding to any two areas of the same area on the earth's surface on the target projection is less than a preset value; the earth's surface can be converted to a plane using a map projection method to obtain the target projection, wherein the map projection method can be an overall projection or a zonal projection. The earth's surface plane of equal longitude and latitude is obtained by reducing and deforming the earth's surface by a certain proportion. On the earth's surface plane of equal longitude and latitude, the reduction ratio of the area far from the equator is much smaller than the reduction ratio of the area close to the equator, that is, the area difference between the two areas corresponding to the area far from the equator and the area close to the equator on the earth's surface plane of equal longitude and latitude is large. In this embodiment, the area difference between the two projection areas corresponding to any two areas of the same area on the earth's surface on the target projection is less than a preset value, that is, any two areas of the same area on the earth's surface are reduced by a certain proportion to obtain the target projection, and the reduction ratio difference between any two areas of the same area on the earth's surface is very small. The present application does not limit the area between the two projection areas corresponding to any two areas of the same area on the earth's surface on the target projection to be completely consistent, but the difference is less than a preset value. According to the target projection, a square plane is determined, and the area where the target projection is located is covered by the square plane; after obtaining the target projection, a square plane can be determined according to the target projection, and the area where the target projection is located is covered by the square plane. In order to completely cover the area where the target projection is located, the area of ​​the square plane should be larger than the target projection, and the side of the square plane should be outside the area of ​​the target projection. It should be understood that the target projection is non-square, so the square plane includes areas overlapping with the target projection and areas not overlapping with the target projection. The square plane is bisected multiple times in the horizontal and vertical directions to obtain multiple square grids; each time the square plane is bisected in the horizontal and vertical directions, a group of square grids can be obtained. When bisecting the square plane in the next bisecting, each square grid in the group of square grids obtained by the previous bisecting can be bisected in the horizontal and vertical directions to obtain a new group of square grids. The multiple square grids are encoded to obtain multiple codes, each code corresponds to a square grid, and each code is used to indicate an area covered by the corresponding square grid.After the square plane is bisected in the horizontal and vertical directions, a group of square grids can be obtained, and then the group of square grids can be encoded to obtain the code of the group of square grids, and then the next bisection can be performed, that is, each square grid in the group of square grids obtained by the previous division can be bisected in the horizontal and vertical directions to obtain a new group of square grids, and then the group of square grids can be encoded to obtain the code of the group of square grids, and so on.

[0007] It should be understood that each code obtained after encoding the square grid is used to indicate the area covered by the corresponding square grid, and the present application does not limit the specific encoding method.

[0008] This embodiment performs grid division on the square plane constructed on the target projection, so that the sizes of the earth surface divisions corresponding to the square grids obtained by the same level of division are basically consistent, which is beneficial to calculations that are sensitive to the size of grid data.

[0009] In a possible implementation, the first target point and the second target point are any two points on the target projection, the first target point is the projection of the first position point on the earth surface, the second target point is the projection of the second position point on the earth surface, the distance between the first target point and the second target point is the first distance, the distance between the first position point and the second position point is the second distance, and the ratio between the first distance and the second distance is within a preset range. The target projection can more accurately characterize the Euclidean distance of each point on the earth surface. Specifically, the distance between any two points on the target projection can characterize the Euclidean distance between two physical position points on the corresponding earth surface. The so-called Euclidean distance refers to the Euclidean metric. In mathematics, the Euclidean distance or Euclidean metric is the "ordinary" (i.e., straight line) distance between two points in the Euclidean space. The target projection is an approximate undistorted projection of the earth surface. The distance between each point on the target projection is obtained by reducing the earth surface by a certain proportion, and the reduction proportion of the distance between each point on the earth surface is basically the same. Therefore, the target projection can more accurately characterize the Euclidean distance of each point on the earth surface.

[0010] In a possible implementation, the earth's surface includes an equatorial latitude and a target longitude perpendicular to the equatorial latitude, the target projection includes a projected equatorial latitude and a projected target longitude, the horizontal axis direction of the square plane is consistent with the direction of the projected equatorial latitude, and the vertical axis direction of the square plane is consistent with the direction of the projected target longitude. First, it is necessary to determine the direction of the square plane, and the so-called direction refers to the direction of the horizontal axis and the vertical axis included in the square plane. The square plane may include a horizontal axis and a vertical axis, wherein the horizontal axis may be parallel to one side of the square plane, and the vertical axis is perpendicular to the horizontal axis. In more detail, the horizontal axis may be an axis parallel to one side of the square plane and passing through the center of the square plane, and the vertical axis is an axis perpendicular to the horizontal axis and passing through the center of the square plane.

[0011] In a possible implementation, determining a square plane according to the target projection includes:

[0012] The side length of the square plane is determined according to the target length of the target projected in the direction of the equatorial latitude after the projection, wherein the side length of the square plane is greater than or equal to the target length, and the side length of the square plane is a preset multiple of an integer power of 2; the reason why the side length of the square plane is a preset multiple of an integer power of 2 is that when performing subsequent encoding, the square plane needs to be divided into two equal parts in the horizontal and vertical directions. In order to obtain an integer power of 2 square grid after the division, it is necessary to ensure that the side length of the square plane is a preset multiple of an integer power of 2. For example, the side length of the square plane can be Wherein, L is the length of the equatorial latitude in the target projection, and S is the scaling factor, that is, the preset multiple mentioned above. According to the side length of the square plane, the square plane is determined so that the area where the target projection is located is covered by the square plane. In this embodiment, the size of the square grid obtained after the subdivision can be controlled by controlling the size of the preset multiple.

[0013] In a possible implementation, determining the side length of the square plane according to the target length of the target projected in the direction where the equatorial latitude line is located after the projection includes:

[0014] The preset multiple is obtained, and the side length of the square plane is determined according to the target length of the target projected in the direction of the equatorial latitude after the projection and the preset multiple. In a possible implementation, the central horizontal axis of the square plane overlaps with the equatorial latitude after the projection.

[0015] In a possible implementation, the square plane is bisected multiple times in the horizontal and vertical directions to obtain multiple square grids, including:

[0016] The square plane is bisected multiple times in the horizontal and vertical directions to obtain multiple groups of square grids, each group of square grids includes multiple square grids, wherein the multiple bisecting includes the Nth bisecting and the N+1th bisecting, and after the Nth bisecting, M square grids are obtained, and after the N+1th bisecting, M*4 square grids are obtained, and the M*4 square grids are obtained by bisecting each square grid in the M square grids in the horizontal and vertical directions;

[0017] The encoding of the plurality of square grids to obtain a plurality of codes includes: encoding a plurality of square grids included in each group of square grids in the plurality of groups of square grids to obtain a plurality of groups of codes, each group of codes including a plurality of codes.

[0018] In a possible implementation, the target projection is obtained by performing an overall projection or a zone projection on the earth's surface.

[0019] In a second aspect, the present application provides a grid coding device, the device comprising:

[0020] An acquisition module, used for acquiring a target projection, where the target projection is a projection of the earth's surface on a plane;

[0021] A determination module, used for determining a square plane according to the target projection, wherein the area where the target projection is located is covered by the square plane;

[0022] A grid division module, used for dividing the square plane into two equal parts multiple times in the horizontal axis direction and the vertical axis direction to obtain multiple square grids;

[0023] The encoding module is used to encode the multiple square grids to obtain multiple codes, each code corresponds to a square grid, and each code is used to indicate the area covered by the corresponding square grid.

[0024] In a possible implementation, the first target point and the second target point are any two points on the target projection, the first target point is the projection of a first position point on the earth's surface, the second target point is the projection of a second position point on the earth's surface, the distance between the first target point and the second target point is a first distance, the distance between the first position point and the second position point is a second distance, and the ratio between the first distance and the second distance is within a preset range.

[0025] In one possible implementation, the earth's surface includes an equatorial latitude and a target longitude perpendicular to the equatorial latitude, the target projection includes a projected equatorial latitude and a projected target longitude, the horizontal axis direction of the square plane is consistent with the direction of the projected equatorial latitude, and the vertical axis direction of the square plane is consistent with the direction of the projected target longitude.

[0026] In a possible implementation, the determination module is used to determine the side length of the square plane according to the target length of the target projected in the direction of the equatorial latitude after the projection, wherein the side length of the square plane is greater than or equal to the target length, and the side length of the square plane is a preset multiple of an integer power of 2;

[0027] The square plane is determined according to the side length of the square plane, so that the area where the target projection is located is covered by the square plane.

[0028] In a possible implementation, the determination module is used to obtain the preset multiple, and determine the side length of the square plane according to the target length of the target projected in the direction of the equatorial latitude after the projection and the preset multiple.

[0029] In a possible implementation, the central horizontal axis of the square plane overlaps with the projected equatorial latitude.

[0030] In a possible implementation, the grid division module is used to perform multiple bisection divisions on the square plane in the horizontal and vertical directions to obtain multiple groups of square grids, each group of square grids includes multiple square grids, wherein the multiple bisection divisions include the Nth bisection division and the N+1th bisection division, and after the Nth bisection division, M square grids are obtained, and after the N+1th bisection division, M*4 square grids are obtained, and the M*4 square grids are obtained by bisectioning each square grid in the M square grids in the horizontal and vertical directions;

[0031] The encoding module is used to encode the multiple square grids included in each group of square grids in the multiple groups of square grids to obtain multiple groups of codes, each group of codes including multiple codes.

[0032] In a possible implementation, the target projection is obtained by performing an overall projection or a zone projection on the earth's surface.

[0033] In a third aspect, the present application provides a computer system, which includes a memory and a processor, wherein the memory is used to store computer-readable instructions (or computer programs), and the processor is used to read the computer-readable instructions to implement the method provided by any of the aforementioned implementations.

[0034] In a fourth aspect, the present application provides a computer storage medium, which may be non-volatile. The computer storage medium stores computer-readable instructions, which, when executed by a processor, implement the method provided by any of the aforementioned implementations.

[0035] In a fifth aspect, the present application provides a computer program product, which includes computer-readable instructions. When the computer-readable instructions are executed by a processor, the method provided by any of the aforementioned implementations is implemented.

[0036] The embodiment of the present application provides a grid encoding method, the method comprising: obtaining a target projection, the target projection being a projection of the earth's surface on a plane, and the area difference between the two projection areas corresponding to any two areas of the same area on the earth's surface on the target projection is less than a preset value; determining a square plane according to the target projection, the area where the target projection is located being covered by the square plane; bisecting the square plane multiple times in the horizontal and vertical directions to obtain a plurality of square grids; encoding the multiple square grids to obtain a plurality of codes, each code corresponding to a square grid, and each code being used to indicate the area covered by the corresponding square grid. In the above manner, grid division is performed on the square plane constructed on the target projection, so that the division sizes of the earth's surface corresponding to the square grids obtained by the same level of division are substantially consistent, which is beneficial to calculations that are sensitive to the size of the grid data. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] Figure 1 A schematic diagram of a GeoSOT grid;

[0038] Figure 2 A schematic diagram of a GeoSOT grid;

[0039] Figure 3 A schematic diagram of a grid coding method provided in an embodiment of the present application;

[0040] Figure 4 A schematic diagram of a target projection provided in an embodiment of the present application;

[0041] Figure 5 A schematic diagram of a square plane provided in an embodiment of the present application;

[0042] Figure 6A schematic diagram of the division of a square grid provided in an embodiment of the present application;

[0043] Figure 7 A schematic diagram of the division of a square grid provided in an embodiment of the present application;

[0044] Figure 8 A schematic diagram of the structure of a grid coding device provided in an embodiment of the present application;

[0045] Fig. 9 A schematic diagram of the structure of a computer system provided in this embodiment;

[0046] Fig.10 A schematic diagram of the structure of an NPU provided in this embodiment. DETAILED DESCRIPTION

[0047] The following describes the embodiments of the present application in conjunction with the accompanying drawings. Obviously, the described embodiments are only embodiments of a part of the present application, rather than all embodiments. It is known to those skilled in the art that with the development of technology and the emergence of new scenarios, the technical solutions provided in the embodiments of the present application are also applicable to similar technical problems.

[0048] The terms "first", "second", etc. in the specification and claims of this application and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchangeable where appropriate, so that the embodiments described here can be implemented in an order other than that illustrated or described here. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or modules is not necessarily limited to those steps or modules clearly listed, but may include other steps or modules that are not clearly listed or inherent to these processes, methods, products or devices. The naming or numbering of steps in this application does not mean that the steps in the method flow must be executed in the time / logical sequence indicated by the naming or numbering. The process steps that have been named or numbered can change the execution order according to the technical purpose to be achieved, as long as the same or similar technical effects can be achieved.

[0049] In addition, in the description of the present invention, unless otherwise specified, "plurality" means two or more. The term "and / or" or the character " / " in this application is only a description of the association relationship of associated objects, indicating that there may be three relationships, for example, A and / or B, or A / B, can mean: A exists alone, A and B exist at the same time, and B exists alone.

[0050] Discrete Global Grid (DGG) is a kind of earth body fitting grid based on sphere (or ellipsoid) that can be infinitely subdivided without changing its shape. When subdivided to a certain extent, it can achieve the purpose of simulating the earth's surface. DGG has the characteristics of hierarchy and global continuity, which overcomes many constraints and uncertainties that limit the application of geographic information systems, so that spatial data of any resolution (different accuracy) obtained at any location on the earth can be expressed and analyzed in a standardized manner, and multi-scale operations can be performed with a certain accuracy.

[0051] According to different methods of constructing grids, the global discrete grid system can be roughly divided into four categories: equal longitude and latitude global grid, variable longitude and latitude global grid, adaptive global grid and regular polyhedron global grid system. For example, GeoSOT is an equal longitude and latitude global grid, the degenerate quadtree grid is a variable longitude and latitude global grid, the levels of detail (LOD) model of the digital elevation model (DEM) data is an adaptive global grid, and the triangular grid, diamond grid, hexagonal grid, etc. are regular polyhedron global grids.

[0052] The GeoSOT grid belongs to the quadtree grid system with equal longitude and latitude. First, the 360°*180° space on the earth's surface is expanded to a 512°*512° plane space. The GeoSOT level 0 grid is defined as a 512°*512° plane space grid centered on the intersection of the equator and the prime meridian. The GeoSOT level 1 grid is evenly divided into four equal parts based on level 0. The size of each level 1 grid is 256°*256°. Figure 1 As shown. GeoSOT's 2-level grid is divided into four equal parts based on level 1. The size of each level 1 grid is 128°*128°. The following subdivision levels are deduced according to the principle of the above quadtree. GeoSOT's 9-level grid is 1°*1° in size. Grids above level 9 are GeoSOT's "degree" level grids. Grids from level 10 to 15 are "sub" level grids. The starting point of the "sub" level patch is the 1° patch of the level 9 grid. The starting numerical space size of the grid extends from 60' to 64'. GeoSOT's 10-level grid is divided into four equal parts of 64'*64'. The size of each level 10 grid is 32'*32'. Grids from level 10 to 15 are "sub" level grids, and the subdivision method is recursive according to the above rules, such as Figure 2As shown. Levels 16 to 21 are "second" level grids. The "second" level grid division method refers to the graded grid, that is, the range of the 15th level 1'*1' face is extended to 64"*64", and the above-mentioned face patches are divided according to the quartering method. The grids of levels 22 to 32 below the second level are strictly divided according to the quartering method, and the size of the 32nd level grid is 1 / 2048"*1 / 2048".

[0053] However, at each level of GeoSOT, the Euclidean distance of the grid is different, and the size of the grid gradually decreases from the equator to the poles. It is more difficult to generalize calculations that are sensitive to grid size. For example, it is difficult to generalize an artificial intelligence (AI) model trained in an area near the equator to Iceland, or the training convergence is difficult or the model suffers from precision loss.

[0054] In order to solve the above problems, the present application embodiment provides a grid coding method, referring to Figure 3 , Figure 3 A schematic diagram of a grid coding method provided in an embodiment of the present application is shown in FIG. Figure 3 As shown, the grid coding method provided in the embodiment of the present application includes:

[0055] 301. Obtain a target projection, where the target projection is a projection of the earth's surface on a plane, and an area difference between two corresponding projection areas of any two areas of the same area on the earth's surface on the target projection is less than a preset value.

[0056] In an embodiment of the present application, a map projection method can be used to convert the earth's surface into a plane to obtain a target projection, wherein the map projection method can be an overall projection or a zone projection, that is, the target projection can be obtained by overall projection or zone projection of the earth's surface.

[0057] In existing implementations, the GeoSOT grid is obtained by directly dividing the grid on a square surface of equal longitude and latitude obtained by deforming the earth's surface, and the Euclidean distances in different grids divided by longitude and latitude are different. In the embodiment of the present application, the target projection is not deformed into a square plane, but a square plane that can cover the target projection is re-determined outside the target projection.

[0058] The plane of the earth's surface of equal longitude and latitude is obtained by reducing and deforming the earth's surface by a certain ratio. On the plane of equal longitude and latitude on the earth's surface, the reduction ratio of the area far from the equator is much smaller than the reduction ratio of the area close to the equator. In other words, the area difference between the two areas corresponding to the area far from the equator and the area close to the equator on the plane of the earth's surface of equal longitude and latitude is very large. In this embodiment, the area difference between the two projection areas corresponding to any two areas of the same area on the earth's surface on the target projection is less than a preset value. In other words, any two areas of the same area on the earth's surface are reduced by a certain ratio to obtain the target projection, and the reduction ratio difference between any two areas of the same area on the earth's surface is very small.

[0059] The target projection can more accurately represent the Euclidean distance of each point on the earth's surface. Specifically, the distance between any two points on the target projection can represent the Euclidean distance between the corresponding two physical locations on the earth's surface. The so-called Euclidean distance refers to the Euclidean metric. In mathematics, the Euclidean distance or Euclidean metric is the "ordinary" (i.e. straight-line) distance between two points in Euclidean space.

[0060] Figure 1 The distance between two points on the plane of equal longitude and latitude on the surface of the earth cannot accurately represent the Euclidean distance of each point on the surface of the earth. Specifically, the plane of equal longitude and latitude on the surface of the earth is obtained by reducing and deforming the surface of the earth by a certain proportion. On the plane of equal longitude and latitude on the surface of the earth, the reduction ratio of the distance between two points far from the equator is much smaller than the reduction ratio of the distance between two points close to the equator. For example, the distance between two points close to the equator on the plane of equal longitude and latitude represents the distance between the corresponding actual physical position points on the surface of the earth, which is longer than the distance between two points with the same distance on the plane of equal longitude and latitude far from the equator.

[0061] In this embodiment, the target projection is an approximate undistorted projection of the earth's surface, and the distance between each point on the target projection is obtained by reducing the earth's surface by a certain ratio, and the reduction ratio of the distance between each point on the earth's surface is basically the same, so the target projection can more accurately represent the Euclidean distance of each point on the earth's surface. Specifically, the first target point and the second target point are any two points on the target projection, the first target point is the projection of the first position point on the earth's surface, the second target point is the projection of the second position point on the earth's surface, the distance between the first target point and the second target point is the first distance, the distance between the first position point and the second position point is the second distance, and the ratio between the first distance and the second distance is within a preset range. The preset range can be selected according to actual conditions, and the preset range can represent the reduction ratio of the earth's surface when performing plane projection. Since the reduction ratio of the distance between any two points on the earth's surface when projecting is within the preset range, the target projection can more accurately represent the Euclidean distance of each point on the earth's surface.

[0062] It should be understood that the present application does not limit the areas of any two regions of the same area on the earth's surface to be completely consistent with each other, but the difference is less than a preset value.

[0063] It should be understood that this embodiment does not limit the reduction ratio of the distance between each point to be completely consistent when the projection of the earth's surface is performed, but is within a preset range.

[0064] Specifically, the target projection is taken as an example in which the surface of the earth is projected using the universal transverse mercator (UTM) projection on a plane.

[0065] In this embodiment, the UTM projection is a transverse angular secant elliptic cylinder projection. The elliptic cylinder secants the earth into two contour circles at 80° south latitude and 84° north latitude. After projection, there is no deformation on the two secant meridians, and the length ratio on the central meridian is 0.9996. The central meridian of the projection is projected as the ordinate axis, and the straight line after the projection of the equatorial latitude is the transverse axis. The 6° zone of the UTM projection divides the world into 60 projection zones, and the longitude difference of each zone is 6°. The first zone is from 180° west longitude to 174° west longitude, and is numbered continuously eastward.

[0066] The relationship between UTM plane coordinates (x, y) and geodetic coordinates (L, B) can be expressed by the following formula:

[0067]

[0068]

[0069] Where: T = tan 2 B; C = e ′2 cos 2 B; A = (L-L0) cosB;

[0070]

[0071]

[0072] Among them: a is the major semi-axis of the earth ellipsoid, b is the minor semi-axis of the earth ellipsoid, e is the first eccentricity, e′ is the second eccentricity, and L0 is the central longitude of the zone.

[0073] like Figure 4 As shown, Figure 4 A schematic diagram of a target projection provided in an embodiment of the present application is shown in FIG. Figure 4 The target projection shown in is obtained by UTM projection of the earth's surface, where the central meridian of the target projection is projected as the ordinate axis and the straight line after the equatorial latitude is projected as the abscissa axis. The UTM projection divides the world into 60 projection zones, each with a longitude span of 6°.

[0074] 302. Determine a square plane according to the target projection, wherein an area where the target projection is located is covered by the square plane.

[0075] In the embodiment of the present application, after obtaining the target projection, a square plane can be determined according to the target projection, and the area where the target projection is located is covered by the square plane. In order to completely cover the area where the target projection is located, the area of ​​the square plane should be larger than the target projection, and the side of the square plane should be outside the area of ​​the target projection.

[0076] It should be understood that the target projection is non-square, so the square plane includes an area overlapping with the target projection and also includes an area not overlapping with the target projection.

[0077] Next, we describe how to determine a square plane that can cover the target projection.

[0078] First, it is necessary to determine the direction of the square plane, which refers to the direction of the horizontal axis and the vertical axis included in the square plane. The square plane may include a horizontal axis and a vertical axis, wherein the horizontal axis may be parallel to one side of the square plane, and the vertical axis is perpendicular to the horizontal axis. More specifically, the horizontal axis may be an axis parallel to one side of the square plane and passing through the center of the square plane, and the vertical axis is an axis perpendicular to the horizontal axis and passing through the center of the square plane.

[0079] In the embodiment of the present application, the earth surface includes an equatorial latitude and a target meridian perpendicular to the equatorial latitude, the target projection includes a projected equatorial latitude and a projected target meridian, the horizontal axis direction of the square plane is consistent with the direction of the projected equatorial latitude, and the vertical axis direction of the square plane is consistent with the direction of the projected target meridian. Figure 5 As shown, Figure 5 A schematic diagram of a square plane provided in an embodiment of the present application is shown in FIG. Figure 5 The horizontal axis of the middle square plane can be axis x(m), and the vertical axis can be y(m), and the horizontal axis direction of the square plane is consistent with the direction of the projected equatorial latitude, and the vertical axis direction of the square plane is consistent with the direction of the projected target longitude.

[0080] Next, the side length of the square plane can be determined. Specifically, in one implementation, the side length of the square plane can be determined according to the target length of the target projection in the direction of the equatorial latitude after the projection, wherein the side length of the square plane is greater than or equal to the target length, and the side length of the square plane is a preset multiple of an integer power of 2, and the square plane is determined according to the side length of the square plane so that the area where the target projection is located is covered by the square plane. In one implementation, the central horizontal axis of the square plane overlaps with the equatorial latitude after the projection.

[0081] The reason why the side length of the square plane is a preset multiple of an integer power of 2 is that when performing subsequent encoding, the square plane needs to be divided into two equal parts in the horizontal and vertical directions. In order to obtain an integer power of 2 square grid after the division, it is necessary to ensure that the side length of the square plane is a preset multiple of an integer power of 2. For example, the side length of the square plane can be Wherein, L is the length of the equatorial latitude line in the target projection, and S is the scaling factor, which is the preset multiple in the above embodiment.

[0082] In an embodiment of the present application, the preset multiple can be obtained, and the side length of the square plane can be determined based on the target length of the target projection in the direction of the equatorial latitude after the projection and the preset multiple. In this embodiment, the size of the square grid obtained after the subdivision can be controlled by controlling the size of the preset multiple.

[0083] In the embodiment of the present application, a plane rectangular coordinate system can be established with the straight line after the projection of the equator latitude in the target projection as the horizontal axis and the projection meridian perpendicular to it as the vertical axis, and the target projection can be divided into four parts, namely the northeast projection facet, the southeast projection facet, the northwest projection facet and the southwest projection facet. Let the length of the northeast projection facet and the northwest projection facet in the horizontal axis direction be the larger one, denoted as L meters, and the northeast projection facet, the southeast projection facet, the northwest projection facet and the southwest projection facet are all expanded to a square plane of size W*W, where, Where s is the scaling factor, thereby obtaining a 2W*2W square plane.

[0084] For example, the target projection is obtained by UTM projection of the earth's surface on a plane. The straight line after the equatorial latitude line projection is taken as the horizontal axis, that is, it coincides with the horizontal axis coordinate of the UTM projection, and the coordinate origin is the prime meridian projection as the vertical axis. A two-dimensional plane rectangular coordinate system is established, and the earth projection plane is divided into four parts, namely the northeast projection facet, the southeast projection facet, the northwest projection facet and the southwest projection facet, such as Figure 5 As shown. Let the scaling factor s = 1, The northeast projection facet, southeast projection facet, northwest projection facet and southwest projection facet are all expanded to a square grid of size W*W.

[0085] 303. Divide the square plane into two equal parts multiple times in the horizontal and vertical directions to obtain multiple square grids.

[0086] In an embodiment of the present application, the square plane may be bisected multiple times in the horizontal and vertical directions to obtain multiple groups of square grids, each group of square grids including multiple square grids, wherein each time the square plane is bisected, a group of square grids may be obtained. Specifically, the square plane may be bisected multiple times in the horizontal and vertical directions to obtain multiple groups of square grids, each group of square grids including multiple square grids, wherein the multiple bisecting includes the Nth bisecting and the N+1th bisecting, and after the Nth bisecting, M square grids are obtained, and after the N+1th bisecting, M*4 square grids are obtained, and the M*4 square grids are obtained by bisecting each of the M square grids in the horizontal and vertical directions.

[0087] In an embodiment of the present application, each time a square plane is bisected in the horizontal and vertical directions, a group of square grids can be obtained. When bisecting the square plane for the next time, each square grid in the group of square grids obtained by the previous bisection can be bisected in the horizontal and vertical directions to obtain a new group of square grids.

[0088] For details, please refer to Figure 6 , after two bisections, we can get Figure 6 The 16 square grids shown in the figure can be used for the next mesh generation process. Figure 6 Each of the 16 square grids shown is bisected in the horizontal and vertical directions to obtain 64 square grids (specifically, as shown in FIG. Figure 7 shown).

[0089] 304. Encode the multiple square grids to obtain multiple codes, each code corresponding to a square grid, and each code is used to indicate an area covered by the corresponding square grid.

[0090] In an embodiment of the present application, after the square plane is bisected multiple times in the horizontal and vertical directions to obtain multiple square grids, the multiple square grids can be encoded to obtain multiple codes, each code corresponding to a square grid, and each code is used to indicate the area covered by the corresponding square grid.

[0091] In an embodiment of the present application, the multiple square grids included in each group of square grids in the multiple groups of square grids may be encoded to obtain multiple groups of codes, each group of codes including multiple codes.

[0092] It should be understood that the timing between step 303 and step 304 is not limited in the embodiment of the present application. Specifically, after the square plane is bisected once in the horizontal and vertical directions, a group of square grids can be obtained, and then the group of square grids can be encoded to obtain the encoding of the group of square grids, and then the next bisection can be performed, that is, each square grid in the group of square grids obtained by the previous division can be bisected in the horizontal and vertical directions to obtain a new group of square grids, and then the group of square grids can be encoded to obtain the encoding of the group of square grids, and so on.

[0093] It should be understood that each code obtained after encoding the square grid is used to indicate the area covered by the corresponding square grid, and the present application does not limit the specific encoding method.

[0094] Exemplarily, the target projection is obtained by UTM projection of the earth's surface on a plane. For the expanded 2W*2W plane, it is divided into two equal parts in the horizontal and vertical directions to form four square grids of equal size, and each square grid generated by the division is encoded according to certain rules. After that, each new square grid generated by the above division is further divided into two equal parts in the horizontal and vertical directions to form four square grids of equal size, and each square grid generated by the division is encoded according to certain rules. This cycle is repeated until the encoding of each square grid generated by the division meets the requirements, such as the side length of the smallest square grid is 0.5 meters.

[0095] In the embodiment of the present application, a square grid is obtained by dividing a square plane, which has the characteristics of a global discrete grid such as global coverage, uniqueness, hierarchy, hierarchical membership and coding operation, and the grid size according to grid sampling or statistical data is a square, and the size of the square can be selected at an appropriate level according to application needs, or a suitable scaling factor s can be selected when constructing the divided grid, so as to meet the application requirements of algorithms that are sensitive to grid size.

[0096] In this embodiment, the grid is divided on the square plane constructed on the target projection, so that the size of the earth surface corresponding to the square grid obtained by the same level of division is basically the same, which is conducive to the calculation sensitive to the size of the grid data, such as the AI ​​model with convolutional neural networks (CNN) as the feature extraction layer. At the same time, the conversion cost of the raster data in the local coordinate system aligned with the local area of ​​the equidistant global discrete grid is low.

[0097] It should be understood that, based on the above-mentioned plane rectangular coordinates, the height dimension can be added to construct a three-dimensional rectangular coordinate, and the height dimension can be equally divided, which can be extended to a 3D equidistant global discrete grid.

[0098] The embodiment of the present application provides a grid encoding method, the method comprising: obtaining a target projection, the target projection being a projection of the earth's surface on a plane, and the area difference between the two projection areas corresponding to any two areas of the same area on the earth's surface on the target projection is less than a preset value; determining a square plane according to the target projection, the area where the target projection is located being covered by the square plane; bisecting the square plane multiple times in the horizontal and vertical directions to obtain a plurality of square grids; encoding the multiple square grids to obtain a plurality of codes, each code corresponding to a square grid, and each code being used to indicate the area covered by the corresponding square grid. In the above manner, grid division is performed on the square plane constructed on the target projection, so that the division sizes of the earth's surface corresponding to the square grids obtained by the same level of division are substantially consistent, which is beneficial to calculations that are sensitive to the size of the grid data.

[0099] Reference Figure 8 , Figure 8 A schematic diagram of the structure of a grid coding device provided in an embodiment of the present application is shown in FIG. Figure 8 As shown, the grid coding device 800 provided in the embodiment of the present application includes:

[0100] The acquisition module 801 is used to acquire a target projection, where the target projection is a projection of the earth's surface on a plane.

[0101] The specific description of the acquisition module 801 can refer to the description of step 301 and the corresponding embodiment, which will not be repeated here.

[0102] A determination module 802 is used to determine a square plane according to the target projection, and the area where the target projection is located is covered by the square plane;

[0103] The specific description of the determination module 802 can refer to the description of step 302 and the corresponding embodiment, which will not be repeated here.

[0104] A grid division module 803 is used to divide the square plane into two equal parts multiple times in the horizontal axis direction and the vertical axis direction to obtain multiple square grids;

[0105] The specific description of the grid generation module 803 can refer to the description of step 303 and the corresponding embodiment, which will not be repeated here.

[0106] The encoding module 804 is used to encode the multiple square grids to obtain multiple codes, each code corresponds to a square grid, and each code is used to indicate the area covered by the corresponding square grid.

[0107] The specific description of the encoding module 804 can refer to the description of step 304 and the corresponding embodiment, which will not be repeated here.

[0108] In a possible implementation, the first target point and the second target point are any two points on the target projection, the first target point is the projection of a first position point on the earth's surface, the second target point is the projection of a second position point on the earth's surface, the distance between the first target point and the second target point is a first distance, the distance between the first position point and the second position point is a second distance, and the ratio between the first distance and the second distance is within a preset range.

[0109] In one possible implementation, the earth's surface includes an equatorial latitude and a target longitude perpendicular to the equatorial latitude, the target projection includes a projected equatorial latitude and a projected target longitude, the horizontal axis direction of the square plane is consistent with the direction of the projected equatorial latitude, and the vertical axis direction of the square plane is consistent with the direction of the projected target longitude.

[0110] In a possible implementation, the determination module is used to determine the side length of the square plane according to the target length of the target projected in the direction of the equatorial latitude after the projection, wherein the side length of the square plane is greater than or equal to the target length, and the side length of the square plane is a preset multiple of an integer power of 2;

[0111] The square plane is determined according to the side length of the square plane, so that the area where the target projection is located is covered by the square plane.

[0112] In a possible implementation, the determination module is used to obtain the preset multiple, and determine the side length of the square plane according to the target length of the target projected in the direction of the equatorial latitude after the projection and the preset multiple.

[0113] In a possible implementation, the central horizontal axis of the square plane overlaps with the projected equatorial latitude.

[0114] In a possible implementation, the grid division module is used to perform multiple bisection divisions on the square plane in the horizontal and vertical directions to obtain multiple groups of square grids, each group of square grids includes multiple square grids, wherein the multiple bisection divisions include the Nth bisection division and the N+1th bisection division, and after the Nth bisection division, M square grids are obtained, and after the N+1th bisection division, M*4 square grids are obtained, and the M*4 square grids are obtained by bisectioning each square grid in the M square grids in the horizontal and vertical directions;

[0115] The encoding module is used to encode the multiple square grids included in each group of square grids in the multiple groups of square grids to obtain multiple groups of codes, each group of codes including multiple codes.

[0116] In a possible implementation, the target projection is obtained by performing an overall projection or a zone projection on the earth's surface.

[0117] The embodiment of the present application provides a grid coding device, including: an acquisition module 801, used to acquire a target projection, the target projection being a projection of the earth's surface on a plane; a determination module 802, used to determine a square plane according to the target projection, the area where the target projection is located being covered by the square plane; a grid division module 803, used to divide the square plane into two equal parts multiple times in the horizontal and vertical directions to obtain a plurality of square grids; an encoding module 804, used to encode the plurality of square grids to obtain a plurality of codes, each code corresponding to a square grid, and each code is used to indicate the area covered by the corresponding square grid. Grid division is performed on the square plane constructed on the target projection, so that the size of the earth's surface corresponding to the square grids obtained by the same level of division is basically the same, which is conducive to calculations sensitive to the size of grid data.

[0118] The present application also provides a non-volatile computer-readable storage medium, which contains computer instructions. When the computer instructions are executed by a computer, the grid coding method in the above embodiment can be implemented.

[0119] Please refer to Fig. 9 , is a schematic diagram of the structure of a computer system provided in this embodiment. The computer system may be a terminal device (or referred to as an intelligent terminal) or a server. As shown in the figure, the computer system includes a communication module 810, a sensor 820, a user input module 830, an output module 840, a processor 850, an audio and video input module 860, a memory 870, and a power supply 880. Furthermore, the computer system provided in this embodiment may also include an NPU 890.

[0120] The communication module 810 may include at least one module that enables the computer system to communicate with a communication system or other computer systems. For example, the communication module 810 may include one or more of a wired network interface, a broadcast receiving module, a mobile communication module, a wireless Internet module, a local area communication module, and a location (or positioning) information module. These various modules are all implemented in various ways in the prior art, and this application will not describe them one by one.

[0121] The sensor 820 may sense the current state of the system, such as an open / closed state, a position, whether there is contact with a user, a direction, and acceleration / deceleration, and the sensor 820 may generate a sensing signal for controlling the operation of the system.

[0122] The user input module 830 is used to receive input digital information, character information or contact touch operation / contactless gesture, and receive signal input related to user settings and function control of the system, etc. The user input module 830 includes a touch panel and / or other input devices.

[0123] The output module 840 includes a display panel for displaying information input by the user, information provided to the user, or various menu interfaces of the system. Optionally, the display panel can be configured in the form of a liquid crystal display (LCD) or an organic light-emitting diode (OLED). In some other embodiments, the touch panel can cover the display panel to form a touch display screen. In addition, the output module 840 can also include an audio output module, an alarm, and a tactile module.

[0124] The audio and video input module 860 is used to input audio signals or video signals. The audio and video input module 860 may include a camera and a microphone.

[0125] The power supply 880 may receive external power and internal power under the control of the processor 850 and provide power required for the operation of various components of the system.

[0126] The processor 850 includes one or more processors. For example, the processor 850 may include a central processing unit and a graphics processing unit. The central processing unit in this application has multiple cores and belongs to a multi-core processor. These multiple cores can be integrated on the same chip or each can be an independent chip.

[0127] The memory 870 stores computer programs, which include an operating system program 872 and an application program 871. Typical operating systems include Microsoft's Windows, Apple's MacOS, and other systems for desktops or notebooks, and systems for mobile terminals such as the Android system developed by Google. The method provided in the above embodiment can be implemented by software and can be considered as a specific implementation of the operating system program 872. The memory 870 can be one or more of the following types: flash memory, hard disk type memory, micro multimedia card type memory, card memory (such as SD or XD memory), random access memory (random access memory, RAM), static random access memory (static RAM, SRAM), read only memory (read only memory, ROM), electrically erasable programmable read-only memory (electrically erasable programmable read-only memory, EEPROM), programmable read-only memory (programmable ROM, PROM), rollback protection storage block (replay protected memory block, RPMB), magnetic storage, disk or optical disk. In some other embodiments, the memory 870 may also be a network storage device on the Internet, and the system may perform operations such as updating or reading on the memory 870 on the Internet.

[0128] The processor 850 is used to read the computer program in the memory 870 and then execute the method defined by the computer program. For example, the processor 850 reads the operating system program 872 to run the operating system and implement various functions of the operating system on the system, or reads one or more application programs 871 to run applications on the system.

[0129] The memory 870 also stores other data 873 in addition to the computer program.

[0130] NPU 890 is mounted on the main processor 850 as a coprocessor to execute the tasks assigned to it by the main processor 850. In this embodiment, NPU 890 can be called by one or more sub-threads of face recognition TA to implement some complex algorithms involved in face recognition. Specifically, the sub-threads of face recognition TA run on multiple cores of the main processor 850, and then the main processor 850 calls NPU 890, and the results implemented by NPU 890 are returned to the main processor 850.

[0131] The connection relationship between the above modules is only an example. The grid coding method provided in any embodiment of the present application can also be applied to terminal devices or servers with other connection methods, for example, all modules are connected through a bus.

[0132] Fig.10 9 is a schematic diagram of the structure of an NPU 900 provided in this embodiment. The NPU 900 is connected to a main processor and an external memory. The core part of the NPU 900 is an operation circuit 903, which is controlled by a controller 904 to extract data from the memory and perform mathematical operations.

[0133] In some implementations, the operation circuit 903 includes multiple processing engines (PEs) inside. In some implementations, the operation circuit 903 is a two-dimensional systolic array. The operation circuit 903 can also be a one-dimensional systolic array or other electronic circuits capable of performing mathematical operations such as multiplication and addition. In other implementations, the operation circuit 903 is a general-purpose matrix processor.

[0134] For example, assume there is an input matrix A, a weight matrix B, and an output matrix C. The operation circuit 903 takes the corresponding data of the matrix B from the weight memory 902 and caches it on each PE of the operation circuit 903. The operation circuit 903 takes the matrix A data from the input memory 901 and performs a matrix operation with the matrix B, and the partial result or the final result of the matrix is ​​stored in the accumulator 908.

[0135] The unified memory 906 is used to store input data and output data. The weight data is directly transferred to the weight memory 902 through the storage unit access controller 905 (eg, direct memory access controller, DMAC).

[0136] The input data is also transferred to the unified memory 906 through the storage unit access controller 905 .

[0137] The bus interface unit 910 (BIU) is used for the interaction between the AXI (advanced extensible interface) bus and the storage unit access controller 905 and the instruction fetch buffer 909 (instruction fetch buffer).

[0138] The bus interface unit 910 is used for the instruction fetch memory 909 to obtain instructions from the external memory, and is also used for the storage unit access controller 905 to obtain the original data of the input matrix A or the weight matrix B from the external memory.

[0139] The storage unit access controller 905 is mainly used to transfer input data in the external memory to the unified memory 906 or to transfer weight data to the weight memory 902 or to transfer input data to the input memory 901 .

[0140] The vector calculation unit 907 generally includes a plurality of operation processing units, and performs further processing on the output of the operation circuit 903 when necessary, such as vector multiplication, vector addition, exponential operation, logarithmic operation, and / or size comparison, etc.

[0141] In some implementations, the vector calculation unit 907 can store the processed vector in the unified memory 906. For example, the vector calculation unit 907 can apply a nonlinear function to the output of the operation circuit 903, such as a vector of accumulated values, to generate an activation value. In some implementations, the vector calculation unit 907 generates a normalized value, a merged value, or both. In some implementations, the processed vector can be used as an activation input of the operation circuit 903.

[0142] The instruction fetch memory 909 connected to the controller 904 is used to store instructions used by the controller 904 .

[0143] The unified memory 906, the input memory 901, the weight memory 902 and the instruction fetch memory 909 are all on-chip memories. The external memory in the figure is independent of the NPU hardware architecture.

[0144] It should be noted that the method for configuring the address translation relationship provided in this embodiment can also be applied to non-terminal computer devices, such as cloud servers.

[0145] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0146] In the several embodiments provided in the present application, it should be understood that the disclosed systems, devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of units is only a logical function division. There may be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be an indirect coupling or communication connection through some interfaces, devices or units, which can be electrical, mechanical or other forms.

[0147] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0148] In addition, each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0149] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present application is essentially or the part that contributes to the prior art or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions for a computer device (which can be a personal computer, server, or other network device, etc.) to execute the present application. Figure 3 All or part of the steps of the method described in the embodiment. The aforementioned storage medium includes: a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk, and other media that can store program codes.

[0150] As described above, the above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them. Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present application.

Claims

1. A grid coding method, characterized in that: The method comprises: Acquire a target projection, where the target projection is a projection of the earth's surface on a plane, and the area difference between two corresponding projection areas of any two areas of the same area on the earth's surface on the target projection is less than a preset value; Determine a square plane according to the target projection, wherein the area where the target projection is located is covered by the square plane; Dividing the square plane into two equal parts multiple times in the horizontal and vertical directions to obtain multiple square grids; The multiple square grids are encoded to obtain multiple codes, each code corresponds to a square grid, and each code is used to indicate an area covered by the corresponding square grid.

2. The method according to claim 1, characterized in that The first target point and the second target point are any two points on the target projection, the first target point is the projection of the first position point on the earth's surface, the second target point is the projection of the second position point on the earth's surface, the distance between the first target point and the second target point is the first distance, the distance between the first position point and the second position point is the second distance, and the ratio between the first distance and the second distance is within a preset range.

3. The method according to claim 1 or 2, characterized in that: The earth's surface includes an equatorial latitude and a target longitude perpendicular to the equatorial latitude, the target projection includes a projected equatorial latitude and a projected target longitude, the horizontal axis direction of the square plane is consistent with the direction of the projected equatorial latitude, and the vertical axis direction of the square plane is consistent with the direction of the projected target longitude.

4. The method according to claim 3, characterized in that The step of determining a square plane according to the target projection comprises: Determine the side length of the square plane according to the target length of the target projected in the direction of the equatorial latitude after the projection, wherein the side length of the square plane is greater than or equal to the target length, and the side length of the square plane is a preset multiple of an integer power of 2; The square plane is determined according to the side length of the square plane, so that the area where the target projection is located is covered by the square plane.

5. The method according to claim 4, characterized in that The determining the side length of the square plane according to the target length of the target projected in the direction where the equatorial latitude line is located after the projection includes: The preset multiple is obtained, and the side length of the square plane is determined according to the target length of the target projected in the direction of the equatorial latitude after the projection and the preset multiple.

6. The method according to claim 3, characterized in that: The central horizontal axis of the square plane overlaps with the projected equatorial latitude.

7. The method according to claim 1 or 2, characterized in that: The step of bisecting the square plane multiple times in the horizontal and vertical directions to obtain multiple square grids includes: The square plane is bisected multiple times in the horizontal and vertical directions to obtain multiple groups of square grids, each group of square grids includes multiple square grids, wherein the multiple bisecting includes the Nth bisecting and the N+1th bisecting, and after the Nth bisecting, M square grids are obtained, and after the N+1th bisecting, M*4 square grids are obtained, and the M*4 square grids are obtained by bisecting each square grid in the M square grids in the horizontal and vertical directions; The encoding of the plurality of square grids to obtain a plurality of codes comprises: The multiple square grids included in each group of square grids in the multiple groups of square grids are encoded to obtain multiple groups of codes, each group of codes including multiple codes.

8. The method according to claim 1 or 2, characterized in that: The target projection is obtained by performing overall projection or zone projection on the earth surface.

9. A grid coding device, characterized in that: The device comprises: An acquisition module, used for acquiring a target projection, where the target projection is a projection of the earth's surface on a plane; A determination module, used for determining a square plane according to the target projection, wherein the area where the target projection is located is covered by the square plane; A grid division module, used for dividing the square plane into two equal parts multiple times in the horizontal axis direction and the vertical axis direction to obtain multiple square grids; The encoding module is used to encode the multiple square grids to obtain multiple codes, each code corresponds to a square grid, and each code is used to indicate the area covered by the corresponding square grid.

10. The device according to claim 9, characterized in that The first target point and the second target point are any two points on the target projection, the first target point is the projection of the first position point on the earth's surface, the second target point is the projection of the second position point on the earth's surface, the distance between the first target point and the second target point is the first distance, the distance between the first position point and the second position point is the second distance, and the ratio between the first distance and the second distance is within a preset range.

11. The device according to claim 9 or 10, characterized in that The earth's surface includes an equatorial latitude and a target longitude perpendicular to the equatorial latitude, the target projection includes a projected equatorial latitude and a projected target longitude, the horizontal axis direction of the square plane is consistent with the direction of the projected equatorial latitude, and the vertical axis direction of the square plane is consistent with the direction of the projected target longitude.

12. The device according to claim 11, characterized in that The determination module is used to determine the side length of the square plane according to the target length of the target projected in the direction of the equatorial latitude after the projection, wherein the side length of the square plane is greater than or equal to the target length, and the side length of the square plane is a preset multiple of an integer power of 2; The square plane is determined according to the side length of the square plane, so that the area where the target projection is located is covered by the square plane.

13. The device according to claim 12, characterized in that The determination module is used to obtain the preset multiple and determine the side length of the square plane according to the target length of the target projected in the direction of the equatorial latitude after the projection and the preset multiple.

14. The device according to any one of claims 12 to 13, characterized in that: The central horizontal axis of the square plane overlaps with the projected equatorial latitude.

15. The device according to claim 9 or 10, characterized in that The grid division module is used to divide the square plane into two equal parts multiple times in the horizontal and vertical directions to obtain multiple groups of square grids, each group of square grids includes multiple square grids, wherein the multiple bisection divisions include the Nth bisection division and the N+1th bisection division, and after the Nth bisection division, M square grids are obtained, and after the N+1th bisection division, M*4 square grids are obtained, and the M*4 square grids are obtained by bisectioning each square grid in the M square grids in the horizontal and vertical directions; The encoding module is used to encode the multiple square grids included in each group of square grids in the multiple groups of square grids to obtain multiple groups of codes, each group of codes including multiple codes.

16. The device according to claim 9 or 10, characterized in that The target projection is obtained by performing overall projection or zone projection on the earth surface.

17. A computer system, characterized in that: comprising a memory and a processor, wherein: The memory is used to store computer-readable instructions; the processor is used to read the computer-readable instructions and implement the method according to any one of claims 1-8.

18. A computer storage medium, characterized in that: Computer-readable instructions are stored, and when the computer-readable instructions are executed by a processor, the method according to any one of claims 1 to 8 is implemented.

19. A computer program product, characterized in that The method comprises codes, which, when executed, are used to implement the method according to any one of claims 1 to 8.

Images