A coverage area calculation method and system, electronic equipment and medium
By using Green's formula and the intersection operation method of sorting, the boundary curve of the intersection of the protection area and the force coverage area is described, which solves the problems of slow calculation accuracy and speed in the existing technology and realizes fast and accurate coverage area calculation.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- ROCKET FORCE UNIV OF ENG
- Filing Date
- 2022-12-22
- Publication Date
- 2026-07-21
AI Technical Summary
Existing technologies are not accurate and slow when calculating the area covered by support forces, especially when evaluating the effectiveness of multiple support force deployments over a large scale, and cannot meet practical needs.
Green's formula is used to calculate the area of the protection force coverage area. The boundary curves of the intersection of the protection area and the force coverage area are described in two ways: the boundary curve segments of the protection force coverage circle are described by angles, while the boundary curves on the protection area are described by the corresponding proportions of the start and end points. By combining the sorted intersection operation method, all the boundary curves on the circumference of the protection circle within the protection area are obtained.
It improves the calculation efficiency of the coverage area of the support force, with short calculation time and high solution quality, and can quickly and accurately calculate the coverage area of multiple support forces.
Smart Images

Figure CN116383548B_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of operations management technology, and in particular to a method, system, electronic device, and medium for calculating coverage area. Background Technology
[0002] During the deployment of support forces, such as communication base stations and emergency support warehouses, simulation methods are often used to evaluate the effectiveness of emergency support warehouses and communication base stations. This involves generating numerous simulations to determine the locations of deployment points and calculating the area covered by these points within a given support area. This involves rapidly calculating the area covered by the support forces given the number of support forces and their radius of influence, serving as a basis for evaluating the deployment plan's effectiveness. The computational load required for this coverage problem increases dramatically with the number of support forces. The simplest method for quickly calculating the coverage area is to grid the area covered by the support forces and statistically analyze the area of the covered grid regions. However, this approach is not very accurate or fast, and is only a simple method for calculating the coverage area within a small area. When using simulation methods to evaluate the deployment effectiveness of multiple support forces over a larger support area, using Green's formula to calculate the coverage area is more accurate and effective than the grid method. However, while Green's formula is effective, automatically determining the boundary curve of the coverage area is a complex problem. Existing methods for calculating the boundary curve of the coverage area are time-consuming and cannot meet the needs of practical work. Therefore, there is an urgent need for a more efficient method for calculating the coverage area of support forces. Summary of the Invention
[0003] The purpose of this invention is to provide a method, system, electronic device, and medium for calculating coverage area, which can improve the calculation efficiency of the coverage area of the protection force.
[0004] To achieve the above objectives, the present invention provides the following solution:
[0005] A method for calculating coverage area, comprising:
[0006] Obtain the coordinates of all vertices on the boundary of the protection area, as well as the center and radius of the protection circle corresponding to each protection force within the protection area; the protection circle corresponding to each protection force is constructed with the deployment position of the protection force as the center and the radius of action of the protection force as the radius;
[0007] For any given protection circle, a first polar angle set corresponding to the protection circle is calculated based on the center and radius of the protection circle and the centers and radii of the remaining protection circles. The first polar angle set includes the polar angles of each intersection point in the first intersection point set corresponding to the protection circle in the target polar coordinate system. The remaining protection circle set includes all protection circles excluding the protection circle. The target polar coordinate system is constructed with the center of the protection circle as its pole, and two intersecting protection circles have two intersection points. The first intersection point set corresponding to the protection circle includes the intersection points of the protection circle with each of the remaining protection circles.
[0008] The second polar angle set of the protection circle is calculated based on the center and radius of the protection circle and the endpoint coordinates of all line segments on the boundary of the protection area. The starting point ratio and ending point ratio of each line segment when the protection circle intersects with all line segments on the boundary of the protection area are also calculated. The second polar angle set includes the starting point polar angle and ending point polar angle of the protection circle in the target polar coordinate system when the protection circle intersects with all line segments on the boundary of the protection area.
[0009] The starting coordinates and ending coordinates of all line segments on the boundary of the protected area are determined based on the starting and ending proportions of each line segment when all the protected circles intersect with all line segments on the boundary of the protected area.
[0010] The first polar angle set and the second polar angle set corresponding to the protection circle are merged and subtracted to obtain the total polar angle set corresponding to the protection circle.
[0011] The coverage area of the protection force over the protection area is calculated based on the radius and center of all protection circles, the total set of polar angles corresponding to all protection circles, and the starting and ending coordinates of all line segments on the boundary of the protection area.
[0012] A system for calculating coverage area, comprising:
[0013] The acquisition module is used to acquire the coordinates of all vertices on the boundary of the protection area and the center and radius of the protection circle corresponding to each protection force within the protection area; the protection circle corresponding to each protection force is constructed with the deployment position of the protection force as the center and the radius of action of the protection force as the radius;
[0014] The first polar angle set determination module is used to calculate, for any given protection circle, the first polar angle set corresponding to the protection circle based on the center and radius of the protection circle and the centers and radii of each protection circle in the remaining protection circle set; the first polar angle set includes the polar angles of each intersection point in the first intersection point set corresponding to the protection circle in the target polar coordinate system; the remaining protection circle set includes all protection circles excluding the protection circle; the target polar coordinate system is constructed with the center of the protection circle as the pole, and two intersecting protection circles have two intersection points; the first intersection point set corresponding to the protection circle includes the intersection points of the protection circle with each protection circle in the remaining protection circle set;
[0015] The second polar angle set and the proportion determination module are used to calculate the second polar angle set of the protection circle and the starting and ending proportions of each line segment when the protection circle intersects with all line segments on the boundary of the protection area, based on the center and radius of the protection circle and the endpoint coordinates of all line segments on the boundary of the protection area; the second polar angle set includes the starting and ending polar angles of the protection circle in the target polar coordinate system when the protection circle intersects with all line segments on the boundary of the protection area respectively;
[0016] The start and end coordinate determination module is used to determine the start and end coordinates of all line segments on the boundary of the protection area based on the start and end ratios of each line segment when all protection circles intersect with all line segments on the boundary of the protection area.
[0017] The deletion module is used to merge and delete the first polar angle set and the second polar angle set corresponding to the protection circle to obtain the total polar angle set corresponding to the protection circle;
[0018] The coverage area calculation module is used to calculate the coverage area of the protection force covering the protection area based on the radius and center of all protection circles, the total set of polar angles corresponding to all protection circles, and the starting and ending coordinates of all line segments on the boundary of the protection area.
[0019] An electronic device, comprising:
[0020] A memory and a processor, the memory for storing a computer program, the processor for running the computer program to cause the electronic device to perform the coverage area calculation method as described above.
[0021] A computer-readable storage medium storing a computer program that, when executed by a processor, implements the method for calculating the coverage area as described above.
[0022] According to specific embodiments provided by the present invention, the present invention discloses the following technical effects: The present invention uses Green's formula to calculate the coverage area of the protection force in the protection area, and uses two methods to describe the boundary curve of the intersection of the protection area and the coverage area of the protection force: the boundary curve segment of the protection force coverage circle is described by angle, while the boundary curve on the protection area is described by the corresponding proportion of the start and end points. The intersection operation method based on sorting is used to obtain all the boundary curves on the protection circle within the protection area, and the union operation method based on sorting is used to obtain all the boundary curves of the protected part on the boundary of the protection area, thereby improving the calculation efficiency of the area covered by the protection force in the protection area. Attached Figure Description
[0023] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.
[0024] Figure 1 A flowchart illustrating a method for calculating the area of a protection force coverage zone, provided as an embodiment of the present invention. Detailed Implementation
[0025] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.
[0026] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.
[0027] This invention provides a method for calculating coverage area, applicable to calculating the area covered by support forces. It specifically addresses methods for calculating the coverage area of mobile communication base stations and emergency support warehouses. When multiple support force locations are given, the coverage area within the coverage area is a key indicator for optimizing support force deployment. The purpose of this invention is to provide a rapid method for calculating the coverage area of support forces, with short calculation time and high solution quality. The concept is as follows: Green's formula is used to calculate the coverage area of support forces within the coverage area. Two methods are used to describe the boundary curves of the intersection of the coverage area and the support force coverage area: the boundary curve segments of the support force coverage circle are described by angles, while the boundary curves on the coverage area are described by the corresponding proportions of the start and end points. A sorted intersection operation method is used to obtain all boundary curves on the circumference of the coverage area within the coverage area. A sorted union operation method is used to obtain all boundary curves of the protected portion on the boundary of the coverage area. Specific steps include:
[0028] Get the coordinates (x, y) of all vertices on the boundary of the protected region. vh ,y vh (h=1,2,…,n) V ), where h is the vertex number, and the total number of vertices is n. V And the center and radius r of the protection circle corresponding to each protection force within the protection area. ci The support circle corresponding to the support force is defined by the deployment location C of the support force. i (x ci ,y ci The circle is constructed with (i = 1, 2, ..., m) as the center and the radius of the protective force as the radius.
[0029] For any given protection circle, a first polar angle set corresponding to the protection circle is calculated based on the center and radius of the protection circle and the centers and radii of each protection circle in the remaining protection circle set. The first polar angle set includes the polar angles of each intersection point in the first intersection point set corresponding to the protection circle in the target polar coordinate system. The remaining protection circle set includes all protection circles excluding the protection circle. The target polar coordinate system is constructed with the center of the protection circle as the pole, and two intersecting protection circles have two intersection points. The first intersection point set corresponding to the protection circle includes the intersection points of the protection circle with each protection circle in the remaining protection circle set.
[0030] The second polar angle set of the protection circle is calculated based on the center and radius of the protection circle and the endpoint coordinates of all line segments on the boundary of the protection area. The starting point ratio and ending point ratio of each line segment when the protection circle intersects with all line segments on the boundary of the protection area are also calculated. The second polar angle set includes the starting point polar angle and ending point polar angle of the protection circle in the target polar coordinate system when the protection circle intersects with all line segments on the boundary of the protection area.
[0031] The starting and ending coordinates of all line segments on the boundary of the protected area are determined based on the starting and ending proportions of each line segment when all the protected circles intersect with all line segments on the boundary of the protected area.
[0032] The first polar angle set and the second polar angle set corresponding to the protection circle are merged and reduced to obtain the total polar angle set corresponding to the protection circle.
[0033] The coverage area of the protection force over the protection area is calculated based on the radius and center of all protection circles, the total set of polar angles corresponding to all protection circles, and the starting and ending coordinates of all line segments on the boundary of the protection area.
[0034] In practical applications, the step of calculating the first polar angle set corresponding to the protection circle based on the center and radius of the protection circle and the centers and radii of each protection circle in the remaining set of protection circles specifically includes:
[0035] For any one of the remaining guarantee circles in the set, it is taken as the target guarantee circle.
[0036] The intersection of the protection circle and the target protection circle is defined as including the start point and the end point; based on the center and radius of the protection circle and the center and radius of the target protection circle, the polar angles of the two intersection points of the protection circle and the target protection circle in the target polar coordinate system are calculated to obtain the start point polar angle and the end point polar angle of the intersection points of the protection circle and the target protection circle in the target polar coordinate system.
[0037] If the starting polar angle of the intersection point of the protection circle and the target protection circle in the target polar coordinate system is less than the ending polar angle of the intersection point of the protection circle and the target protection circle in the target polar coordinate system, then the polar angle of the starting point of the intersection point of the protection circle and the target protection circle in the target polar coordinate system is determined to be a Class 1 angle, and the polar angle of the ending point of the intersection point of the protection circle and the target protection circle in the target polar coordinate system is determined to be a Class 2 angle.
[0038] If the starting polar angle of the intersection point of the protection circle and the target protection circle in the target polar coordinate system is greater than the ending polar angle of the intersection point of the protection circle and the target protection circle in the target polar coordinate system, then the polar angle of the starting point of the intersection point of the protection circle and the target protection circle in the target polar coordinate system is determined to be a type 3 angle, and the polar angle of the ending point of the intersection point of the protection circle and the target protection circle in the target polar coordinate system is determined to be a type 4 angle.
[0039] Based on the category of the starting polar angle and the category of the ending polar angle of each intersection point in the first intersection point set corresponding to the protection circle in the target polar coordinate system, the polar angles of each intersection point in the first intersection point set corresponding to the protection circle in the target polar coordinate system are reduced to obtain the first polar angle set corresponding to the protection circle.
[0040] In practical applications, based on the center and radius of the protection circle and the center and radius of the target protection circle, the polar angles of the two intersection points of the protection circle and the target protection circle in the target polar coordinate system are calculated to obtain the starting polar angle and ending polar angle of the intersection points of the protection circle and the target protection circle in the target polar coordinate system. Specifically, this includes:
[0041] Let s be the initial number of pairs of intersections (start and end points) between the two protective circles. ci =0, calculate the intersection point of the two protective circles i1, i2, where i1, i2 = 1, 2, ..., m, i1 ≠ i2, and the protective circles i1, i2 have two disjoint intersection points, which satisfies formula (2)-formula (3). According to formula (3) and formula (4), determine the center of the protective circle i1. Let be the polar angle between the two intersection points of the pole in the polar coordinate system, and call it the starting polar angle on the protective circle i1. and the endpoint polar angle This represents the radius of the guaranteed circle i1. This represents the radius of the guaranteed circle i2.
[0042] The distance between the centers is
[0043]
[0044]
[0045]
[0046]
[0047]
[0048] in, Let i represent the logarithm of the starting and ending points corresponding to the guaranteed circle i1. Save to In this context, these two variables represent the polar angles of the starting and ending points of the boundary curve segment of the communication coverage area on the protection circle i1 after the protection circle i1 intersects with the protection circle i2. After obtaining the polar angles of the intersection points of the protection circle i1 with other circles, the values are calculated according to formula (5). and Update.
[0049] Based on formulas (6) and (7), determine the center of circle i2. Let i1 and i2 be the polar angles at the two intersection points of the protective circles i1 and i2 in the polar coordinate system with the pole as the pole, and call them the starting polar angles on the protective circle i2. and the endpoint polar angle
[0050]
[0051]
[0052]
[0053] in, Let i represent the logarithm of the starting and ending points corresponding to the guaranteed circle i2. Save to In this context, these two variables represent the polar angles of the starting and ending points of the arc segment of the communication coverage area boundary curve on the protection circle i2 after the protection circle i1 intersects with the protection circle i2. After obtaining the polar angles of the intersection points of the protection circle i2 with other circles, the polar angles are calculated according to formula (8). and Update.
[0054] Repeat this step until the intersection of any two safeguard circles has been calculated. If the intersection of the two safeguard circles does not satisfy formula (2), there is no intersection that needs to be saved.
[0055] In practical applications, the step of reducing the polar angles of each intersection point in the first intersection point set corresponding to the protection circle in the target polar coordinate system based on the category of the starting polar angle and the category of the ending polar angle of each intersection point in the target polar coordinate system, to obtain the first polar angle set corresponding to the protection circle, specifically includes:
[0056] The starting polar angle with the largest value among the Class 1 angles corresponding to the first intersection point set corresponding to the protection circle in the target polar coordinate system is determined as the first circle threshold.
[0057] The endpoint polar angle with the smallest value among the two types of angles corresponding to the first intersection point set corresponding to the protection circle in the target polar coordinate system is determined as the second circle threshold.
[0058] The starting polar angle with the largest value among the three types of angles corresponding to the first intersection point set corresponding to the protection circle in the target polar coordinate system is determined as the third circle threshold.
[0059] The endpoint polar angle with the smallest value among the four types of angles corresponding to the first intersection point set corresponding to the protection circle in the target polar coordinate system is determined as the fourth circle threshold.
[0060] The first set of polar angles corresponding to the protection circle is obtained by deleting the polar angles of each intersection point in the first set of intersection points corresponding to the protection circle that are less than the first circle threshold, greater than the second circle threshold, or satisfy the first condition in the polar coordinate system; the first condition is that the polar angle is greater than the fourth circle threshold and less than the third circle threshold.
[0061] In practical applications, the first polar angle set corresponding to the protection circle is obtained by deleting polar angles from the first set of intersection points corresponding to the protection circle that are less than the first circle threshold, greater than the second circle threshold, or satisfy the first condition in the polar coordinate system. Specifically, this includes:
[0062] Angles that are less than the first circle threshold, greater than the second circle threshold, or meet the first condition are deleted. Among the remaining angles, the first polar angle set is obtained by taking the angles corresponding to the first type of angle or the third type of angle as the starting point and the angles corresponding to the end point of the remaining angles.
[0063] In practical applications, the calculation of the second polar angle set of the protection circle and the starting and ending proportions of each line segment when the protection circle intersects with all line segments on the boundary of the protection area, based on the center and radius of the protection circle and the endpoint coordinates of all line segments on the boundary of the protection area, specifically includes:
[0064] Take any line segment on the boundary of the protected area as the target line segment.
[0065] Calculate the starting point ratio and ending point ratio of the target line segment when the protective circle intersects the target line segment based on the coordinates of the center of the protective circle and the endpoints of the target line segment.
[0066] Based on the center and radius of the protective circle and the coordinates of the endpoints of the target line segment, calculate the starting polar angle and ending polar angle of the protective circle in the target polar coordinate system when the protective circle intersects the target line segment.
[0067] If the starting polar angle of the protective circle in the target polar coordinate system when it intersects with the target line segment is less than the ending polar angle of the protective circle in the target polar coordinate system when it intersects with the target line segment, then the starting polar angle of the protective circle in the target polar coordinate system when it intersects with the target line segment is determined to be a Class 1 angle, and the ending polar angle of the protective circle in the target polar coordinate system when it intersects with the target line segment is determined to be a Class 2 angle.
[0068] If the starting polar angle of the protective circle in the target polar coordinate system when it intersects with the target line segment is greater than the ending polar angle of the protective circle in the target polar coordinate system when it intersects with the target line segment, then the starting polar angle of the protective circle in the target polar coordinate system when it intersects with the target line segment is determined to be a type 3 angle, and the ending polar angle of the protective circle in the target polar coordinate system when it intersects with the target line segment is determined to be a type 4 angle.
[0069] The second set of polar angles corresponding to the protection circle is determined based on the category of the starting polar angle and the category of the ending polar angle of the protection circle in the target polar coordinate system when the protection circle intersects with all line segments on the boundary of the protection area.
[0070] In practical applications, the starting point ratio and ending point ratio of the target line segment when the protective circle intersects with the target line segment are calculated based on the center and endpoint coordinates of the protective circle and the target line segment. Specifically, the starting point polar angle and ending point polar angle of the protective circle in the target polar coordinate system when the protective circle intersects with the target line segment are calculated based on the center and radius of the protective circle and the endpoint coordinates of the target line segment.
[0071] Step 1.4: Calculate the intersection of the protection circle and the boundary of the protection area.
[0072] Let s be the initial number of pairs (start and end points) of intersections between the protection circle and the boundary segment of the protection area. ai =0, the initial number of intersection pairs (start and end points) between the protection circle and the boundary segment of the protected area on the boundary segment is s. vh =0(h=1,2,…,n) V Calculate the intersection point of the circumference of the protection circle i (i = 1, 2, ..., m) and the boundary segment PQ of the protection area, with the center of the protection circle i (x... ci ,y ci Let P(x) be the polar angle in the polar coordinate system of the pole, where P(x) is the polar angle of the pole. vh ,y vh ),Q(x vh+1 ,y vh+1 h = 1, ..., n Vh represents the vertex index, n V This represents the total number of vertices. The direction of line segment PQ is the same as the positive direction of the boundary curve of the protected area. ci This represents the radius of the guaranteed circle i.
[0073] For ease of description, remember
[0074] a=(x vh+1 -x vh ) 2 +(y vh+1 -y vh ) 2
[0075] b = 2(x) vh+1 -x vh (x) vh -x ci )+2(y vh+1 -y vh )(y vh -y ci )
[0076]
[0077]
[0078] Δ=b 2 -4ac1
[0079] (1) When Δ>0, -2a<b<0, a+b+c1≥0, min{c1,c2}≥0, the circle and line segment PQ have two intersection points D1 and D2. Intersection point D1 is closer to the endpoint P of the line segment, and intersection point D2 is closer to the endpoint Q of the line segment. The polar angles γ and θ corresponding to the starting and ending points on circle i are determined by formulas (9) and (10) respectively, and are saved to [the appropriate database]. The intersection point can also be viewed as a point on the boundary line segment with a fixed ratio; the starting and ending points are stored proportionally respectively. middle.
[0080]
[0081]
[0082]
[0083]
[0084] (2) When c1 < 0 < c2, the circle and the line segment have only one intersection point D2. This intersection point is the starting point. The polar angle γ corresponding to the starting point and the ending point on circle i is determined by formula (9).
[0085]
[0086]
[0087] (3) When c2 < 0 < c1, the circle and the line segment have only one intersection point D1. This intersection point is the endpoint. The polar angle θ between the starting point and the endpoint on circle i is determined by formula (10).
[0088]
[0089]
[0090] In practical applications, determining the start and end coordinates of all line segments on the boundary of the protected area based on the start and end ratios of each line segment when all protected circles intersect with all line segments on the boundary of the protected area specifically includes:
[0091] The sorting sequence is obtained by sorting the starting point ratio and ending point ratio of each line segment when all the protection circles intersect with all line segments on the boundary of the protection area in ascending order.
[0092] If two adjacent values in the sorted sequence are both starting proportions, then the larger starting proportion of the two adjacent values is deleted to obtain the sorted sequence after deletion. If two adjacent values in the sorted sequence are both ending proportions, then the smaller ending proportion of the two adjacent values is deleted to obtain the sorted sequence after deletion.
[0093] The starting and ending coordinates of all line segments on the boundary of the protected area are obtained based on the sorted sequence after deletion.
[0094] In practical applications, the starting and ending coordinates of all line segments on the boundary of the protected area are obtained based on the sorted sequence after deletion, specifically including:
[0095] The intersection information of each edge is summarized together, and the starting coordinates (x, y) of each edge are also included. a2j-1 y a2j-1 The coordinates of the endpoint (x) on each edge are determined by formula (17); a2j y a2j ) is determined by formula (18).
[0096] (x a2j-1 ,y a2j-1 )=(r shk (x vh+1 -x vh )+x vh ,r shk (y vh+1 -y vh )+y vh (17)
[0097] (xa2j ,y a2j )=(r thk (x vh+1 -x vh )+x vh ,r thk (y vh+1 -y vh )+y vh (18) where h = 1, 2, ..., n V j = 1, 2, ..., n A ; Following the positive direction of the boundary curve, h represents the starting point V. h With the endpoint V h+1 The corresponding edge number, k represents the number of the start-end pair on this edge, for a given h, k = 1, 2, ..., s vh j represents the unique number of all start-end pairs on the boundary line, which increases with the increase of h and k.
[0098] In practical applications, determining the second polar angle set corresponding to the protection circle based on the category of the starting polar angle and the category of the ending polar angle of the protection circle in the target polar coordinate system when the protection circle intersects with all line segments on the boundary of the protection area specifically includes:
[0099] When the protection circle intersects with the line segment on the boundary of the protection area, the starting polar angle and ending polar angle of the protection circle in the target polar coordinate system are sorted in ascending order to obtain a sorting sequence. The starting polar angle and ending polar angle of the protection circle that are less than 0 in the target polar coordinate system are deleted to obtain a set of deleted polar angles.
[0100] Iterate through each polar angle in the set of deleted polar angles;
[0101] If the polar angle currently being traversed is a type 1 angle or a type 3 angle, then the polar angle with the smallest value among the polar angle currently being traversed and the type 2 or type 4 angle that is larger than the polar angle currently being traversed is determined as the starting polar angle and the ending polar angle.
[0102] If the polar angle currently traversed is the largest polar angle in the deleted polar angle set, and the currently traversed polar angle is a Class 3 angle, then the polar angle currently traversed and the smallest polar angle among all Class 2 and Class 4 angles in the deleted polar angle set are called the starting polar angle and the ending polar angle. That is, the starting polar angle and ending polar angle corresponding to the intersection points of all line segments on the boundary of the protection circle and the protection area in the target polar coordinate system are arranged in ascending order, and polar angles with values less than 0 are deleted. A Class 1 angle (Class 3 angle) and a Class 2 angle (Class 4 angle) with a value larger than that angle form a new starting angle-ending angle pair. If the maximum value of all angles is a Class 3 angle, that angle is the starting angle, and its corresponding ending angle is the smallest polar angle among all Class 2 and Class 4 angles.
[0103] In practical applications, the first polar angle set and the second polar angle set corresponding to the protection circle are merged and subtracted to obtain the total polar angle set corresponding to the protection circle, specifically including:
[0104] Let k denote the index of the start-end pair on the guaranteed circle, and let α represent the start-end polar angle and the end-end polar angle in the second polar angle set. aik ,β aik (k = 1, 2, ..., s) ai The starting and ending polar angles in the first polar angle set are represented by α. cik ,β cik (k = 1, 2, ..., s) ci ) indicates that when α aik <β aik At that time, α aik For type 1 angles, β aik It is a type 2 angle. α aik >β aik At that time, α aik It is a type 3 angle, β aik It is a type 4 angle. α cik <β cik At that time, α cik For type 1 angles, β cik It is a type 2 angle. α cik >β cik At that time, α cik For three types of angles, β cikThere are four types of angles. All polar angles in the first and second polar angle sets corresponding to the protection circle are arranged in ascending order. Angles smaller than the first circle threshold (the largest value among all type 1 angles) are deleted; angles larger than the second circle threshold (the smallest value among all type 2 angles) are deleted; and angles between the fourth circle threshold (the smallest value among all type 4 angles) and the third circle threshold (the largest value among all type 3 angles) are deleted. Of the remaining angles, angles corresponding to either type 1 or type 3 angles are sequentially assigned as starting points and ending points (if the first data point is a type 4 angle, it can be used as the last starting point and ending point angle), denoted as α. il ,β il (l=1,2,…,s i ), where l represents the number of the start-end pair after removing redundant angles, s i This represents the number of start-end pairs after removing redundant angles. If a guaranteed circle's center C... i Within the protected area, and if it intersects with no other circles or the boundary of the protected area, then α il =0,β il =2π.
[0105] In practical applications, the coverage area of the protection force over the protection area is calculated based on the radius and center of all protection circles, the total set of polar angles corresponding to all protection circles, and the start and end coordinates of all line segments on the boundary of the protection area. Specifically, this includes:
[0106] Step 2.1: Input the boundary curve of the coverage area when deploying support forces;
[0107] Assume there are m support forces in total, and the coordinates of the base station (or emergency support warehouse) deployment location are C. i (x ci ,y ci (i = 1, 2, ..., m), effective radius of action r ci , with C i With the origin as the origin and the positive x-axis as the polar axis, the effective working circle (or emergency support warehouse) of this base station (or emergency support warehouse) within the support area has polar angles corresponding to the starting and ending points of the positive boundary curve, respectively. ik β ik (k = 1, 2, ..., s) i This refers to the total set of polar angles corresponding to all protected circles; input the starting coordinates (x, y) of the positive boundary segment of the protected area. a2j-1 ,y a2j-1 ) and endpoint coordinates (x a2j ,y a2j The starting and ending points are numbered j = 1, 2, ..., n AThe starting and ending points are both on the same straight line segment and are arranged along the positive boundary curve of the protected area.
[0108] Step 2.2: The area A of the support force covering the support zone is...
[0109]
[0110] This invention also provides a system for calculating coverage area in accordance with the above method, comprising:
[0111] The acquisition module is used to acquire the coordinates of all vertices on the boundary of the protection area and the center and radius of the protection circle corresponding to each protection force in the protection area; the protection circle corresponding to each protection force is constructed with the deployment position of the protection force as the center and the radius of action of the protection force as the radius.
[0112] The first polar angle set determination module is used to calculate, for any given protection circle, the first polar angle set corresponding to the protection circle based on the center and radius of the protection circle and the centers and radii of each protection circle in the remaining protection circle set; the first polar angle set includes the polar angles of each intersection point in the first intersection point set corresponding to the protection circle in the target polar coordinate system; the remaining protection circle set includes all protection circles excluding the protection circle; the target polar coordinate system is constructed with the center of the protection circle as the pole, and two intersecting protection circles have two intersection points; the first intersection point set corresponding to the protection circle includes the intersection points of the protection circle with each protection circle in the remaining protection circle set.
[0113] The second polar angle set and the proportion determination module are used to calculate the second polar angle set of the protection circle and the starting point proportion and ending point proportion of each line segment when the protection circle intersects with all line segments on the boundary of the protection area, based on the center and radius of the protection circle and the endpoint coordinates of all line segments on the boundary of the protection area; the second polar angle set includes the starting point polar angle and ending point polar angle of the protection circle in the target polar coordinate system when the protection circle intersects with all line segments on the boundary of the protection area respectively.
[0114] The start and end coordinate determination module is used to determine the start and end coordinates of all line segments on the boundary of the protection area based on the start and end ratios of each line segment when all protection circles intersect with all line segments on the boundary of the protection area.
[0115] The deletion module is used to merge and delete the first polar angle set and the second polar angle set corresponding to the protection circle to obtain the total polar angle set corresponding to the protection circle.
[0116] The coverage area calculation module is used to calculate the coverage area of the protection force covering the protection area based on the radius and center of all protection circles, the total set of polar angles corresponding to all protection circles, and the starting and ending coordinates of all line segments on the boundary of the protection area.
[0117] This invention also provides an electronic device, comprising:
[0118] A memory and a processor, the memory for storing a computer program, the processor for running the computer program to cause the electronic device to perform the coverage area calculation method according to the above embodiments.
[0119] This invention also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the method for calculating the coverage area as described in the above embodiments.
[0120] like Figure 1 As shown in the figure, this embodiment of the invention illustrates the above method using a specific example of eight emergency support warehouses protecting a rectangular area (100 km long and 50 km wide). In this embodiment, each warehouse has a protection radius of 20 km, and their deployment locations are shown in Table 1. The vertex coordinates of the rectangular area where the warehouses are located are (0,0), (100,0), (100,50), and (0,50). Calculate the protection area of these eight warehouses.
[0121] Table 1. Emergency Support Warehouse Deployment Information (Unit: km)
[0122]
[0123]
[0124] Based on the warehouse location information given in Table 1, the coverage area is 5000 km². 2 .
[0125] Step 1: Method for determining the start and end points of the boundary segment of the coverage area when deploying support forces:
[0126] Step 1.1: Input the vertices (x, y) of the polygon protected region. vh ,y vh (h=1,2,…,n) V ), n V =4, the effective radius r of the emergency support warehouse ci and the deployment location C i (x ci ,y ci (i = 1, 2, ..., m), m = 8. The initial number of pairs of intersections between the guarantee circles is s. ci=0, the initial number of pairs of intersections between the protection circle and the boundary segment of the protection area is s. ai =0, the initial number of intersection points between the protection circle and the boundary segment of the protected area on the boundary segment is s. vh =0(h=1,2,…,n) V ).
[0127] Step 1.2: Calculate the polar angles of the intersection points of the two protective circles i1 and i2 in a polar coordinate system with the centers of the two circles as poles, where i1, i2 = 1, 2, ..., m, i1 ≠ i2, and the circles i1 and i2 have two distinct intersection points (satisfying formulas (2)-(3)). The polar angles corresponding to the starting and ending points on circle i1. The polar angles corresponding to the starting and ending points on circle i2 are determined by formulas (3) and (4) respectively. The polar angle and the number of intersection points are determined by formulas (6) and (7) respectively, and the polar angle and the number of intersection points are saved according to formulas (5) and (8). Step 1.2 is repeated until the intersection points of any two protection circles are calculated. If formula (2) is not satisfied, there are no intersection points of the two protection circles that need to be saved. According to step 1.2, the intersection points on the circumference of the protection circles of the emergency warehouse are shown in Table 2.
[0128] Table 2. Polar angles of the intersection points of the circumferences of the emergency warehouse's protective circle (unit: radians)
[0129]
[0130]
[0131] *When the starting value is greater than the ending value, the starting value is equivalent to -2π.
[0132] Step 1.3: Eliminate redundant information about the intersection points of different protection circles on the circumference of the protection circle.
[0133] Define α cik ,β cik (k = 1, 2, ..., s) ci Category: α cik <β cik At that time, α cik For type 1 angles, β cik It is a type 2 angle. α cik >β cik At that time, α cik For three types of angles, β cik There are four types of angles.
[0134] α cik ,β cik (k = 1, 2, ..., s) ci Arrange them in ascending order.
[0135] Delete angles smaller than the first circle threshold (the largest among class 1 angles), delete angles larger than the second circle threshold (the smallest among class 2 angles), and delete angles between the fourth circle threshold (the smallest among class 4 angles) and the third circle threshold (the largest among class 3 angles).
[0136] Of the retained angles, those starting with either type 1 or type 3 angles are designated as the starting angles, and those ending with subsequent angles. For convenience, the angles after eliminating redundant intersection information on the circumference of the guaranteed circle are still denoted as α. cik ,β cik (k = 1, 2, ..., s) ci ).
[0137] Following step 1.3, redundant information about the intersection points of different protection circles on the circumference of the protection circle is eliminated. For example, the sorted results of the intersection points of protection circle 1 are shown in Table 2.1.
[0138] Table 2.1 Results of sorting the intersection points of circle 1
[0139] category 1 1 1 1 4 4 3 2 2 3 2 2
[0140] To delete angles less than Class 1, angles 0.0019, 0.1657, and 0.1795 should be deleted; to delete angles greater than Class 2, angles 5.3012, 5.4143, 5.9338, and 6.0187 should be deleted; to delete angles between Class 4 and Class 3, angles 3.6949, 4.3675, 4.5776, and 5.3012 should be deleted. After deleting redundant information, only angles 2.0271 and 2.8224 remain on circle 1. The results are shown in Table 3.
[0141] Table 3. Polar angles of the intersection points of the support circles of the emergency warehouse after eliminating redundant information (unit: radians)
[0142] 1 (2.0271,2.8224) 2 (0.1143,2.8771) 3 (5.3879,1.9883) 4 (3.4764,5.7404) 5 (4.5102,5.0665) 6 (1.2259,2.5033) (4.5424,4.8508) 7 (4.3477,6.1684) 8 (1.4008,5.6449)
[0143] Step 1.4: Calculate the intersection of the protection circle and the boundary of the protection area.
[0144] Calculate the intersection point of the circumference of the protection circle i (i = 1, 2, ..., m) and the boundary segment PQ of the protection area at the center (x) of the protection circle. ci ,y ci Let P(x) be the polar angle in the polar coordinate system of the pole, where P(x) is the polar angle of the pole. vh ,y vh ),Q(x vh+1 ,y vh+1 h = 1, ..., n V The direction of line segment PQ is the same as the positive direction of the boundary curve of the protected area.
[0145] (1) When the circle and line segment PQ have two intersection points D1 and D2, intersection point D1 is closer to the endpoint P of the line segment, and intersection point D2 is closer to the endpoint Q of the line segment. The polar angles γ and θ corresponding to the starting and ending points on circle i are determined by formulas (9) and (10), respectively. s ai , and s vh Determined according to formulas (11) and (12).
[0146] (2) When the circle and the line segment have only one intersection point D2, the polar angle γ between the starting point and the ending point on circle i is determined by formula (9). s ai , and s vh Determined according to formulas (13) and (14).
[0147]
[0148]
[0149] (3) When the circle and the line segment have only one intersection point D1, the polar angle θ between the starting point and the ending point on circle i is determined by formula (10). s ai , and s vh Determined according to formulas (15) and (16).
[0150] Step 1.5: Eliminate redundant intersection information of the protection circle and boundary line segment on the polygon protection area.
[0151] r shk ,r thk Sort them in ascending order. If two adjacent data points in the sort are both starting points on the boundary line segment r, then... shk If the larger starting data is deleted, then delete the data from the end point on the boundary line segment. thk If the smaller endpoint data is not found, then the smaller endpoint data will be deleted. The proportions of the start and end points in the final result are denoted as r. shk r thk .
[0152] The intersection information of each edge is summarized together, and the starting point coordinates are determined by formula (17); the ending point coordinates are determined by formula (18).
[0153] Following steps 1.4-1.5, the intersection points on the boundary line segments of the protected area are obtained. The results after eliminating redundant information are shown in Table 4.
[0154] Table 4 shows the intersection points of the protection circle and the protection area boundary on the protection area boundary segment after eliminating redundant information.
[0155] Table 4. Intersections of boundary segments (start and end points are shown proportionally).
[0156]
[0157] Step 1.6: Eliminate redundant intersection information between the protection circle and the boundary line segment on the protection circle circumference.
[0158] Define α aik ,β aik (k = 1, 2, ..., s) ai Category: α aik <β aik At that time, α aik For type 1 angles, β aik It is a type 2 angle. α aik >β aik At that time, α aik For three types of angles, β aik There are four types of angles.
[0159] α aik ,β aik (k = 1, 2, ..., s) ai Arrange the angles in ascending order. Delete polar angles with values less than 0. Among the remaining angles, use either a Class 1 or Class 3 angle as the starting point and subsequent angles as the ending point (if the first angle is a Class 4 angle, it can be used as the last starting point to ending point). For convenience, the angle after eliminating redundant intersection information on the circumference of the circle is still denoted as α. aik β aik (k = 1, 2, ..., s) ai ).
[0160] Following step 1.4, the intersection point between the protection circle and the boundary of the protection area is obtained. After eliminating redundant information in step 1.6, the intersection points on the circumference are shown in Table 5.
[0161] Table 5. Intersection points of the boundary segment of the protected area on the protected circle and the protected circle after eliminating redundant information (unit: radians)
[0162] 1 (2.2838,0.8578) 2 (3.0212,0.1204) 3 (1.9731,5.2416) (1.0416,1.1685) 4 (5.8113,3.6134) 5 (5.1210,4.3038) 6 7 (5.1302,5.7736) (0.5096,4.2945) 8 (5.4699,1.8879)
[0163] Step 1.7: Eliminate redundant intersection information on the circumference of the protection circle.
[0164] Define α aik ,β aik (k = 1, 2, ..., s) ai ) and α cik ,β cik (k = 1, 2, ..., s) ci Category: α aik <β aik At that time, αaik For type 1 angles, β aik It is a type 2 angle. α aik >β aik At that time, α aik For three types of angles, β aik It is a type 4 angle. α cik <β cik At that time, α cik For type 1 angles, β cik It is a type 2 angle. α cik >β cik At that time, α cik For three types of angles, β cik There are four types of angles.
[0165] Arrange the above data in ascending order.
[0166] Delete angles that are less than the first circle threshold (the largest value among all Class 1 angles), delete angles that are greater than the second circle threshold (the smallest value among all Class 2 angles), and delete angles that are between the fourth circle threshold (the smallest value among all Class 4 angles) and the third circle threshold (the largest value among all Class 3 angles).
[0167] In the retained data, the angles corresponding to the starting point are successively defined as either Class 1 or Class 3 angles, and the angles corresponding to the ending point are successively defined as the following angles (if the first angle is a Class 4 angle, it can be used as the angle corresponding to the last starting point and ending point), denoted as α. ik ,β ik (k = 1, 2, ..., s) i If a guaranteed circle's center C... i Within the protected area, and if it intersects with no other circles or the boundary of the protected area, then α ik =0,β ik =2π.
[0168] Following step 1.7, redundant intersection information on the circumference of the protection circle is eliminated, and the results are shown in Table 6.
[0169] Table 6. Intersection Points on the Guarantee Circle After Redundancy Information Removal (Unit: radians)
[0170]
[0171] Step 2: Mathematical model for the coverage area of the emergency support warehouse during deployment:
[0172] Step 2.1: Input the boundary curve of the coverage area when deploying the emergency support warehouse.
[0173] Assume there are m emergency supply warehouses in total, with an effective radius of r. ci The emergency support warehouse is located at coordinates C. i (x ci,y ci ), i = 1, 2, ..., m. Let C i With the origin as the origin and the positive x-axis as the polar axis, the polar angles corresponding to the starting and ending points of the positive boundary curve of the effective function circle of this emergency support warehouse within the support area are α and α, respectively. ik β ik (k = 1, 2, ..., s) i ); Input the starting coordinates (x, y) of the positive boundary segment of the protected area. a2j-1 ,y a2j-1 ) and endpoint coordinates (x a2j ,y a2j (j=1,2,…,n) A The starting and ending points are both on the same straight line segment and are arranged along the positive boundary curve of the protected area;
[0174] Step 2.2: Calculate the area covered by all emergency support warehouses according to formula (19).
[0175] Following step 2, the area covered by the eight emergency warehouses was calculated to be A = 4559.7522 km². 2 .
[0176] As can be seen from the solution process in this example, when calculating the actual area of the emergency support warehouse within the support area, this invention, in order to simplify the problem, successively calculates the intersection points of the support circles on the support circle, the intersection points of the support circle and the boundary curve of the support area, and the intersection points of the support circle and the boundary of the support area on the boundary of the support area. When redundant information needs to be eliminated, a sorting method is used to quickly obtain the calculation results. For the problem of finding the intersection (or union) of n sets, the number of comparisons required by the traditional method is O(n^2). 2 The number of comparisons required by the sorting method is O(nlnn), which improves the computation speed of eliminating redundant information.
[0177] For cases where the protected area is circular or elliptical, a similar method can be used to quickly and accurately determine the area covered by the protected area within the protected region.
[0178] Advantages of this invention:
[0179] 1. It enables automatic and accurate calculation of the coverage area for issues such as mobile communication base station and warehouse site selection.
[0180] 2. A sorting method is used to eliminate redundant intersections on the circumference of the guarantee circle, enabling fast calculation of intersection sets of multiple sets.
[0181] 3. A sorting method is used to eliminate redundant intersections on the boundary segments of the protected area, enabling fast calculation of the union of multiple sets.
[0182] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For the systems disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the descriptions are relatively simple; relevant parts can be referred to the method section.
[0183] This document uses specific examples to illustrate the principles and implementation methods of the present invention. The descriptions of the above embodiments are only for the purpose of helping to understand the method and core ideas of the present invention. Furthermore, those skilled in the art will recognize that, based on the ideas of the present invention, there will be changes in the specific implementation methods and application scope. Therefore, the content of this specification should not be construed as a limitation of the present invention.
Claims
1. A method for calculating coverage area, characterized in that, include: Obtain the coordinates of all vertices on the boundary of the protection area, as well as the center and radius of the protection circle corresponding to each protection force within the protection area; the protection circle corresponding to each protection force is constructed with the deployment position of the protection force as the center and the radius of action of the protection force as the radius; For any given protection circle, a first polar angle set corresponding to the protection circle is calculated based on the center and radius of the protection circle and the centers and radii of the remaining protection circles. The first polar angle set includes the polar angles of each intersection point in the first intersection point set corresponding to the protection circle in the target polar coordinate system. The remaining protection circle set includes all protection circles excluding the protection circle. The target polar coordinate system is constructed with the center of the protection circle as its pole, and two intersecting protection circles have two intersection points. The first intersection point set corresponding to the protection circle includes the intersection points of the protection circle with each of the remaining protection circles. Based on the center and radius of the protection circle and the centers and radii of each protection circle in the remaining set of protection circles, calculate the first polar angle set corresponding to the protection circle, specifically including: For any one of the remaining protection circles in the set, it is taken as the target protection circle; The intersection of the protection circle and the target protection circle is defined to include a start point and an end point; based on the center and radius of the protection circle and the center and radius of the target protection circle, the polar angles of the two intersection points of the protection circle and the target protection circle in the target polar coordinate system are calculated to obtain the start point polar angle and the end point polar angle of the intersection points of the protection circle and the target protection circle in the target polar coordinate system. If the starting polar angle of the intersection of the protection circle and the target protection circle in the target polar coordinate system is less than the ending polar angle of the intersection of the protection circle and the target protection circle in the target polar coordinate system, then the starting polar angle of the intersection of the protection circle and the target protection circle in the target polar coordinate system is determined to be a Class 1 angle, and the ending polar angle of the intersection of the protection circle and the target protection circle in the target polar coordinate system is determined to be a Class 2 angle. If the starting polar angle of the intersection of the protection circle and the target protection circle in the target polar coordinate system is greater than the ending polar angle of the intersection of the protection circle and the target protection circle in the target polar coordinate system, then the starting polar angle of the intersection of the protection circle and the target protection circle in the target polar coordinate system is determined to be a type 3 angle, and the ending polar angle of the intersection of the protection circle and the target protection circle in the target polar coordinate system is determined to be a type 4 angle. Based on the category of the starting polar angle and the category of the ending polar angle of each intersection point in the first intersection point set corresponding to the protection circle in the target polar coordinate system, the polar angles of each intersection point in the first intersection point set corresponding to the protection circle in the target polar coordinate system are reduced to obtain the first polar angle set corresponding to the protection circle. The second polar angle set of the protection circle is calculated based on the center and radius of the protection circle and the endpoint coordinates of all line segments on the boundary of the protection area. The starting point ratio and ending point ratio of each line segment when the protection circle intersects with all line segments on the boundary of the protection area are also calculated. The second polar angle set includes the starting point polar angle and ending point polar angle of the protection circle in the target polar coordinate system when the protection circle intersects with all line segments on the boundary of the protection area. The starting coordinates and ending coordinates of all line segments on the boundary of the protected area are determined based on the starting and ending proportions of each line segment when all the protected circles intersect with all line segments on the boundary of the protected area. The first polar angle set and the second polar angle set corresponding to the protection circle are merged and subtracted to obtain the total polar angle set corresponding to the protection circle. The coverage area of the protection force over the protection area is calculated based on the radius and center of all protection circles, the total set of polar angles corresponding to all protection circles, and the starting and ending coordinates of all line segments on the boundary of the protection area. Specifically, the coverage area is calculated according to the formula... The calculation is performed, where A represents the coverage area of the support force. Indicates the total number of line segments. This represents the coordinates of the endpoint of the j-th line segment. This represents the coordinates of the starting point of the j-th line segment. Indicates the total number of support personnel. This represents the radius of the protection circle corresponding to the i-th protection force. This represents the number of start-end pairs in the total set of polar angles corresponding to the i-th support force's support circle. This represents the k-th starting polar angle in the total set of polar angles corresponding to the polar circle of the i-th support force. This represents the k-th endpoint polar angle in the total set of polar angles corresponding to the polar circles of the i-th support force. This represents the coordinates of the center of the protection circle corresponding to the i-th protection force.
2. The method for calculating coverage area according to claim 1, characterized in that, The calculation of the second polar angle set of the protection circle and the starting and ending proportions of each line segment when the protection circle intersects with all line segments on the boundary of the protection area, based on the center and radius of the protection circle and the endpoint coordinates of all line segments on the boundary of the protection area, specifically includes: Take any line segment on the positive boundary of the protected area as the target line segment; Calculate the starting point ratio and ending point ratio of the target line segment when the protective circle intersects with the target line segment based on the center of the protective circle and the endpoint coordinates of the target line segment; Based on the center and radius of the protective circle and the coordinates of the endpoints of the target line segment, calculate the starting polar angle and ending polar angle of the protective circle in the target polar coordinate system when the protective circle intersects with the target line segment; If the starting polar angle of the protective circle in the target polar coordinate system when it intersects with the target line segment is less than the ending polar angle of the protective circle in the target polar coordinate system when it intersects with the target line segment, then the starting polar angle of the protective circle in the target polar coordinate system when it intersects with the target line segment is determined to be a Class 1 angle, and the ending polar angle of the protective circle in the target polar coordinate system when it intersects with the target line segment is determined to be a Class 2 angle. If the starting polar angle of the protective circle in the target polar coordinate system when it intersects with the target line segment is greater than the ending polar angle of the protective circle in the target polar coordinate system when it intersects with the target line segment, then the starting polar angle of the protective circle in the target polar coordinate system when it intersects with the target line segment is determined to be a type 3 angle, and the ending polar angle of the protective circle in the target polar coordinate system when it intersects with the target line segment is determined to be a type 4 angle. The second set of polar angles corresponding to the protection circle is determined based on the category of the starting polar angle and the category of the ending polar angle of the protection circle in the target polar coordinate system when the protection circle intersects with all line segments on the boundary of the protection area.
3. The method for calculating the coverage area according to claim 1, characterized in that, The step of reducing the polar angles of each intersection point in the first intersection point set corresponding to the protection circle in the target polar coordinate system according to the category of the starting polar angle and the category of the ending polar angle of each intersection point in the target polar coordinate system, to obtain the first polar angle set corresponding to the protection circle, specifically includes: The first circle threshold is determined by identifying the angle with the largest value among the starting polar angles of each intersection point in the first intersection point set corresponding to the protection circle in the target polar coordinate system. The second circle threshold is determined by identifying the type 2 angle with the smallest value among the endpoint polar angles corresponding to the first intersection point set of the protection circle in the target polar coordinate system. The three types of angles with the largest numerical value among the starting polar angles of each intersection point in the first intersection point set corresponding to the protection circle in the target polar coordinate system are determined as the third circle threshold. The four types of angles with the smallest values among the endpoint polar angles corresponding to the first intersection points in the target polar coordinate system corresponding to the protection circle are determined as the fourth circle threshold. The first set of polar angles corresponding to the protection circle is obtained by deleting the polar angles of each intersection point in the first set of intersection points corresponding to the protection circle that are less than the first circle threshold, greater than the second circle threshold, or satisfy the first condition in the polar coordinate system; the first condition is that the polar angle is greater than the fourth circle threshold and less than the third circle threshold.
4. The method for calculating the coverage area according to claim 1, characterized in that, The step of determining the start and end coordinates of all line segments on the boundary of the protected area based on the start and end proportions of each line segment when all protected circles intersect with all line segments on the boundary of the protected area specifically includes: When all the protection circles intersect with all line segments on the boundary of the protection area, the starting point ratio and ending point ratio of each line segment are sorted in ascending order to obtain the sorting sequence; If two adjacent values in the sorting sequence are both starting proportions, then the larger starting proportion of the two adjacent values is deleted to obtain the sorting sequence after deletion. If two adjacent values in the sorting sequence are both ending proportions, then the smaller ending proportion of the two adjacent values is deleted to obtain the sorting sequence after deletion. The starting and ending coordinates of all line segments on the boundary of the protected area are obtained based on the sorted sequence after deletion.
5. The method for calculating the coverage area according to claim 2, characterized in that, The step of determining the second polar angle set corresponding to the protection circle based on the category of the starting polar angle and the category of the ending polar angle of the protection circle in the target polar coordinate system when the protection circle intersects with all line segments on the boundary of the protection area specifically includes: When the protection circle intersects with all line segments on the boundary of the protection area, delete the starting polar angle and ending polar angle of the protection circle that are less than 0 in the target polar coordinate system to obtain the set of deleted polar angles; Iterate through each polar angle in the set of deleted polar angles; If the polar angle currently being traversed is a type 1 angle or a type 3 angle, then the polar angle currently being traversed and the smallest polar angle among type 2 or type 4 angles that are larger than the polar angle currently being traversed are determined to be the starting polar angle and the ending polar angle, respectively, forming a starting polar angle-ending polar angle pair; If the polar angle currently being traversed is the polar angle with the largest value in the set of polar angles after deletion, and the polar angle currently being traversed is a type 3 angle, then the polar angle currently being traversed and the polar angle with the smallest value among all type 2 and type 4 angles in the set of polar angles after deletion are determined to be the starting polar angle and the ending polar angle, respectively, forming a starting polar angle-ending polar angle pair.
6. A system for calculating coverage area, characterized in that, include: The acquisition module is used to acquire the coordinates of all vertices on the boundary of the protection area and the center and radius of the protection circle corresponding to each protection force within the protection area; the protection circle corresponding to each protection force is constructed with the deployment position of the protection force as the center and the radius of action of the protection force as the radius; The first polar angle set determination module is used to calculate, for any given protection circle, the first polar angle set corresponding to the protection circle based on the center and radius of the protection circle and the centers and radii of each protection circle in the remaining protection circle set. The first polar angle set includes the polar angles of each intersection point in the first intersection point set corresponding to the protection circle in the target polar coordinate system. The remaining protection circle set includes all protection circles excluding the protection circle. The target polar coordinate system is constructed with the center of the protection circle as its pole, and two intersecting protection circles have two intersection points. The first intersection point set corresponding to the protection circle includes the intersection points of the protection circle with each protection circle in the remaining protection circle set. Calculating the first polar angle set corresponding to the protection circle based on the center and radius of the protection circle and the centers and radii of each protection circle in the remaining protection circle set specifically includes: For any one of the remaining protection circles in the set, it is taken as the target protection circle; The intersection of the protection circle and the target protection circle is defined to include a start point and an end point; based on the center and radius of the protection circle and the center and radius of the target protection circle, the polar angles of the two intersection points of the protection circle and the target protection circle in the target polar coordinate system are calculated to obtain the start point polar angle and the end point polar angle of the intersection points of the protection circle and the target protection circle in the target polar coordinate system. If the starting polar angle of the intersection of the protection circle and the target protection circle in the target polar coordinate system is less than the ending polar angle of the intersection of the protection circle and the target protection circle in the target polar coordinate system, then the starting polar angle of the intersection of the protection circle and the target protection circle in the target polar coordinate system is determined to be a Class 1 angle, and the ending polar angle of the intersection of the protection circle and the target protection circle in the target polar coordinate system is determined to be a Class 2 angle. If the starting polar angle of the intersection of the protection circle and the target protection circle in the target polar coordinate system is greater than the ending polar angle of the intersection of the protection circle and the target protection circle in the target polar coordinate system, then the starting polar angle of the intersection of the protection circle and the target protection circle in the target polar coordinate system is determined to be a type 3 angle, and the ending polar angle of the intersection of the protection circle and the target protection circle in the target polar coordinate system is determined to be a type 4 angle. Based on the category of the starting polar angle and the category of the ending polar angle of each intersection point in the first intersection point set corresponding to the protection circle in the target polar coordinate system, the polar angles of each intersection point in the first intersection point set corresponding to the protection circle in the target polar coordinate system are reduced to obtain the first polar angle set corresponding to the protection circle. The second polar angle set and the proportion determination module are used to calculate the second polar angle set of the protection circle and the starting and ending proportions of each line segment when the protection circle intersects with all line segments on the boundary of the protection area, based on the center and radius of the protection circle and the endpoint coordinates of all line segments on the boundary of the protection area; the second polar angle set includes the starting and ending polar angles of the protection circle in the target polar coordinate system when the protection circle intersects with all line segments on the boundary of the protection area respectively; The start and end coordinate determination module is used to determine the start and end coordinates of all line segments on the boundary of the protection area based on the start and end ratios of each line segment when all protection circles intersect with all line segments on the boundary of the protection area. The deletion module is used to merge and delete the first polar angle set and the second polar angle set corresponding to the protection circle to obtain the total polar angle set corresponding to the protection circle; The coverage area calculation module is used to calculate the coverage area of the protection force covering the protection area based on the radius and center of all protection circles, the total set of polar angles corresponding to all protection circles, and the start and end coordinates of all line segments on the boundary of the protection area. Specifically, it calculates the coverage area of the protection force covering the protection area according to the formula. The calculation is performed, where A represents the coverage area of the support force. Indicates the total number of line segments. This represents the coordinates of the endpoint of the j-th line segment. This represents the coordinates of the starting point of the j-th line segment. Indicates the total number of support personnel. This represents the radius of the protection circle corresponding to the i-th protection force. This represents the number of start-end pairs in the total set of polar angles corresponding to the i-th support force's support circle. This represents the k-th starting polar angle in the total set of polar angles corresponding to the polar circle of the i-th support force. This represents the k-th endpoint polar angle in the total set of polar angles corresponding to the polar circles of the i-th support force. This represents the coordinates of the center of the protection circle corresponding to the i-th protection force.
7. An electronic device, characterized in that, include: A memory and a processor, the memory being used to store a computer program, the processor running the computer program to cause the electronic device to perform the method for calculating the coverage area according to any one of claims 1 to 5.
8. A computer-readable storage medium, characterized in that, It stores a computer program that, when executed by a processor, implements the method for calculating the coverage area as described in any one of claims 1 to 5.