CN108601031B - Allocation mapping method and device for communication cells and geographic areas - Google Patents

Allocation mapping method and device for communication cells and geographic areas Download PDF

Info

Publication number
CN108601031B
CN108601031B CN201810273229.0A CN201810273229A CN108601031B CN 108601031 B CN108601031 B CN 108601031B CN 201810273229 A CN201810273229 A CN 201810273229A CN 108601031 B CN108601031 B CN 108601031B
Authority
CN
China
Prior art keywords
cell
area
sub
thiessen
latitude
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810273229.0A
Other languages
Chinese (zh)
Other versions
CN108601031A (en
Inventor
蔡子龙
安瑞虹
于飞
傅俊峰
王题
王一骢
韦广林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China United Network Communications Group Co Ltd
China Information Technology Designing and Consulting Institute Co Ltd
Original Assignee
China United Network Communications Group Co Ltd
China Information Technology Designing and Consulting Institute Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China United Network Communications Group Co Ltd, China Information Technology Designing and Consulting Institute Co Ltd filed Critical China United Network Communications Group Co Ltd
Priority to CN201810273229.0A priority Critical patent/CN108601031B/en
Publication of CN108601031A publication Critical patent/CN108601031A/en
Application granted granted Critical
Publication of CN108601031B publication Critical patent/CN108601031B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W16/00Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
    • H04W16/18Network planning tools

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

The embodiment of the invention provides a method and a device for apportionment mapping of a communication cell and a geographic area, relates to the field of communication geographic data processing, and can more accurately obtain the matching relationship between the communication cell and the geographic area. The method comprises the following steps: carrying out bias processing on the longitude and latitude of the base station according to the geographic data of the communication cell to obtain the longitude and latitude of the communication cell; building a Thiessen cell corresponding to the communication cell according to the longitude and latitude of the communication cell; dividing the grid into at least one sub-area according to boundary data in the geographic data of the Thiessen cell and the geographic data of the grid; calculating the area of the sub-region and the longitude and latitude of the centroid of the sub-region according to the geographic data of the sub-region; matching the sub-area with the Thiessen cell according to the longitude and latitude of the centroid of the sub-area and the geographic data of the Thiessen cell; calculating the area of the Thiessen cell according to the geographic data of the Thiessen cell; calculating the ratio of the area of each subregion matched with the Thiessen cell to the area of the Thiessen cell.

Description

Allocation mapping method and device for communication cells and geographic areas
Technical Field
The invention relates to the field of communication geographic data processing, in particular to a method and a device for apportionment mapping of communication cells and geographic areas.
Background
Currently, network planning of a single city of an operator mainly adopts a form of dividing a geographic area (a grid and a micro grid) for planning. The geographical area planning is mainly to closely relate the service requirement, the current network situation and the planning scheme through the geographical position, so as to realize the accuracy of the network planning construction. Each micro-grid or mesh has a corresponding dynamically expanded geographic information base, the information base content comprising: basic information of the micro-grid or grid, internal user structure, service distribution, network structure, etc. In order to maximize the investment benefit, operators propose investment planning guidance and introduce charging system data as a consideration factor. However, since the user data in the charging system does not have geographical location information, investment planning cannot be accurately performed according to the user traffic of different geographical areas (grids or microgrids), so that the communication cell generated by the service needs to be indirectly associated with the geographical areas, and then the approximate traffic of different geographical areas needs to be obtained according to the association relationship, so that the mapping relationship between the communication cell and the geographical areas (grids or microgrids) needs to be known first in the analysis process.
The mapping method of the communication cell and the geographic area which is commonly used at present mainly judges the geographic area according to the longitude and latitude of the communication cell, taking a grid as an example, the grid to which the longitude and latitude of the communication cell belongs is the corresponding grid; however, in such a case, a small part of the area where the longitude and latitude corresponding position of the communication cell is located is in a certain grid, and the final matching result is that the cell corresponds to the grid, but the grid where the large part of the cell is located is not the grid obtained by matching, which causes inaccuracy of matching between the communication cell and the grid, and accurate traffic in the grid cannot be obtained by using the matching result.
Disclosure of Invention
The embodiment of the invention provides a method and a device for apportionment mapping of a communication cell and a geographic area, which can more accurately obtain the matching relationship between the communication cell and the geographic area.
In order to achieve the above purpose, the embodiment of the invention adopts the following technical scheme:
in a first aspect, a method for apportionment mapping of communication cells and geographic areas is provided, which includes:
acquiring geographic data of at least one communication cell, and performing bias processing on the longitude and latitude of a base station according to the longitude and latitude of the base station in the geographic data of the communication cell and the azimuth angle of the communication cell to acquire the longitude and latitude of the communication cell;
building a Thiessen cell corresponding to the communication cell based on the Thiessen polygon principle according to the longitude and latitude of the communication cell;
acquiring geographic data of at least one grid in a geographic area, and dividing the grid into at least one sub-area according to boundary data in the geographic data of the Thiessen cell and the geographic data of the grid;
calculating the area of the sub-region and the longitude and latitude of the centroid of the sub-region according to the geographic data of the sub-region;
matching the sub-area with the Thiessen cell according to the longitude and latitude of the centroid of the sub-area and the geographic data of the Thiessen cell;
calculating the area of the Thiessen cell according to the geographic data of the Thiessen cell;
calculating the ratio of the area of each subregion matched with the Thiessen cell to the area of the Thiessen cell.
Optionally, before matching the sub-area with the thieson cell according to the longitude and latitude of the centroid of the sub-area and the geographic data of the thieson cell, the method further includes:
judging whether a first target subregion exists in the subregions or not; the centroid of the first target sub-region lies outside or on the boundary of the first target sub-region;
and when the first target sub-area exists in the sub-area, modifying the longitude and latitude of the mass center of the first target sub-area.
Optionally, before calculating the area ratio of each sub-area matched with the tesson cell to the tesson cell, the method further includes:
calculating the ratio of the sum of the areas of all the sub-areas successfully matched with the Thiessen cell to the area of the Thiessen cell;
when the ratio of the sum of the areas of all sub-areas successfully matched with the Thiessen cell to the area of the Thiessen cell is determined to be larger than a preset percentage, determining that all the sub-areas successfully matched with the Thiessen cell are second target sub-areas, and correcting the longitude and latitude of the mass center of the second target sub-areas;
matching the sub-area with the Thiessen cell according to the longitude and latitude of the centroid of the sub-area and the geographic data of the Thiessen cell;
and recalculating the ratio of the sum of the areas of all the sub-areas successfully matched with the Thiessen cell to the area of the Thiessen cell.
Optionally, when it is determined that the ratio of the sum of the areas of all sub-areas successfully matched with the thieson cell to the area of the thieson cell is not greater than the preset percentage, calculating the ratio of the area of each sub-area matched with the thieson cell to the area of the thieson cell.
Optionally, before calculating the area ratio of each sub-area matched with the tesson cell to the tesson cell, the method further includes:
calculating a target ratio of the area of each sub-region matched with the Thiessen cell to the sum of the areas of all the sub-regions matched with the Thiessen cell;
and modifying the area of each sub-area matched with the Thiessen cell according to a preset formula according to the target ratio and the area of the Thiessen cell.
Optionally, when the first target sub-region is a triangle, modifying the longitude and latitude of the centroid of the first target sub-region includes:
and solving the longitude and latitude of any middle point of the first target sub-region according to the geographic data of the first target sub-region, and taking the longitude and latitude of any middle point of the first target sub-region as the longitude and latitude of the corrected centroid of the first target sub-region.
Optionally, when the first target sub-region is a non-triangle, modifying the longitude and latitude of the centroid of the first target sub-region includes:
sequentially selecting vertexes from all vertexes of the first target sub-area until the longitude and latitude of the midpoint of the target are positioned in the first target sub-area; the target midpoint is the midpoint of a connecting line of intersection points of a circle drawn by taking the selected vertex as the center of a circle and two adjacent edges of the selected vertex, wherein the preset length is the radius;
and taking the longitude and latitude of the midpoint of the target as the longitude and latitude of the centroid corrected by the first target subregion.
Optionally, when the second target sub-region is a triangle, modifying the longitude and latitude of the centroid of the second target sub-region includes:
and solving the longitude and latitude of any middle point of the second target sub-region according to the geographic data of the second target sub-region, and taking the longitude and latitude of the middle point of the second target sub-region as the longitude and latitude of the corrected centroid of the second target sub-region.
Optionally, when the second target sub-region is a non-triangular shape, modifying the longitude and latitude of the centroid of the second target sub-region includes:
sequentially selecting vertexes from all vertexes of the second target subregion, and stopping until the longitude and latitude of the midpoint of the target are positioned in the second target subregion; the target midpoint is the midpoint of a connecting line of intersection points of a circle drawn by taking the selected vertex as the center of a circle and two adjacent edges of the selected vertex, wherein the preset length is the radius;
and taking the longitude and latitude of the midpoint of the target as the longitude and latitude of the centroid corrected by the second target subarea.
In a second aspect, an apparatus for apportionment mapping of communication cells and geographic areas is provided, comprising: the system comprises an acquisition module, a Taisen cell module, a segmentation module, a calculation module and a matching module;
the acquisition module is used for acquiring the geographic data of at least one communication cell and carrying out bias processing on the longitude and latitude of the base station according to the longitude and latitude of the base station in the geographic data of the communication cell and the azimuth angle of the communication cell so as to acquire the longitude and latitude of the communication cell;
the Thiessen cell module is used for constructing a Thiessen cell corresponding to the communication cell based on the Thiessen polygon principle according to the longitude and latitude of the communication cell acquired by the acquisition module;
the obtaining module is further used for obtaining geographic data of at least one grid in the geographic area, and the dividing module is used for dividing the grid into at least one sub-area according to boundary data in the geographic data of the Thiessen cell constructed by the Thiessen cell module and the geographic data of the grid obtained by the obtaining module;
the calculation module is used for calculating the area of the subarea and the longitude and latitude of the centroid of the subarea according to the geographical data of the subarea divided by the division module;
the matching module is used for matching the sub-region with the Thiessen cell according to the longitude and latitude of the centroid of the sub-region calculated by the calculating module and the geographic data of the Thiessen cell constructed by the Thiessen cell module;
the calculation module is also used for calculating the area of the Thiessen cell according to the geographic data of the Thiessen cell constructed by the Thiessen cell module;
the calculation module is further used for calculating the ratio of the area of each sub-area matched for the Thiessen cell by the matching module to the area of the Thiessen cell.
Optionally, the system further comprises a judging module and a correcting module;
the judging module is used for judging whether a first target subregion exists in the subregions divided by the dividing module before the matching module matches the subregions with the Thiessen cells according to the longitude and latitude of the centroid of the subregions calculated by the calculating module and the geographic data of the Thiessen cells constructed by the Thiessen cell module; the centroid of the first target sub-region lies outside or on the boundary of the first target sub-region;
when the judgment module determines that the first target sub-region exists in the sub-region, the correction module is used for correcting the longitude and latitude of the mass center of the first target sub-region.
Optionally, the calculating module is further configured to calculate a ratio of a sum of areas of all sub-regions successfully matched with the thieson cell to the area of the thieson cell before calculating the ratio of the area of each sub-region matched with the thieson cell by the matching module to the area of the thieson cell;
when the judging module determines that the ratio of the sum of the areas of all sub-regions successfully matched with the Thiessen cell to the area of the Thiessen cell is greater than the preset percentage, determining that all the sub-regions successfully matched with the Thiessen cell are second target sub-regions, and the correcting module is also used for correcting the second target sub-regions;
the matching module is further used for matching the sub-area with the Thiessen cell again according to the longitude and latitude of the centroid of the sub-area and the geographic data of the Thiessen cell;
the calculation module is further used for recalculating the ratio of the sum of the areas of all the sub-areas successfully matched with the Thiessen cell to the area of the Thiessen cell.
Optionally, the calculating module is further configured to calculate an area ratio of the area of each sub-region matched with the thieson cell to the thieson cell when the determining module determines that the area ratio of the sum of the areas of all sub-regions successfully matched with the thieson cell to the area of the thieson cell is not greater than a preset percentage.
Optionally, the system further comprises an apportioning module;
the calculation module is further used for calculating a target occupation ratio of the area of each sub-region matched with the Thiessen cell to the sum of the areas of all sub-regions matched with the Thiessen cell before calculating the area occupation ratio of each sub-region matched with the Thiessen cell to the Thiessen cell;
and the allocation module is used for modifying the area of each sub-region matched with the Thiessen cell according to a preset formula according to the target ratio calculated by the calculation module and the area of the Thiessen cell.
Optionally, when the first target sub-region is a triangle, the modification module is specifically configured to:
and solving the longitude and latitude of any central line midpoint of the first target sub-region according to the geographic data of the first target sub-region, and taking the longitude and latitude of any central line midpoint of the sub-region as the longitude and latitude of the corrected centroid of the sub-region.
Optionally, when the first target sub-region is a non-triangle, the modification module is specifically configured to:
sequentially selecting vertexes from all vertexes of the first target subarea, and stopping until the longitude and latitude of the midpoint of the target are positioned in the subarea; the target midpoint is the midpoint of a connecting line of intersection points of a circle drawn by taking the selected vertex as the center of a circle and two adjacent edges of the selected vertex, wherein the preset length is the radius;
and taking the longitude and latitude of the midpoint of the target as the longitude and latitude of the centroid corrected by the first target subregion.
Optionally, when the second target sub-region is a triangle, the modification module is specifically configured to:
and solving the longitude and latitude of any central line midpoint of the second target sub-region according to the geographic data of the second target sub-region, and taking the longitude and latitude of any central line midpoint of the sub-region as the longitude and latitude of the corrected centroid of the sub-region.
Optionally, when the second target sub-region is a non-triangle, the modification module is specifically configured to:
sequentially selecting vertexes from all vertexes of the second target subarea, and stopping until the longitude and latitude of the midpoint of the target are positioned in the subarea; the target midpoint is the midpoint of a connecting line of intersection points of a circle drawn by taking the selected vertex as the center of a circle and two adjacent edges of the selected vertex, wherein the preset length is the radius;
and taking the longitude and latitude of the midpoint of the target as the longitude and latitude of the centroid corrected by the second target subarea.
The embodiment of the invention provides a method and a device for apportionment mapping of communication cells and geographic areas, wherein the method comprises the following steps: acquiring geographic data of at least one communication cell, and performing bias processing on the longitude and latitude of a base station according to the longitude and latitude of the base station in the geographic data of the communication cell and the azimuth angle of the communication cell to acquire the longitude and latitude of the communication cell; building a Thiessen cell corresponding to the communication cell based on the Thiessen polygon principle according to the longitude and latitude of the communication cell; acquiring geographic data of at least one grid in a geographic area, and dividing the grid into at least one sub-area according to boundary data in the geographic data of the Thiessen cell and the geographic data of the grid; calculating the area of the sub-region and the longitude and latitude of the centroid of the sub-region according to the geographic data of the sub-region; matching the sub-area with the Thiessen cell according to the longitude and latitude of the centroid of the sub-area and the geographic data of the Thiessen cell; calculating the area of the Thiessen cell according to the geographic data of the Thiessen cell; calculating the ratio of the area of each subregion matched with the Thiessen cell to the area of the Thiessen cell. When the communication cell and the geographic area are matched, the longitude and the latitude of the communication cell are obtained according to the longitude and the latitude and the azimuth of the base station, the Thiessen cell layer of the communication cell is obtained according to the longitude and the latitude of the communication cell, the Thiessen cell layer and the geographic area layer are overlapped and cut, so that the geographic area layer is divided into a plurality of sub-areas by the boundary in the Thiessen cell layer, the area occupation ratio values of the sub-areas in different Thiessen cell areas which are correspondingly matched with each other are calculated, the area matching relation with more than one area between the Thiessen cell layer and the geographic area layer can be obtained, and then an operator can obtain the service volume of different geographic areas according to the area matching relation with more than one area between the Thiessen cell layer and the geographic area layer, namely the sharing mapping result. In the technical scheme provided by the implementation of the invention, the geographical area map layer is divided by using the area boundary in the Thiessen cell map layer of the communication cell, and the ratio of each divided sub-area to the area of the respectively matched Thiessen cell is obtained, so that the matching relationship between the Thiessen cell and the geographical area grid of each communication cell is accurately obtained from the level of the area ratio, the problem of inaccuracy caused by matching only the longitude and the latitude of the communication cell and the geographical area in the prior art is avoided, and the more accurate area matching relationship between the Thiessen cell and the geographical area corresponding to the communication cell can be obtained.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a schematic flowchart of a method for apportionment mapping of communication cells and geographic areas according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a Delaunay triangle formed by discrete points;
FIG. 3 is a schematic diagram of a Delaunay triangle with the same vertex;
fig. 4 is a flowchart illustrating a method for apportionment mapping of communication cells and geographic areas according to another embodiment of the present invention;
FIG. 5 is a schematic diagram illustrating a reason for the first latitude and longitude correction according to an embodiment of the present invention;
fig. 6 is a schematic diagram of selecting modified longitude and latitude of a triangular sub-region according to an embodiment of the present invention;
fig. 7 is a schematic diagram of selecting modified longitude and latitude of a sub-area in a non-triangular shape according to an embodiment of the present invention;
fig. 8 is a schematic diagram illustrating a reason why a gap is generated when the taison cell layer is cut into the grid layer according to the embodiment of the present invention;
fig. 9 is a schematic diagram of a reason for second longitude and latitude correction according to the embodiment of the present invention;
fig. 10 is a schematic structural diagram of an apparatus for apportionment mapping of communication cells and geographic areas according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It should be noted that, in the embodiments of the present invention, words such as "exemplary" or "for example" are used to indicate examples, illustrations or explanations. Any embodiment or design described as "exemplary" or "e.g.," an embodiment of the present invention is not necessarily to be construed as preferred or advantageous over other embodiments or designs. Rather, use of the word "exemplary" or "such as" is intended to present concepts related in a concrete fashion.
It should be noted that, in the embodiments of the present invention, "of", "corresponding" and "corresponding" may be sometimes used in combination, and it should be noted that, when the difference is not emphasized, the intended meaning is consistent.
For better setting investment plans in different geographic areas, operators need to know the size of the traffic corresponding to each geographic area, and the size of the traffic exists only in the traffic data of the communication cell, so that the communication cell and the geographic area need to be matched firstly. When the currently and generally used communication cell is matched with a geographic area, the geographic area is judged mainly according to the longitude and latitude of the communication cell, taking a grid as an example, the grid to which the longitude and latitude of the communication cell belongs is a corresponding grid; however, in such a case, a small part of the area where the longitude and latitude corresponding position of the communication cell is located is in a certain grid, and the final matching result is that the cell corresponds to the grid, but the grid where the large part of the cell is located is not the grid obtained by matching, which causes inaccuracy of matching between the communication cell and the grid, and accurate traffic in the grid cannot be obtained by using the matching result.
In view of the above problem, referring to fig. 1, an embodiment of the present invention provides a method for apportionment mapping of communication cells and geographic areas, including:
101. geographic data of at least one communication cell is obtained.
102. And carrying out bias processing on the longitude and latitude of the base station according to the longitude and latitude of the base station in the geographic data of the communication cell and the azimuth angle of the communication cell so as to obtain the longitude and latitude of the communication cell.
Specifically, because the longitude and latitude in the geographic data of the communication cell are the longitude and latitude of the base station to which the communication cell belongs, and one base station corresponds to a plurality of cells, if the subsequent processing is not accurate enough by taking the longitude and latitude of the base station as the longitude and latitude of the communication cell, in order to simulate the coverage range of different communication cells under the same base station, the invention adopts the form of the longitude and latitude of the base station plus the direction angle to set the longitude and latitude bias for different communication cells under the same base station to obtain the longitude and latitude of the communication cell, and exemplarily, the longitude and latitude of the base station are the longitude and latitude of the communication cell by biasing the longitude and latitude of the base station according to a certain distance along the direction corresponding to the azimuth angle of the communication cell.
103. And constructing the Thiessen cell corresponding to the communication cell based on the Thiessen polygon principle according to the longitude and latitude of the communication cell.
Specifically, the main steps of the Thiessen polygon are as follows:
(1) construction of Delaulay triangles
In the construction of the Thiessen polygon, discrete points are first constructed into a triangular mesh. Such a triangulation network is called a Delaunay triangulation network.
The construction of the Delaunay triangulation network is also called the construction of an irregular triangulation network, that is, the triangulation network is constructed by discrete data points, as shown in fig. 2, that is, which three data points form a triangle is determined, and the construction is also called the automatic connection triangulation network. I.e. for n discrete points on the plane, the plane coordinate is (x)i,yi) And i is 1,2, …, n, and the three similar points form an optimal triangle, so that each discrete point becomes the vertex of the triangle.
The result of automatically joining the triangulation is the labels of the three vertices of all triangles, e.g. 1,2, 8; 2,8, 3; 3,8, 7; … …
In order to obtain the optimal triangle, when constructing the triangulation, the triangle's interior three corners should be as sharp as possible, i.e. meet the criteria generated by the Delaunay triangle:
any other discrete point cannot be contained within the circumscribed circle of any one Delaunay triangle.
Two adjacent Delaunay triangles form a convex quadrangle, and after the diagonals of the convex quadrangle are exchanged, the minimum of six internal angles is not increased, and the property is the minimum angle maximum criterion.
(2) The numbers of all triangles adjacent to each discrete point are found and recorded. This step requires finding all triangles in the constructed triangulation network that have one and the same vertex.
The triangles adjacent to each discrete point are sorted in a clockwise or counterclockwise direction for next step connection to generate a Thiessen polygon. The method of sorting may be as shown in fig. 3. Let the discrete point be o. Finding out a triangle with o as a vertex and setting the triangle as A; taking another vertex of the triangle A except o as a, and finding out another vertex as f; the next triangle must be bounded by of, which is triangle F; the other vertex of the triangle F is e, and the next triangle takes oe as the side; this is repeated until the oa edge is reached.
(3) And calculating the center of a circumscribed circle of each triangle and recording the center.
(4) And connecting the centers of the circumscribed circles of the adjacent triangles according to the adjacent triangles of each discrete point to obtain the Thiessen polygon.
For example, based on the above principle, in the embodiment of the present invention, the longitude and latitude of all the communication cells correspond to one discrete point, the discrete point corresponding to the longitude and latitude of the communication cell of the thieson cell needs to be established is taken as a vertex, then the vertex and the discrete points around the vertex form a plurality of triangles, a thieson polygon is established based on the distance, and the range corresponding to the thieson polygon is the thieson cell corresponding to the communication cell.
104. Geographic data for at least one grid in a geographic area is obtained.
105. The grid is divided into at least one sub-area according to the boundary data in the geographical data of the Thiessen cell and the geographical data of the grid.
Specifically, it is desirable that both the grid layer and the thiessen cell layer are possible regardless of who is the boundary layer or the bottom map layer. However, in practice, the thiessen cell map is a thiessen polygon automatically generated by MapInfo (desktop geographic information system software), and the format is more standard compared with the case that the grid map layers are manually drawn, so that a large gap and an overlapping area cannot be generated, and non-polygon objects (such as a broken line object and a character object) cannot be contained. Because when these irregular objects are used as cutting boundaries, abnormal errors can occur in the segmentation process. Therefore, in this regard, in the embodiment of the present invention, the grid layer is used as the base map, and the tesson cell layer is used as the cutting boundary layer.
106. And calculating the area of the sub-region and the longitude and latitude of the centroid of the sub-region according to the geographic data of the sub-region.
107. And matching the sub-area with the Thiessen cell according to the longitude and latitude of the centroid of the sub-area and the geographic data of the Thiessen cell.
Because the specific process of matching in the embodiment of the invention is performed in the existing map processing software, and the existing software (MapInfo) only uses the boundary data of the thieson cells to cut the grid when cutting the grid, although it can be seen on the visual interface according to the superposition of the grid layer and the layer of the thieson cells that the sub-region specifically corresponds to, the software itself needs to be identified again, and in order to reduce the complexity of data processing, the longitude and latitude of the centroid of the sub-region are generally firstly found, and then the sub-region and the thieson cells are matched according to the longitude and latitude of the centroid of the sub-region and the geographic data of the thieson cells.
108. And calculating the area of the Thiessen cell according to the geographic data of the Thiessen cell.
109. Calculating the ratio of the area of each subregion matched with the Thiessen cell to the area of the Thiessen cell.
Specifically, the area of each sub-area matched with the thieson cell accounts for the area of the thieson cell, and the area is matched with the communication cell and the grid, and when the traffic in the grid needs to be calculated, the following formula is adopted to substitute the area for the area:
the traffic volume in the target grid is the ratio of the area of the Thiessen cell 1 traffic volume subregion 1 to the area of the Thiessen cell 1 + the ratio of the area of the Thiessen cell 2 traffic volume subregion 2 to the area of the Thiessen cell 2+ … … the ratio of the area of the Thiessen cell n traffic volume subregion n to the area of the Thiessen cell n;
wherein, subregion 1 to subregion n all belong to the target grid, subregion 1 matches with Thiessen district 1, subregion 2 matches with Thiessen district 2, … …, subregion n matches with Thiessen district n.
The method for apportionment mapping of communication cells and geographic areas provided by the above embodiments is that when the communication cells and the geographic areas are matched in each part, firstly, the longitude and latitude of a communication cell are obtained according to the longitude and latitude and the azimuth angle of a base station, then a Thiessen cell image layer of the communication cell is obtained according to the longitude and latitude of the communication cell, then mutually overlapping and cutting the Thiessen cell layer and the geographical area layer so as to divide the geographical area layer into a plurality of sub-areas by the boundary in the Thiessen cell layer, then calculating the area ratio of the plurality of sub-areas to the area of different Thiessen cells which are respectively correspondingly matched, so as to obtain the area matching relationship with more than one area between the Thiessen cell layer and the geographical area layer, then, the operator can obtain the traffic of different geographical areas according to the area matching relationship of more than one area between the Thiessen cell layer and the geographical area layer. In the technical scheme provided by the implementation of the invention, the geographical area map layer is divided by using the area boundary in the Thiessen cell map layer of the communication cell, and the ratio of each divided sub-area to the area of the respectively matched Thiessen cell is obtained, so that the matching relationship between the Thiessen cell and the geographical area grid of each communication cell is accurately obtained from the level of the area ratio, the problem of inaccuracy caused by matching only the longitude and the latitude of the communication cell and the geographical area in the prior art is avoided, and the more accurate area matching relationship between the Thiessen cell and the geographical area corresponding to the communication cell can be obtained.
Referring to fig. 4, an embodiment of the present invention further provides another method for apportionment mapping of communication cells and geographic areas, which is a supplement to the technical solution provided by the foregoing embodiment, and specifically includes:
401. geographic data of at least one communication cell is obtained.
402. And carrying out bias processing on the longitude and latitude of the base station according to the longitude and latitude of the base station in the geographic data of the communication cell and the azimuth angle of the communication cell so as to obtain the longitude and latitude of the communication cell.
403. And constructing the Thiessen cell corresponding to the communication cell based on the Thiessen polygon principle according to the longitude and latitude of the communication cell.
404. Geographic data for at least one grid in a geographic area is obtained.
405. The grid is divided into at least one sub-area according to the boundary data in the geographical data of the Thiessen cell and the geographical data of the grid.
406. And calculating the area of the sub-region and the longitude and latitude of the centroid of the sub-region according to the geographic data of the sub-region.
407. And judging whether a first target sub-area exists in the sub-areas, wherein the center of mass of the first target sub-area is positioned outside or on the boundary of the first target sub-area.
If so, execute 408; if not, 409 is executed;
specifically, for some special polygons, the calculated centroid may be located exactly on the edge line of the polygon, and then when the relationship between the centroid and the thieson cell is matched, not only the points inside the polygon but also the points on the edge line are included. At this time, two adjacent sub-regions on the edge correspondingly contain the particle, and a matching condition that one particle corresponds to two sub-regions occurs. As shown in fig. 5, the dotted line is a cell boundary, the solid line is a grid boundary, the dots are longitude and latitude points of the centroid of the sub-region, and the "sub-region 1 centroid" appears at the boundary line between the "sub-region 1" and the "sub-region 3" in the graph, so that the software is not easily identified;
therefore, it is necessary to determine whether the sub-region has the condition, that is, whether the longitude and latitude of the first time is wrong, and correct the longitude and latitude of the centroid of the sub-region in the condition.
408. And correcting the longitude and latitude of the mass center of the first target sub-area.
Step 407 is executed after step 408;
optionally, as shown in fig. 6, when the first target sub-region is a triangle, modifying the longitude and latitude of the centroid of the first target sub-region includes: and solving the longitude and latitude of any middle point of the first target sub-region, namely the black point in the figure 6, according to the geographic data of the first target sub-region, and taking the longitude and latitude of any middle point of the first target sub-region as the longitude and latitude of the corrected centroid of the first target sub-region.
Optionally, as shown in fig. 7, when the first target sub-region is non-triangular, modifying the longitude and latitude of the centroid of the first target sub-region includes:
sequentially selecting vertexes from all vertexes of the first target sub-area until the longitude and latitude of the midpoint of the target are positioned in the first target sub-area; the target midpoint is a midpoint of a connecting line between a circle which takes the selected vertex as a center and takes the preset length as a radius and an intersection point of two adjacent edges of the selected vertex, namely black points ((x2+ x3)/2, (y2+ y3)/2) in the figure; taking the longitude and latitude of the midpoint of the target as the longitude and latitude of the centroid corrected by the first target subregion; illustratively, the preset length is 5 meters.
409. And matching the sub-area with the Thiessen cell according to the longitude and latitude of the centroid of the sub-area and the geographic data of the Thiessen cell.
410. And calculating the ratio of the sum of the areas of all the sub-areas successfully matched with the Thiessen cell to the area of the Thiessen cell.
411. And judging whether the ratio of the sum of the areas of all the sub-areas successfully matched with the Thiessen cell to the area of the Thiessen cell is greater than a preset percentage.
If yes, go to 412; if not, then 413 is executed; preferably, the preset percentage is 101%;
specifically, in practice, all points corresponding to the longitude and latitude are arranged in longitude and latitude coordinates, and a longitude and latitude coordinate system is measured by taking the latitude as a unit, and only 6 effective digits are reserved after decimal points, so that in the division calculation, vertex coordinates of a boundary formed after division are passivated, a new boundary after division is not completely superposed with the division boundary, a small gap is reserved between the new boundary and the division boundary, as shown in an enlarged diagram of a coordinate lattice shown in fig. 8, black points are longitude and latitude coordinate lattices, the horizontal and vertical intervals are all 0.000001 degrees which is the minimum unit of the longitude and latitude coordinate system, and the calculated longitude and latitude are fitted on the points in the closest longitude and latitude coordinate lattice due to limited precision; the dotted line is the thiessen cell border; the thick solid line is the grid layer boundary (there will still be a small gap because the grid lines are drawn by hand); the thin solid lines are the boundaries of the partitioned sub-regions. The circle is the intersection of the two boundary lines of the dotted line and the thick line, and the divided sub-regions should generate 1 new vertex at the intersection point (the actual intersection point position in the figure). However, due to the fact that the intersection points are passivated in the precision problem, one point (the intersection point position is calculated in the graph) in the dot matrix is selected nearby to serve as a new vertex, and finally, the thin solid line, the dotted line and the thick solid line are not completely overlapped, and a small gap is generated;
when the centroid longitude and latitude of the sub-area are just located in the gap area, when the relation between the centroid point and the Thiessen cell is matched, attribution errors are caused. As shown in fig. 9, the dotted line is the boundary of the thiessen cell, the solid line is the boundary of the mesh, and the black point is the centroid of the sub-region after segmentation. The left figure shows clearly that "sub-area 2" should belong within the range of "cell B". The centroid is at the cell boundary, and as can be seen from the enlarged right image at the circle, the centroid point of the "subregion 2" is at the left side of the dotted line and falls within the range of the "cell a". The method should belong to the range of the 'cell B', but fall into the range of the 'cell A', and the situation of wrong attribution occurs at gaps;
therefore, it is necessary to perform a determination step 412 for the above situation that the assignment error is easily generated, and when the calculated percentage of the sum of the areas of all sub-areas successfully matched with the thieson cell to the area of the thieson cell is greater than the preset percentage (101%), it indicates that the sub-area included in the adjacent thieson cell belongs to the local thieson cell, and at this time, it is determined that the longitude and latitude are determined to be an error for the second time.
412. And determining that all the sub-areas successfully matched with the Thiessen cell are second target sub-areas, and correcting the longitude and latitude of the mass center of the second target sub-areas.
Step 412 is followed by step 409;
optionally, referring to fig. 6, when the second target sub-region is a triangle, modifying the longitude and latitude of the centroid of the second target sub-region includes:
and solving the longitude and latitude of any middle point of the second target sub-area, namely the black point in the figure 6, according to the geographic data of the second target sub-area, and taking the longitude and latitude of the middle point of the second target sub-area as the longitude and latitude of the corrected centroid of the second target sub-area.
Optionally, referring to fig. 7, when the second target sub-region is non-triangular, modifying the longitude and latitude of the centroid of the second target sub-region includes:
sequentially selecting vertexes from all vertexes of the second target subregion, and stopping until the longitude and latitude of the midpoint of the target are positioned in the second target subregion; the target midpoint is a midpoint of a connecting line of a circle drawn by taking the selected vertex as a circle center and a preset length as a radius and an intersection point of two adjacent edges of the selected vertex, namely a black point in the graph 6; and taking the longitude and latitude of the midpoint of the target as the longitude and latitude of the centroid corrected by the second target subarea.
413. Calculating a target ratio of the area of each sub-region matched with the Thiessen cell to the sum of the areas of all the sub-regions matched with the Thiessen cell;
414. and modifying the area of each sub-area matched with the Thiessen cell according to a preset formula according to the target ratio and the area of the Thiessen cell.
Specifically, in practical applications, in some areas, such as mountains, rivers, lakes, and the like, where there are few users, in order to focus on the distribution of users more intensively, the areas are defined as invalid areas in a grid layer, the graph of the part of the area in the grid layer is hollow, while the tesson cells are seamless and not hollow, so that after the grid is used as a base map for segmentation, the tesson cells causing the hollow part to have no grid corresponding to the grids, that is, the sum of the areas of the sub-areas included in the tesson cells near the hollow areas cannot reach 100%, in order to avoid this problem, we need to perform further processing, and redistribute the area of the cell polygon according to the ratio of the sub-areas in the valid area, and the formula is as follows:
subregion 1 area modification ═ target thieson cell area (subregion 1 area/(subregion 1 area + subregion 2 area + … + subregion n area));
the target Thiessen cell is a Thiessen cell matched with the subarea 1 to the subarea n;
the above formula is specifically implemented by the steps 413 and 414.
415. And calculating the area ratio of each sub-area matched with the Thiessen cell to the area ratio of the Thiessen cell.
The embodiment of the invention provides a method for apportionment mapping of communication cells and geographic areas, which comprises the following steps: acquiring geographic data of at least one communication cell, and performing bias processing on the longitude and latitude of a base station according to the longitude and latitude of the base station in the geographic data of the communication cell and the azimuth angle of the communication cell to acquire the longitude and latitude of the communication cell; building a Thiessen cell corresponding to the communication cell based on the Thiessen polygon principle according to the longitude and latitude of the communication cell; acquiring geographic data of at least one grid in a geographic area, and dividing the grid into at least one sub-area according to boundary data in the geographic data of the Thiessen cell and the geographic data of the grid; calculating the area of the sub-region and the longitude and latitude of the centroid of the sub-region according to the geographic data of the sub-region; matching the sub-area with the Thiessen cell according to the longitude and latitude of the centroid of the sub-area and the geographic data of the Thiessen cell; calculating the area of the Thiessen cell according to the geographic data of the Thiessen cell; calculating the ratio of the area of each subregion matched with the Thiessen cell to the area of the Thiessen cell. When the communication cell and the geographic area are shared and mapped, the longitude and the latitude of the communication cell are obtained according to the longitude and the latitude and the azimuth of the base station, the Thiessen cell layer of the communication cell is obtained according to the longitude and the latitude of the communication cell, the Thiessen cell layer and the geographic area layer are overlapped and cut, so that the geographic area layer is divided into a plurality of sub-areas by the boundary in the Thiessen cell layer, the area occupation ratio values of the sub-areas in different Thiessen cell areas which are correspondingly matched with each other are calculated, the area matching relation with more than one area between the Thiessen cell layer and the geographic area layer can be obtained, and then an operator can obtain the service volume of different geographic areas according to the area matching relation with more than one area between the Thiessen cell layer and the geographic area layer. In the technical scheme provided by the implementation of the invention, the geographical area layer is divided by using the area boundary in the Thiessen cell layer of the communication cell, and the ratio of each divided sub-area to the area of the respectively matched Thiessen cell is obtained, so that the apportionment mapping relation between the Thiessen cell and the geographical area grid of each communication cell is accurately obtained from the level of the area ratio, the problem of inaccuracy caused by matching the communication cell longitude and latitude and the geographical area in the prior art is avoided, and the more accurate area matching relation between the Thiessen cell and the geographical area corresponding to the communication cell can be obtained.
Referring to fig. 10, an embodiment of the present invention further provides an apparatus 01 for apportionment mapping of communication cells and geographic areas, including:
an acquisition module 011, a Thiessen cell module 012, a segmentation module 013, a calculation module 014 and a matching module 015;
the acquisition module 011 is used for acquiring the geographic data of at least one communication cell and offsetting the longitude and latitude of the base station according to the longitude and latitude of the base station in the geographic data of the communication cell and the azimuth angle of the communication cell so as to acquire the longitude and latitude of the communication cell;
the thiessen cell module 012 is configured to construct a thiessen cell corresponding to the communication cell based on the thiessen polygon principle according to the longitude and latitude of the communication cell acquired by the acquisition module 011;
the obtaining module 011 is further configured to obtain geographic data of at least one grid in the geographic region, and the dividing module 013 is configured to divide the grid into at least one sub-region according to the boundary data in the geographic data of the thiessen cells constructed by the thiessen cell module 012 and the geographic data of the grid obtained by the obtaining module 011;
a calculating module 014, configured to calculate the area of the sub-region and the longitude and latitude of the centroid of the sub-region according to the geographic data of the sub-region divided by the dividing module 013;
a matching module 015, configured to match the sub-area with the thieson cell according to the longitude and latitude of the centroid of the sub-area calculated by the calculation module 014 and the geographic data of the thieson cell constructed by the thieson cell module 012;
the calculation module 014 is further configured to calculate the area of the thiessen cell according to the geographical data of the thiessen cell constructed by the thiessen cell module 012;
calculation module 014 is also used to calculate the ratio of the area of each sub-area matched by matching module 015 for a thieson cell to the area of the thieson cell.
It should be noted that, in the embodiment of the present invention, the apportionment mapping apparatus uses MapInfo written in mapsubac language to work on the basis of the existing MapInfo software, so as to implement the work of the modules and implement information flow interaction work among the modules.
Optionally, the apparatus further includes a determining module 016 and a modifying module 017;
the judging module 016 is configured to judge whether the sub-region divided by the dividing module 013 has a first target sub-region before the matching module 015 matches the sub-region with the thiessen cell according to the longitude and latitude of the centroid of the sub-region calculated by the calculating module 014 and the geographic data of the thiessen cell constructed by the thiessen cell module 012; the centroid of the first target sub-region lies outside or on the boundary of the first target sub-region;
when the determining module 016 determines that the first target sub-region exists in the sub-region, the correcting module 017 is configured to correct the longitude and latitude of the centroid of the first target sub-region.
Optionally, the calculating module 014 is further configured to calculate a ratio of the sum of the areas of all sub-areas successfully matched with the thieson cell to the area of the thieson cell before the area ratio of each sub-area matched with the thieson cell by the matching module 015 to the area of the thieson cell;
when the determining module 016 determines that the ratio of the sum of the areas of all sub-regions successfully matched with the thieson cell to the area of the thieson cell is greater than the preset percentage, it is determined that all sub-regions successfully matched with the thieson cell are the second target sub-region, and the correcting module 017 is further used for correcting the second target sub-region;
the matching module 015 is further configured to match the sub-area with the thieson cell again according to the longitude and latitude of the centroid of the sub-area and the geographic data of the thieson cell;
the calculating module 014 is further configured to recalculate the ratio of the sum of the areas of all sub-areas successfully matched with the thieson cell to the area of the thieson cell.
Optionally, the calculating module 014 is further configured to calculate an area ratio of the area of each sub-region matched with the thieson cell to the thieson cell when the determining module 016 determines that the area ratio of the sum of the areas of all sub-regions successfully matched with the thieson cell to the area of the thieson cell is not greater than a preset percentage.
Optionally, an apportioning module 018 is further included;
the calculating module 014 is further configured to calculate a target occupation ratio value of the area of each sub-region matching the thieson cell to the sum of the areas of all sub-regions matching the thieson cell before calculating the area occupation ratio value of each sub-region matching the thieson cell;
an apportioning module 018 is configured to modify the area of each sub-region matched with the thieson cell according to a preset formula according to the target occupancy value calculated by the calculating module 014 and the area of the thieson cell.
Optionally, when the first target sub-region is a triangle, the modification module 017 is specifically configured to:
and solving the longitude and latitude of any central line midpoint of the first target sub-region according to the geographic data of the first target sub-region, and taking the longitude and latitude of any central line midpoint of the sub-region as the longitude and latitude of the corrected centroid of the sub-region.
Optionally, when the first target sub-region is non-triangular, the modification module 017 is specifically configured to:
sequentially selecting vertexes from all vertexes of the first target subarea, and stopping until the longitude and latitude of the midpoint of the target are positioned in the subarea; the target midpoint is the midpoint of a connecting line of intersection points of a circle drawn by taking the selected vertex as the center of a circle and two adjacent edges of the selected vertex, wherein the preset length is the radius;
and taking the longitude and latitude of the midpoint of the target as the longitude and latitude of the centroid corrected by the first target subregion.
Optionally, when the second target sub-region is a triangle, the modification module 017 is specifically configured to:
and solving the longitude and latitude of any central line midpoint of the second target sub-region according to the geographic data of the second target sub-region, and taking the longitude and latitude of any central line midpoint of the sub-region as the longitude and latitude of the corrected centroid of the sub-region.
Optionally, when the second target sub-region is non-triangular, the modification module 017 is specifically configured to:
sequentially selecting vertexes from all vertexes of the second target subarea, and stopping until the longitude and latitude of the midpoint of the target are positioned in the subarea; the target midpoint is the midpoint of a connecting line of intersection points of a circle drawn by taking the selected vertex as the center of a circle and two adjacent edges of the selected vertex, wherein the preset length is the radius;
and taking the longitude and latitude of the midpoint of the target as the longitude and latitude of the centroid corrected by the second target subarea.
For example, taking a grid in the city XX and a thieson cell as an example, through matching by the matching device, the area occupied by each grid in the thieson cell is as shown in the following table 1:
Figure GDA0003093200660000181
TABLE 1
As can be seen from the above table, the area of the cell "768624-51" is 76.78% located within the grid "SX-WL-XA-CA-potential grid-YY street-0986" and 23.22% located within the grid "office-XX tobacco Corp cigarette distribution center", thereby determining that a majority of the users of the cell are distributed within the grid "office-XX tobacco Corp cigarette distribution center".
Further, taking the "office building-XX tobacco company cigarette logistics distribution center" grid as an example, through the calculation of the matching device, the matching conditions of the thiessen cells matched with each sub-area of the grid can be obtained as shown in the following table 2:
Figure GDA0003093200660000182
Figure GDA0003093200660000191
TABLE 2
Then, the business volume of the grid of the cigarette logistics distribution center of office building-XX tobacco company can be obtained according to the contents in table 2, and the specific formula is as follows:
traffic within the grid, 769828-50 cell traffic 33.34% +768624-51 cell traffic 23.22% +768273-50 cell traffic 16.10% +769828-49 cell traffic 0.54% +768624-49 cell traffic 5.02% +769828-51 cell traffic 74.63%;
then the operator can make corresponding investment planning according to the grid traffic data.
The device for apportionment mapping of communication cells and geographic areas provided by the embodiment of the invention comprises: the acquisition module is used for acquiring the geographic data of at least one communication cell and carrying out bias processing on the longitude and latitude of the base station according to the longitude and latitude of the base station in the geographic data of the communication cell and the azimuth angle of the communication cell so as to acquire the longitude and latitude of the communication cell; the Thiessen cell module is used for constructing a Thiessen cell corresponding to the communication cell based on the Thiessen polygon principle according to the longitude and latitude of the communication cell acquired by the acquisition module; the obtaining module is further used for obtaining geographic data of at least one grid in the geographic area, and the dividing module is used for dividing the grid into at least one sub-area according to boundary data in the geographic data of the Thiessen cell constructed by the Thiessen cell module and the geographic data of the grid obtained by the obtaining module; the calculation module is used for calculating the area of the subarea and the longitude and latitude of the centroid of the subarea according to the geographical data of the subarea divided by the division module; the matching module is used for matching the sub-region with the Thiessen cell according to the longitude and latitude of the centroid of the sub-region calculated by the calculating module and the geographic data of the Thiessen cell constructed by the Thiessen cell module; the calculation module is also used for calculating the area of the Thiessen cell according to the geographic data of the Thiessen cell constructed by the Thiessen cell module; the calculation module is further used for calculating the ratio of the area of each sub-area matched for the Thiessen cell by the matching module to the area of the Thiessen cell. Therefore, when the matching device provided by the embodiment of the invention performs apportionment mapping on the communication cell and the geographic area, firstly, the longitude and latitude of the communication cell are obtained according to the longitude and latitude and the azimuth angle of the base station, then, the thieson cell layer of the communication cell is obtained according to the longitude and latitude of the communication cell, then, the thieson cell layer and the geographic area layer are overlapped and cut, so that the geographic area layer is divided into a plurality of sub-areas by the boundary in the thieson cell layer, then, the occupation ratio of the areas of the sub-areas to the areas of different thieson cells which are respectively correspondingly matched is calculated, the area matching relationship with more than one area between the thieson cell layer and the geographic area layer can be obtained, and then, an operator can obtain the service volume of different geographic areas according to the area matching relationship with more than one area between the thieson cell layer and the geographic area layer. In the technical scheme provided by the implementation of the invention, the geographical area layer is divided by using the area boundary in the Thiessen cell layer of the communication cell, and the ratio of each divided sub-area to the area of the respectively matched Thiessen cell is obtained, so that the apportionment mapping result between the Thiessen cell and the geographical area grid of each communication cell is accurately obtained from the level of the area ratio, the problem of inaccuracy in the prior art when the communication cell is matched with the geographical area only by the longitude and latitude of the communication cell is solved, and the more accurate area matching relationship between the Thiessen cell and the geographical area corresponding to the communication cell can be obtained.
The steps of a method or algorithm described in connection with the disclosure herein may be embodied in hardware or in software instructions executed by a processor. Embodiments of the present invention also provide a storage medium, which may include a memory for storing computer software instructions for a method for split mapping of communication cells and geographic areas, including program code designed to perform the method for split mapping of communication cells and geographic areas. Specifically, the software instructions may be composed of corresponding software modules, and the software modules may be stored in a Random Access Memory (RAM), a flash Memory, a Read Only Memory (ROM), an Erasable Programmable Read Only Memory (EPROM), an Electrically Erasable Programmable Read Only Memory (EEPROM), a register, a hard disk, a removable hard disk, a compact disc Read Only Memory (CD-ROM), or any other form of storage medium known in the art. An exemplary storage medium is coupled to the processor such the processor can read information from, and write information to, the storage medium. Of course, the storage medium may also be integral to the processor. The processor and the storage medium may reside in an ASIC. Additionally, the ASIC may reside in a core network interface device. Of course, the processor and the storage medium may reside as discrete components in a core network interface device.
The embodiment of the present invention further provides a computer program, which can be directly loaded into the memory and contains software codes, and after the computer program is loaded and executed by the computer, the method for apportionment mapping of the communication cells and the geographic areas can be implemented.
Those skilled in the art will recognize that, in one or more of the examples described above, the functions described in this invention may be implemented in hardware, software, firmware, or any combination thereof. When implemented in software, the functions may be stored on or transmitted over as one or more instructions or code on a computer-readable medium. Computer-readable media includes both computer storage media and communication media including any medium that facilitates transfer of a computer program from one place to another. A storage media may be any available media that can be accessed by a general purpose or special purpose computer.
The above description is only for the specific embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (16)

1. A method for apportionment mapping of communication cells and geographic areas, comprising:
acquiring geographic data of at least one communication cell, and performing bias processing on the longitude and latitude of a base station according to the longitude and latitude of the base station and the azimuth angle of the communication cell in the geographic data of the communication cell to acquire the longitude and latitude of the communication cell;
building a Thiessen cell corresponding to the communication cell based on the Thiessen polygon principle according to the longitude and latitude of the communication cell;
acquiring geographic data of at least one grid in a geographic area, and dividing the grid into at least one sub-area according to boundary data in the geographic data of the Thiessen cell and the geographic data of the grid;
calculating the area of the sub-region and the longitude and latitude of the centroid of the sub-region according to the geographic data of the sub-region;
judging whether a first target subregion exists in the subregions or not; the centroid of the first target sub-region lies outside or on a boundary of the first target sub-region;
when the first target sub-area exists in the sub-area, correcting the longitude and latitude of the center of mass of the first target sub-area;
matching the sub-area with the Thiessen cell according to the longitude and latitude of the centroid of the sub-area and the geographic data of the Thiessen cell;
calculating the area of the Thiessen cell according to the geographic data of the Thiessen cell;
calculating a ratio of an area of each sub-region matched with the Thiessen cell to an area of the Thiessen cell.
2. The method of claim 1, wherein the calculating an area fraction value of each sub-area matching the Thiessen cell to the area fraction value of the Thiessen cell further comprises:
calculating the ratio of the sum of the areas of all the sub-areas successfully matched with the Thiessen cell to the area of the Thiessen cell;
when the proportion of the sum of the areas of all sub-areas successfully matched with the Thiessen cell to the area of the Thiessen cell is determined to be larger than a preset percentage, determining that all the sub-areas successfully matched with the Thiessen cell are second target sub-areas, and correcting the longitude and latitude of the mass center of the second target sub-areas;
matching the sub-area with the Thiessen cell again according to the longitude and latitude of the centroid of the sub-area and the geographic data of the Thiessen cell;
and recalculating the ratio of the sum of the areas of all the sub-areas successfully matched with the Thiessen cell to the area of the Thiessen cell.
3. The method according to claim 2, wherein the area ratio of each subregion matching the Thiessen cell to the area of the Thiessen cell is calculated when it is determined that the ratio of the sum of the areas of all subregions matching the Thiessen cell successfully to the area of the Thiessen cell is not greater than a preset percentage.
4. The method of claim 3, wherein the calculating an area fraction value of each sub-area matching the Thiessen cell to the area fraction value of the Thiessen cell further comprises:
calculating a target ratio of the area of each sub-area matched with the Thiessen cell to the sum of the areas of all sub-areas matched with the Thiessen cell;
and according to the target occupation value and the area of the Thiessen cell, modifying the area of each sub-area matched with the Thiessen cell according to a preset formula.
5. The method of claim 1, wherein when the first target sub-region is triangular, the modifying the longitude and latitude of the centroid of the first target sub-region comprises:
and solving the longitude and latitude of any middle point of the first target sub-area according to the geographic data of the first target sub-area, and taking the longitude and latitude of any middle point of the first target sub-area as the longitude and latitude of the corrected centroid of the first target sub-area.
6. The method of claim 1, wherein when the first target sub-region is non-triangular, the modifying the longitude and latitude of the centroid of the first target sub-region comprises:
sequentially selecting vertexes from all vertexes of the first target sub-area until the longitude and latitude of the target midpoint are positioned in the first target sub-area; the target midpoint is the midpoint of a connecting line of intersection points of a circle drawn by taking the selected vertex as the center of a circle and two adjacent edges of the selected vertex, wherein the preset length is the radius;
and taking the longitude and latitude of the midpoint of the target as the longitude and latitude of the centroid corrected by the first target subregion.
7. The method of claim 2, wherein when the second target sub-region is triangular, the modifying the longitude and latitude of the centroid of the second target sub-region comprises:
and solving the longitude and latitude of any middle point of the second target sub-area according to the geographic data of the second target sub-area, and taking the longitude and latitude of the middle point of the second target sub-area as the longitude and latitude of the corrected centroid of the second target sub-area.
8. The method of claim 2, wherein when the second target sub-region is non-triangular, the modifying the longitude and latitude of the centroid of the second target sub-region comprises:
sequentially selecting vertexes from all vertexes of the second target subregion, and stopping until the longitude and latitude of the target midpoint are positioned in the second target subregion; the target midpoint is the midpoint of a connecting line of intersection points of a circle drawn by taking the selected vertex as the center of a circle and two adjacent edges of the selected vertex, wherein the preset length is the radius;
and taking the longitude and latitude of the midpoint of the target as the longitude and latitude of the centroid corrected by the second target subregion.
9. An apparatus for apportionment mapping of communication cells and geographic areas, comprising: the system comprises an acquisition module, a Taisen cell module, a segmentation module, a calculation module, a matching module, a judgment module and a correction module;
the acquisition module is used for acquiring the geographic data of at least one communication cell and carrying out bias processing on the longitude and latitude of the base station according to the longitude and latitude of the base station in the geographic data of the communication cell and the azimuth angle of the communication cell so as to acquire the longitude and latitude of the communication cell;
the Thiessen cell module is used for constructing a Thiessen cell corresponding to the communication cell based on the Thiessen polygon principle according to the longitude and latitude of the communication cell acquired by the acquisition module;
the obtaining module is further configured to obtain geographic data of at least one grid in a geographic area, and the dividing module is configured to divide the grid into at least one sub-area according to boundary data in the geographic data of the Thiessen cell constructed by the Thiessen cell module and the geographic data of the grid obtained by the obtaining module;
the calculation module is used for calculating the area of the sub-region and the longitude and latitude of the centroid of the sub-region according to the geographic data of the sub-region divided by the division module;
the judging module is used for judging whether a first target sub-region exists in the sub-regions divided by the dividing module; the centroid of the first target sub-region lies outside or on a boundary of the first target sub-region;
when the judgment module determines that the first target sub-region exists in the sub-region, the correction module is used for correcting the longitude and latitude of the mass center of the first target sub-region;
the matching module is used for matching the sub-area with the Thiessen cell according to the longitude and latitude of the centroid of the sub-area calculated by the calculating module and the geographic data of the Thiessen cell constructed by the Thiessen cell module;
the calculation module is further used for calculating the area of the Thiessen cell according to the geographic data of the Thiessen cell constructed by the Thiessen cell module;
the calculating module is further used for calculating the ratio of the area of each sub-area matched for the Thiessen cell by the matching module to the area of the Thiessen cell.
10. The apparatus of claim 9,
before the calculation module calculates the area of each sub-region matched with the Thiessen cell by the matching module to the area of the Thiessen cell, the calculation module is further used for calculating the area of the sum of the areas of all the sub-regions successfully matched with the Thiessen cell to the area of the Thiessen cell;
when the judging module determines that the ratio of the sum of the areas of all sub-areas successfully matched with the Thiessen cell to the area of the Thiessen cell is greater than a preset percentage, determining that all the sub-areas successfully matched with the Thiessen cell are second target sub-areas, and the correcting module is further used for correcting the second target sub-areas;
the matching module is further used for matching the sub-area with the Thiessen cell again according to the longitude and latitude of the centroid of the sub-area and the geographic data of the Thiessen cell;
the calculation module is further used for recalculating the ratio of the sum of the areas of all the sub-areas successfully matched with the Thiessen cell to the area of the Thiessen cell.
11. The apparatus of claim 10, wherein the calculating module is further configured to calculate an area ratio of each sub-region matching the Thiessen cell to the Thiessen cell when the determining module determines that a ratio of a sum of areas of all sub-regions matching the Thiessen cell to the area of the Thiessen cell is not greater than a preset percentage.
12. The apparatus of claim 11, further comprising a apportioning module;
the calculation module is further configured to calculate a target ratio of the area of each sub-region matched with the Thiessen cell to the sum of the areas of all sub-regions matched with the Thiessen cell before calculating the area ratio of each sub-region matched with the Thiessen cell to the area ratio of the Thiessen cell;
the allocation module is used for modifying the area of each sub-area matched with the Thiessen cell according to a preset formula according to the target occupation value calculated by the calculation module and the area of the Thiessen cell.
13. The apparatus according to claim 10, wherein when the first target sub-region is triangular, the modification module is specifically configured to:
and solving the longitude and latitude of any central line midpoint of the first target sub-region according to the geographic data of the first target sub-region, and taking the longitude and latitude of any central line midpoint of the sub-region as the longitude and latitude of the corrected centroid of the sub-region.
14. The apparatus of claim 10, wherein when the first target sub-region is non-triangular, the modification module is specifically configured to:
sequentially selecting vertexes from all vertexes of the first target sub-area until the longitude and latitude of the target midpoint are positioned in the sub-area; the target midpoint is the midpoint of a connecting line of intersection points of a circle drawn by taking the selected vertex as the center of a circle and two adjacent edges of the selected vertex, wherein the preset length is the radius;
and taking the longitude and latitude of the midpoint of the target as the longitude and latitude of the centroid corrected by the first target subregion.
15. The apparatus according to claim 11, wherein when the second target sub-region is triangular, the modification module is specifically configured to:
and solving the longitude and latitude of any central line midpoint of the second target subarea according to the geographic data of the second target subarea, and taking the longitude and latitude of any central line midpoint of the subarea as the longitude and latitude of the corrected centroid of the subarea.
16. The apparatus according to claim 11, wherein when the second target sub-region is non-triangular, the modification module is specifically configured to:
sequentially selecting vertexes from all vertexes of the second target subarea until the longitude and latitude of the target midpoint are positioned in the subarea; the target midpoint is the midpoint of a connecting line of intersection points of a circle drawn by taking the selected vertex as the center of a circle and two adjacent edges of the selected vertex, wherein the preset length is the radius;
and taking the longitude and latitude of the midpoint of the target as the longitude and latitude of the centroid corrected by the second target subregion.
CN201810273229.0A 2018-03-29 2018-03-29 Allocation mapping method and device for communication cells and geographic areas Active CN108601031B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810273229.0A CN108601031B (en) 2018-03-29 2018-03-29 Allocation mapping method and device for communication cells and geographic areas

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810273229.0A CN108601031B (en) 2018-03-29 2018-03-29 Allocation mapping method and device for communication cells and geographic areas

Publications (2)

Publication Number Publication Date
CN108601031A CN108601031A (en) 2018-09-28
CN108601031B true CN108601031B (en) 2021-07-27

Family

ID=63624999

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810273229.0A Active CN108601031B (en) 2018-03-29 2018-03-29 Allocation mapping method and device for communication cells and geographic areas

Country Status (1)

Country Link
CN (1) CN108601031B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109918469B (en) * 2019-03-27 2021-02-23 中国联合网络通信集团有限公司 Gridding processing method and device
CN110300415A (en) * 2019-07-27 2019-10-01 北京智数时空科技有限公司 A method of determining wireless network coverage in bounded domain
CN112686483A (en) * 2019-10-17 2021-04-20 ***通信集团陕西有限公司 Early warning area identification method and device, computing equipment and computer storage medium
CN114360255B (en) * 2022-03-21 2022-08-30 北京闪马智建科技有限公司 Flow determination method and device, storage medium and electronic device
CN115952427B (en) * 2023-03-14 2023-05-30 山东美航天天能源技术有限公司 Industrial garden digital operation management method and system

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101203015A (en) * 2006-12-14 2008-06-18 北京亿阳信通软件研究院有限公司 Method and apparatus for determining base station wireless field strength overlay area range
WO2011146180A1 (en) * 2010-05-19 2011-11-24 Alcatel-Lucent Usa Inc. Methods of determining coverage areas
CN103024756A (en) * 2011-09-20 2013-04-03 北京电旗通讯技术股份有限公司 Method and system for planning entire network community coverage
EP2869622A1 (en) * 2013-11-05 2015-05-06 Telefonica Digital España, S.L.U. Method and device for locating network activity in cellular communication networks
CN106470427A (en) * 2015-08-20 2017-03-01 ***通信集团黑龙江有限公司 A kind of partitioning method and device of cell scenario
CN106998556A (en) * 2017-03-17 2017-08-01 北京邮电大学 A kind of system and implementation method for building wireless network traffic distribution map
CN107396374A (en) * 2017-07-07 2017-11-24 江苏奥斯威尔信息科技有限公司 A kind of covering method based on fictitious force and Thiessen polygon

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101203015A (en) * 2006-12-14 2008-06-18 北京亿阳信通软件研究院有限公司 Method and apparatus for determining base station wireless field strength overlay area range
WO2011146180A1 (en) * 2010-05-19 2011-11-24 Alcatel-Lucent Usa Inc. Methods of determining coverage areas
CN103024756A (en) * 2011-09-20 2013-04-03 北京电旗通讯技术股份有限公司 Method and system for planning entire network community coverage
EP2869622A1 (en) * 2013-11-05 2015-05-06 Telefonica Digital España, S.L.U. Method and device for locating network activity in cellular communication networks
CN106470427A (en) * 2015-08-20 2017-03-01 ***通信集团黑龙江有限公司 A kind of partitioning method and device of cell scenario
CN106998556A (en) * 2017-03-17 2017-08-01 北京邮电大学 A kind of system and implementation method for building wireless network traffic distribution map
CN107396374A (en) * 2017-07-07 2017-11-24 江苏奥斯威尔信息科技有限公司 A kind of covering method based on fictitious force and Thiessen polygon

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
A grid-based energy saving scheme with traffic map in heterogeneous dense 5G network;Fangmin Xu等;《 2017 17th International Symposium on Communications and Information Technologies》;20180118;全文 *
Voronoi在4G网络容量规划中的应用及研究;黄小光等;《邮电设计技术》;20170731;参见第2节,图1 *

Also Published As

Publication number Publication date
CN108601031A (en) 2018-09-28

Similar Documents

Publication Publication Date Title
CN108601031B (en) Allocation mapping method and device for communication cells and geographic areas
US5796634A (en) System and method for identifying the geographic region of a geographic area which contains a geographic zone associated with a location
US5961572A (en) System and method for identifying the geographic region of a geographic area which contains a geographic point associated with a location
US20030074132A1 (en) Polar coordinate-based isochrone generation
CN109714784B (en) Antenna azimuth angle optimization method and device
CN115379464B (en) Planning method, device and equipment of base station and computer readable storage medium
CN111935733B (en) Method, device and equipment for checking base station position information and computer storage medium
US8874375B2 (en) Digital map projection
CN110136174A (en) A kind of target object tracking and device
US20180199154A1 (en) Retrospective path analysis
CN114067077A (en) Method and system for accurately measuring average rainfall of regional government regions in water conservancy industry
CN115388873A (en) Positioning detection method and device, electronic equipment and storage medium
CN113670338A (en) Error measurement method and device, electronic equipment and readable storage medium
CN115474206A (en) Real-time people number determination method and device, electronic equipment and storage medium
CN112020091B (en) Weak coverage processing method and device based on digitized geographic information
CN112184900A (en) Method and device for determining elevation data and storage medium
CN110989886A (en) Three-dimensional space grid selection method and device based on space map
CN116030116A (en) Landfill volume analysis method and device, electronic equipment and storage medium
CN114387419B (en) Three-dimensional geographic entity coding method and device based on multi-dimensional space subdivision
CN110334311A (en) The area computation method of Subgraph spot on a kind of earth ellipsoid face
CN115935561A (en) Three-dimensional pipe network modeling method and device, electronic equipment and storage medium
CN108731648A (en) 2000 coordinate system parameter acquiring methods, device and computer readable storage medium
CN114333433B (en) Method for determining sector of aircraft based on geodetic coordinate system
CN114238537B (en) House general survey data grouping processing method and system
CN114173363B (en) Method and device for detecting TAC planning abnormity and storage medium

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant