CN116032823A - Method and system for rapidly calculating scattering communication link based on digital elevation map - Google Patents
Method and system for rapidly calculating scattering communication link based on digital elevation map Download PDFInfo
- Publication number
- CN116032823A CN116032823A CN202310001061.9A CN202310001061A CN116032823A CN 116032823 A CN116032823 A CN 116032823A CN 202310001061 A CN202310001061 A CN 202310001061A CN 116032823 A CN116032823 A CN 116032823A
- Authority
- CN
- China
- Prior art keywords
- point
- determining
- slope
- communication link
- straight line
- 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.)
- Granted
Links
- 238000004891 communication Methods 0.000 title claims abstract description 65
- 238000000034 method Methods 0.000 title claims description 34
- 238000004364 calculation method Methods 0.000 claims abstract description 81
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 22
- 230000005540 biological transmission Effects 0.000 claims abstract description 12
- 238000004590 computer program Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 4
- 108010034145 Helminth Proteins Proteins 0.000 claims 1
- 244000000013 helminth Species 0.000 claims 1
- 230000006870 function Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 3
- 238000006467 substitution reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000012850 discrimination method Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention relates to a rapid calculation method and a rapid calculation system for a scattering communication link based on a digital elevation map, which belong to the field of scattering communication, and comprise the following steps: acquiring a digital elevation grid map and longitude and latitude coordinates of a receiving end and a transmitting end; determining integer coordinates and a linear equation corresponding to the receiving and transmitting ends on the grid map according to the grid map and the longitude and latitude coordinates; extracting a profile coordinate point by adopting a Blisenham algorithm according to the integer coordinate and the linear equation; determining the position and the elevation value of the highest obstacle relative to the receiving and transmitting ends through ratio calculation according to the profile coordinate points; determining a scattering angle according to the highest obstacle position and the elevation value; and determining scattering communication transmission loss according to the scattering angle, and completing scattering communication link calculation. The invention can improve the calculation efficiency and provides a good foundation for simultaneously carrying out large-scale scattered communication link calculation.
Description
Technical Field
The invention relates to the field of scattered communication, in particular to a method and a system for rapidly calculating a scattered communication link based on a digital elevation map.
Background
In the current scatter communication link calculation method, in the step of extracting profile elevation data, a straight line equation where two points are located is generally calculated, grid points passing through the straight line equation are searched to form a profile data set, and the profile data set comprises a large number of floating point number multiplication and division calculation. When performing link calculation, a large number of link calculation cannot be performed at the same time, and certain limitation exists. In the process of calculating the scattering angle after extracting the profile data, most of the current methods are to search the height of the highest obstacle by calculating the angle, so that a large amount of trigonometric function calculation is needed, and a lot of calculation resources are consumed.
Disclosure of Invention
The invention aims to provide a rapid calculation method of a scattering communication link based on a digital elevation map, which is characterized in that a Blisenham algorithm is introduced into the calculation of the scattering communication link to extract profile data, floating point multiplication and division calculation is replaced by integer addition and subtraction calculation, and a ratio substitution angle value method is used in the process of calculating a scattering angle, so that trigonometric function calculation is reduced, the whole algorithm optimizes the calculation thought, improves the calculation efficiency, and provides a good foundation for simultaneously carrying out large-scale scattering communication link calculation.
In order to achieve the above object, the present invention provides the following solutions:
in a first aspect, the present invention provides a method for rapidly calculating a scatter communication link based on a digital elevation map, comprising:
acquiring a digital elevation grid map and longitude and latitude coordinates of a receiving end and a transmitting end;
determining integer coordinates and a linear equation corresponding to the receiving and transmitting ends on the grid map according to the grid map and the longitude and latitude coordinates;
extracting a profile coordinate point by adopting a Blisenham algorithm according to the integer coordinate and the linear equation;
determining the position and the elevation value of the highest obstacle relative to the receiving and transmitting ends through ratio calculation according to the profile coordinate points;
determining a scattering angle according to the highest obstacle position and the elevation value;
and determining scattering communication transmission loss according to the scattering angle, and completing scattering communication link calculation.
Optionally, after the step of determining a scattering angle according to the highest obstacle position and the elevation value, before the step of determining a scattering communication transmission loss according to the scattering angle and completing the step of calculating a scattering communication link, the method further includes:
the distance and the communication frequency between the transmitting and receiving ends are determined.
Optionally, after the step of determining the integer coordinate and the linear equation corresponding to the transceiver end on the grid map according to the grid map and the longitude and latitude coordinate, before the step of extracting the profile coordinate point by using the brisen ham algorithm according to the integer coordinate and the linear equation, the method further includes:
determining the relative position relationship between the slope of the linear equation and the integer coordinates of the two receiving and transmitting ends, wherein the integer coordinates of the two receiving and transmitting ends are respectively marked as (x) 0 ,y 0 ) And (x) k ,y k )。
Optionally, the relative positional relationship of the integer coordinates of the two transceiver ends includes:
the slope of the straight line equation y=mx+b is 0.ltoreq.m < 1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Right side;
the slope of the straight line equation y=mx+b is m.gtoreq.1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Right side;
the slope m < -1 of the straight line equation y=mx+b, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Left side;
the slope of the straight line equation y=mx+b is-1.ltoreq.m < 0, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Left side;
the slope of the straight line equation y=mx+b is 0.ltoreq.m < 1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Left side;
the slope of the straight line equation y=mx+b is m.gtoreq.1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Left side;
the slope of the straight line equation y=mx+b is m < -1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Right side;
the slope of the straight line equation y=mx+b is-1.ltoreq.m < 0, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Right side.
In a second aspect, the present invention provides a digital elevation map based scatter communication link fast computing system comprising:
the data acquisition module is used for acquiring the digital elevation grid map and longitude and latitude coordinates of the receiving and transmitting ends;
the integer coordinate and linear equation determining module is used for determining integer coordinates and linear equations corresponding to the receiving and transmitting ends on the grid map according to the grid map and the longitude and latitude coordinates;
the profile coordinate point extraction module is used for extracting profile coordinate points by adopting a Blisenham algorithm according to the integer coordinates and the linear equation;
the highest obstacle position and elevation value determining module is used for determining the highest obstacle position and elevation value relative to the receiving and transmitting ends through ratio calculation according to the profile coordinate points;
the scattering angle determining module is used for determining a scattering angle according to the highest obstacle position and the elevation value;
and the scattering communication link calculation module is used for determining scattering communication transmission loss according to the scattering angle and finishing the scattering communication link calculation.
Optionally, between the scattering angle determining module and the scattering communication link calculating module, the method further includes:
and the distance and communication frequency determining module is used for determining the distance and the communication frequency between the two receiving and transmitting ends.
Optionally, between the integer coordinate and linear equation determining module and the section coordinate point extracting module, the method further includes:
a slope and relative position relation determining module for determining the relative position relation between the slope of the linear equation and the integer coordinates of the two receiving ends, wherein the integer coordinates of the two receiving ends are respectively marked as (x) 0 ,y 0 ) And (x) k ,y k )。
Optionally, the relative positional relationship of the integer coordinates of the two transceiver ends includes:
the slope of the straight line equation y=mx+b is 0.ltoreq.m < 1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Right side;
the slope of the straight line equation y=mx+b is m.gtoreq.1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Right side;
the slope m < -1 of the straight line equation y=mx+b, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Left side;
the slope of the straight line equation y=mx+b is-1.ltoreq.m < 0, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Left side;
the slope of the straight line equation y=mx+b is 0.ltoreq.m < 1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Left side;
the slope of the straight line equation y=mx+b is m.gtoreq.1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Left side;
the slope of the straight line equation y=mx+b is m < -1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Right side;
the slope of the straight line equation y=mx+b is-1.ltoreq.m < 0, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Right side.
In a third aspect, the present invention provides an electronic device comprising a memory for storing a computer program and a processor that runs the computer program to cause the electronic device to perform the digital elevation map based scatter communication link fast calculation method.
In a fourth aspect, the present invention provides a computer readable storage medium storing a computer program which when executed by a processor implements the digital elevation map based scatter communication link fast calculation method.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
according to the invention, the Blaston Ham algorithm is introduced into the scattering communication link calculation for extracting the profile data, floating point number multiplication and division calculation is replaced by integer addition and subtraction calculation, and a ratio substitution angle value method is used in the process of calculating the scattering angle, so that trigonometric function calculation is reduced, the whole algorithm optimizes the calculation thought, improves the calculation efficiency, and provides a good foundation for simultaneously carrying out large-scale scattering communication link calculation.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions of the prior art, the drawings that are needed in the embodiments will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method for rapidly calculating a scatter communication link based on a digital elevation map according to the present invention;
FIG. 2 is an extracted cross-sectional view of the Blaston Ham algorithm;
FIG. 3 is a diagram of a ratio search algorithm;
fig. 4 is a view showing a positional relationship between the transmitting and receiving ends.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The invention aims to provide a rapid calculation method of a scattering communication link based on a digital elevation map, which is characterized in that a Blisenham algorithm is introduced into the calculation of the scattering communication link to extract profile data, floating point multiplication and division calculation is replaced by integer addition and subtraction calculation, and a ratio substitution angle value method is used in the process of calculating a scattering angle, so that trigonometric function calculation is reduced, the whole algorithm optimizes the calculation thought, improves the calculation efficiency, and provides a good foundation for simultaneously carrying out large-scale scattering communication link calculation.
In order to achieve the aim of the optimization calculation, the invention firstly determines the coordinates (x 0 ,y 0 ),(x k ,y k ) The integer coordinate data set L= [ (x) between the receiving and transmitting ends is rapidly extracted by utilizing Blisenham algorithm 0 ,y 0 ),...,(x i ,y i ),...,(x k ,y k )](0 < i < k), extracting an elevation dataset H= [ H ] located on the coordinate dataset 0 ,…,h k ]Then the highest obstacle position is found out rapidly by using the ratio to replace the trigonometric function calculation method, and the horizontal angle alpha of the emission point is obtained t The same process obtains the horizontal angle alpha of the receiving point r The scattering angle θ can be found and the scattering link transmission loss value can be found using ITU R-p.617-5 recommendations.
In order that the above-recited objects, features and advantages of the present invention will become more readily apparent, a more particular description of the invention will be rendered by reference to the appended drawings and appended detailed description.
FIG. 1 is a flow chart of a method for rapidly calculating a scatter communication link based on a digital elevation map according to the present invention; as shown in fig. 1, a method for rapidly calculating a scattering communication link based on a digital elevation map includes:
step 1: and acquiring a digital elevation grid map and longitude and latitude coordinates of the receiving and transmitting ends.
Step 2: and determining integer coordinates and a linear equation corresponding to the receiving and transmitting ends on the grid map according to the grid map and the longitude and latitude coordinates.
The calculation process of the steps 1 to 2 is as follows:
according to longitude and latitude (x) t ,y t )，(x r ,y r ) To determine the coordinates (x 0 ,y 0 ),(x k ,y k ). The number of map grids is M, representing a degree of each grid of 1/M. The integer coordinate point (x) of the receiving and transmitting ends in the grid elevation map can be calculated by simple numerical calculation 0 ,y 0 ),(x k ,y k )。
Then two points (x 0 ,y 0 ),(x k ,y k ) The relative positional relationship is divided into 8 cases: (1) when the slope of the straight line y=mx+b is 0.ltoreq.m < 1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) On the right side, (2) when the slope of the straight line is m.gtoreq.1, and point (x) k ,y k ) At the point (x) 0 ,y 0 ) In the case of the right side, (3) when the slope of the straight line is m < -1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) In the left case, (4) when the slope of the straight line is-1.ltoreq.m < 0, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) In the left case, (5) when the slope of the straight line is 0.ltoreq.m < 1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) In the left case, (6) when the slope of the straight line is m.gtoreq.1, and point (x) k ,y k ) At the point (x) 0 ,y 0 ) In the left case, (7) when the slope of the straight line is m < -1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) In the case of the right side, (8) when the slope of the straight line is-1.ltoreq.m < 0, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Right hand case.
Step 3: and extracting a profile coordinate point by adopting a Blisenham algorithm according to the integer coordinate and the linear equation.
Specifically, the Blisenham algorithm is adopted to extract the profile elevation data, so that a great amount of floating point number multiplication and division calculation in the conventional algorithm for calculating the intersection point of the straight line and the grid can be converted into addition and subtraction calculation of the integral point coordinates, and meanwhile, fitting errors are not introduced, and the calculation efficiency is improved on the premise of ensuring the calculation accuracy.
Step 4: and determining the highest obstacle position and elevation value relative to the receiving and transmitting ends through ratio calculation according to the profile coordinate points.
Specifically, by utilizing the characteristic that the distance between coordinate points extracted by Blisenham is a constant and the characteristic that the tangent value corresponds to the angle, the highest obstacle position information and elevation value relative to a receiving and transmitting site are quickly searched by adopting a ratio calculation to replace a trigonometric function calculation method. The process of calculating the geometric parameters of the scattering link is simplified, and the calculation efficiency is improved.
The calculation process of the steps 3 to 4 is as follows:
first, consider (1) that, in order to obtain the coordinate data set between two points of the transceiver station, consider a certain point (x i ,y i ) From the properties of the grid elevation digital map, it is known that the next profile data point at the transmitting station is at (x i +1,y i )，(x i +1,y i +1) these two points. At x=x i At +1, the point on the straight line is y=m (x i +1) +b, from (x) i +1,y i )，(x i +1,y i +1) distance between two points d 1 ＝m(x i +1)+b-y i And d 2 ＝(y i +1)-m(x i +1) -b, and Δd=d 1 -d 2 ＝2m(x i +1)+2b-(2y i +1), when Δd > 0, is closer to point (x i +1,y i ) Therefore, the next point (x) i +1,y i ) When Δd < 0, the point is closer (x i +1,y i +1), so take the next point (x i +1,y i +1). Therefore, only the sign of Δd needs to be determined to determine the value of the next point.
To facilitate determination of the Δd symbol, an expression p is constructed i ＝Δx·Δd＝2Δy·x i -2Δx·y i +c, where Δy=y k -y 0 M=Δy/Δx, c=2Δy—Δx (2 b-1). Setting the starting point to the left and the receiving point to the right, (x) k -x 0 ) > 0, so p i With the same sign as Δd, further process p i+1 ＝2Δy·x i+1 -2Δx·y i+1 +c by subtracting two formulae while using x i+1 ＝x i +1 gives p i+1 ＝p i +2Δy-2Δx(y i+1 -y i ) Let the starting point of the straight line be (x 0 ,y 0 ) I.e. satisfy y 0 ＝mx 0 +b, thus obtaining p 0 By using the error discrimination method, the initial value p is obtained 0 =2Δy- Δx, when p i Y is equal to or greater than 0 i+1 ＝y i +1,x i+1 ＝x i +1,p i+1 ＝p i +2 (Δy- Δx), when p i When < 0, y i+1 ＝y i ,x i+1 ＝x i +1,p i+1 ＝p i +2Δy. It can be seen that step i+1 is only compared with the discrimination variable p of step i i The coordinate component differences delta y and delta x of the two endpoints of the correlation and receiving station are related, the operation only comprises the operations of integer addition and multiplication 2, and the multiplication 2 can be completed by shifting one bit left by the computer number, so that the calculation efficiency is greatly improved, and a coordinate point data set L= [ (x) between the two points of the receiving and receiving station is obtained 0 ,y 0 ),...,(x i ,y i ),...,(x k ,y k )](0 < i < k) and the elevation data set h= [ H ] at the coordinate point 1 ,...,h k ]And a distance set j= [ J ] from the emission point 1 ,...,j k ]。
When two points (x 0 ,y 0 ),(x k ,y k ) When the relative position relation is in other areas, the distance set and the elevation value set can be acquired by parameter symbol adjustment and adaptation to different conditions.
When acquiring the section data Gao Chengji H of the connecting line between the two points and the distance set J from the transmitting point, calculating the highest obstacle height H t . Known distance set j= [ J ] 1 ,...,j k ]Is an arithmetic series, and thus has a tangent value of
Can be replaced by the ratio +.>
Without affecting the order of size. The angular magnitude relationship also corresponds one-to-one with the tangential magnitude relationship. Finding the maximum value in the ratio set R to obtain the highest obstacle height h t At this time, the distance from the emission point is j t The integer point coordinates are (x t ,y t ). The highest obstacle height of the receiving station which can be obtained by the same method is h r Distance from receiving station is j r The integer point coordinates are (x r ,y r )。
Step 5: and determining a scattering angle according to the highest obstacle position and the elevation value.
Step 6: and determining scattering communication transmission loss according to the scattering angle, and completing scattering communication link calculation.
The specific calculation process of the step 5 to the step 6 is as follows:
when the highest obstacle information between the two receiving and transmitting ends is obtained, the scattering angle theta can be calculated according to the trigonometric function, and the basic transmission loss is proposed to be L=F+22logf+35 logtheta+17logd+L according to ITU_R-P.617-5 c -Y p Where F is frequency, θ is scattering angle (radian), d is distance between transceiver stations, and y=f+l is set c -Y p . At this time, the formula may be simplified to l=22 logf+35log θ+17 log+y, where y is related to factors such as communication region climate, antenna gain, and the like.
In addition, fig. 2 is a schematic diagram of the british algorithm for extracting a scattering communication link profile in the present invention, specifically including the following steps:
(1) According to the longitude and latitude coordinates (x) of two points of the transceiver station t ,y t )，(x r ,y r ) Calculating the integer coordinate point (x) 0 ,y 0 )，(x k ,y k ). The elevation data of the grid map is read as D (i, j), wherein i, j are integers. (x) 0 ,y 0 )，(x k ,y k ) The line between two points is y=mx+b, and the elevation value of the grid map in one grid is the same value, so that the line between the points (x i ,mx i The elevation value at +b) is point (x i ,y i ) Elevation value h of (2) i ＝D(x i ,y i )。
(2) Setting (x) 0 ,y 0 )，(x k ,y k ) The two-point position relationship is shown in the region (1) in fig. 4, at the point (x i ,y i ) The next point takes a value of (x i +1,y i )，(x i +1,y i +1) one of the two points. Point (x) i +1,m(x i +1) +b) and (x) i +1,y i ) The distance of the points is d 1 And (x) i +1,y i +1) Point distance d 2 . When d 1 ＞d 2 In this case, the intersection distance (x i +1,y i +1) point is near, so the next point takes on the value h i+1 ＝D(x i +1,y i +1), when d 1 ＜d 2 In this case, the intersection distance (x i +1,y i ) The next point is therefore given a value of h i+1 ＝D(x i +1,y i )。
(3) Judgment d 1 -d 2 The sign of (a) may construct an expression p i ＝Δx·Δd＝2Δy·x i -2Δx·y i +c, where Δy=y k -y 0 M=Δy/Δx, c=2Δy—Δx (2 b-1), assuming that the value advance direction is the positive x-axis direction, Δx is positive, indicating p i The same symbols as Δd.
(4) Calculation of p i+1 ＝2Δy·x i+1 -2Δx·y i+1 +c and p i ＝2Δy·x i -2Δx·y i +c, subtracting x i+1 ＝x i +1 is taken into the formula to give p i+1 ＝p i +2Δy-2Δx(y i+1 -y i ) The straight line passes through an integer coordinate point (x 0 ,y 0 ) I.e. y 0 ＝mx 0 +b, the initial value p can be obtained 0 ＝2Δy-Δx。
(5) When p is i When 0 is not less than 0, Δd is not less than 0, and represents a straight line on-line x=x i At +1, a point closer to (x i +1,y i +1), thus the next point takes on the value h i+1 ＝D(x i +1,y i +1), then p is calculated i+1 ＝p i +2Δy; when p is i When < 0, Δd < 0 is represented, and the straight line is represented as on-line x=x i At +1, a point closer to (x i +1,y i ) Therefore, the next point takes the value h i+1 ＝D(x i +1,y i ) At this time, p is then calculated i+1 ＝p i +2Δy。
(6) Continuing to judge that the straight line is x=x i Δd symbol at +2, selecting elevation value point, calculating p i+2 . And the like, until the integer coordinates of the receiving points are valued, obtaining a coordinate point data set L= [ (x) 0 ,y 0 ),...,(x i ,y i ),…,(x k ,y k )]. At the same time, an elevation data set H= [ H ] 1 ,...,h k ]And a distance set j= [ J ] from the transmitting station 1 ,...,j k ]。
FIG. 3 is a schematic diagram of a method for quickly searching the highest obstacle of a transmitting station by adopting a comparison ratio in the patent of the invention, which specifically comprises the following steps:
(1) The highest obstacle position required to block the transmitting station is the obstacle with the highest elevation angle, i.e. the maximum value of alpha. It is understood that the larger the angle, the larger the tangent value is, within the [0,90] degree interval.
(2) Tangent data set is
From the characteristics of the grid map, the sampling interval is a fixed value representingDistance set between two points j= [ J ] 1 ,...,j k ]Is an arithmetic series. Thus R is w Size order of (2)
In agreement, the comparison gives a maximum value of +.>
Tangent value here->
I.e. the maximum tangent, i.e. representing the angle here +>
Is the largest angle.
(3) Obtaining a point (x) t ,y t ) A coordinate point h which is the highest obstacle of the emission points t Is the elevation value of the highest obstacle, alpha t Is the horizontal angle of the emission point. The highest obstacle coordinate point (x r ,y r )，h r Is the elevation value of the highest obstacle, alpha r Is the horizontal angle of the receiving point.
(4) H is known to be t ，α t ，h r ，α r And the distance j between the two transmitting and receiving points k The transmission loss L is rapidly calculated by setting the scattering communication frequency f according to the ITU R-p.617-5 recommendation l=22 logf+35logθ+17 logd+y formula.
Based on the method, the invention also provides a rapid computing system of the scattering communication link based on the digital elevation map, which comprises the following steps:
and the data acquisition module is used for acquiring the digital elevation grid map and longitude and latitude coordinates of the receiving and transmitting ends.
And the integer coordinate and linear equation determining module is used for determining integer coordinates and linear equations corresponding to the receiving and transmitting ends on the grid map according to the grid map and the longitude and latitude coordinates.
And the profile coordinate point extraction module is used for extracting profile coordinate points by adopting a Blisenham algorithm according to the integer coordinates and the linear equation.
And the highest obstacle position and elevation value determining module is used for determining the highest obstacle position and elevation value relative to the receiving and transmitting ends through ratio calculation according to the profile coordinate points.
And the scattering angle determining module is used for determining the scattering angle according to the highest obstacle position and the elevation value.
And the scattering communication link calculation module is used for determining scattering communication transmission loss according to the scattering angle and finishing the scattering communication link calculation.
The invention also discloses the following technical effects:
the invention adopts the Blisenham algorithm to extract the profile elevation data, can convert a great amount of floating point number multiplication and division calculation in the common algorithm for calculating the intersection point of the straight line and the grid into addition and subtraction calculation of the integral point coordinates, does not introduce fitting errors, and improves the calculation efficiency on the premise of ensuring the calculation accuracy.
The invention uses the characteristic that the distance between coordinate points extracted by Blisen Hanm is a constant and the characteristic that the tangent value corresponds to the angle, adopts the ratio calculation to replace the trigonometric function calculation method, and rapidly searches the highest obstacle position information and elevation value relative to the receiving and transmitting ends. The process of calculating the geometric parameters of the scattering link is simplified, and the calculation efficiency is improved.
In the present specification, each embodiment is described in a progressive manner, and each embodiment is mainly described in a different point from other embodiments, and identical and similar parts between the embodiments are all enough to refer to each other. For the system disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant points refer to the description of the method section.
The principles and embodiments of the present invention have been described herein with reference to specific examples, the description of which is intended only to assist in understanding the methods of the present invention and the core ideas thereof; also, it is within the scope of the present invention to be modified by those of ordinary skill in the art in light of the present teachings. In view of the foregoing, this description should not be construed as limiting the invention.
Claims (10)
1. A method for rapidly calculating a scatter communication link based on a digital elevation map, comprising:
acquiring a digital elevation grid map and longitude and latitude coordinates of a receiving end and a transmitting end;
determining integer coordinates and a linear equation corresponding to the receiving and transmitting ends on the grid map according to the grid map and the longitude and latitude coordinates;
extracting a profile coordinate point by adopting a Blisenham algorithm according to the integer coordinate and the linear equation;
determining the position and the elevation value of the highest obstacle relative to the receiving and transmitting ends through ratio calculation according to the profile coordinate points;
determining a scattering angle according to the highest obstacle position and the elevation value;
and determining scattering communication transmission loss according to the scattering angle, and completing scattering communication link calculation.
2. The rapid calculation method of a scattered communication link based on a digital elevation map according to claim 1, wherein after said step of determining a scattering angle based on said highest obstacle position and elevation value, determining a scattered communication transmission loss based on said scattering angle, and before said step of completing the calculation of the scattered communication link, further comprises:
the distance and the communication frequency between the transmitting and receiving ends are determined.
3. The method according to claim 1, wherein after the step of determining the integer coordinates and the linear equation corresponding to the transceiver end on the grid map according to the grid map and the longitude and latitude coordinates, the step of extracting the profile coordinate point according to the integer coordinates and the linear equation by using the british helminth algorithm, further comprises:
determining the slope sum of the linear equationThe relative position relationship of the integer coordinates of the two ends is marked as (x) 0 ,y 0 ) And (x) k ,y k )。
4. The rapid computing method of scattering communication link based on digital elevation map according to claim 3, wherein the relative positional relationship of the integer coordinates of the transmitting and receiving ends comprises:
the slope of the straight line equation y=mx+b is 0.ltoreq.m < 1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Right side;
the slope of the straight line equation y=mx+b is m.gtoreq.1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Right side;
the slope m < -1 of the straight line equation y=mx+b, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Left side;
the slope of the straight line equation y=mx+b is-1.ltoreq.m < 0, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Left side;
the slope of the straight line equation y=mx+b is 0.ltoreq.m < 1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Left side;
the slope of the straight line equation y=mx+b is m.gtoreq.1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Left side;
the slope of the straight line equation y=mx+b is m < -1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Right side;
the slope of the straight line equation y=mx+b is-1.ltoreq.m < 0, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Right side.
5. A digital elevation map based scatter communication link fast computing system, comprising:
the data acquisition module is used for acquiring the digital elevation grid map and longitude and latitude coordinates of the receiving and transmitting ends;
the integer coordinate and linear equation determining module is used for determining integer coordinates and linear equations corresponding to the receiving and transmitting ends on the grid map according to the grid map and the longitude and latitude coordinates;
the profile coordinate point extraction module is used for extracting profile coordinate points by adopting a Blisenham algorithm according to the integer coordinates and the linear equation;
the highest obstacle position and elevation value determining module is used for determining the highest obstacle position and elevation value relative to the receiving and transmitting ends through ratio calculation according to the profile coordinate points;
the scattering angle determining module is used for determining a scattering angle according to the highest obstacle position and the elevation value;
and the scattering communication link calculation module is used for determining scattering communication transmission loss according to the scattering angle and finishing the scattering communication link calculation.
6. The digital elevation map based scatter communication link fast calculation system of claim 5, further comprising, between said scatter angle determination module and scatter communication link calculation module:
and the distance and communication frequency determining module is used for determining the distance and the communication frequency between the two receiving and transmitting ends.
7. The rapid computing system of a digital elevation map based scatter communication link of claim 5, wherein between the integer coordinate and linear equation determination module and the profile coordinate point extraction module, further comprising:
a slope and relative position relation determining module for determining the relative position relation between the slope of the linear equation and the integer coordinates of the two receiving ends, wherein the integer coordinates of the two receiving ends are respectively marked as (x) 0 ,y 0 ) And (x) k ,y k )。
8. The rapid computing system of a digital elevation map based scatter communication link of claim 7, wherein the relative positional relationship of the integer coordinates of the transceiver ends comprises:
the slope of the straight line equation y=mx+b is 0.ltoreq.m < 1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Right side;
the slope of the straight line equation y=mx+b is m.gtoreq.1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Right side;
the slope m < -1 of the straight line equation y=mx+b, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Left side;
the slope of the straight line equation y=mx+b is-1.ltoreq.m < 0, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Left side;
the slope of the straight line equation y=mx+b is 0.ltoreq.m < 1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Left side;
the slope of the straight line equation y=mx+b is m.gtoreq.1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Left side;
the slope of the straight line equation y=mx+b is m < -1, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Right side;
the slope of the straight line equation y=mx+b is-1.ltoreq.m < 0, and the point (x k ,y k ) At the point (x) 0 ,y 0 ) Right side.
9. An electronic device comprising a memory for storing a computer program and a processor that runs the computer program to cause the electronic device to perform the digital elevation map based scatter communication link fast calculation method according to any one of claims 1-4.
10. A computer-readable storage medium, characterized in that it stores a computer program which, when executed by a processor, implements the digital elevation map based scatter communication link fast calculation method according to any one of claims 1-4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310001061.9A CN116032823B (en) | 2023-01-03 | 2023-01-03 | Method and system for rapidly calculating scattering communication link based on digital elevation map |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310001061.9A CN116032823B (en) | 2023-01-03 | 2023-01-03 | Method and system for rapidly calculating scattering communication link based on digital elevation map |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116032823A true CN116032823A (en) | 2023-04-28 |
| CN116032823B CN116032823B (en) | 2023-11-28 |
Family
ID=86070225
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310001061.9A Active CN116032823B (en) | 2023-01-03 | 2023-01-03 | Method and system for rapidly calculating scattering communication link based on digital elevation map |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116032823B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116634474A (en) * | 2023-07-05 | 2023-08-22 | 中国人民解放军61905部队 | Method, system, equipment and medium for predicting communication coverage of large-area cluster |
| CN117135639A (en) * | 2023-08-09 | 2023-11-28 | 中国人民解放军61905部队 | Wireless communication coverage prediction method, system and equipment based on neural network |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103559356A (en) * | 2013-11-11 | 2014-02-05 | 中国电子科技集团公司第五十四研究所 | GIS-based (Geographic Information System-based) scatter link planning method |
| US20170062900A1 (en) * | 2015-08-26 | 2017-03-02 | Telecommunication Systems, Inc. | Troposcatter antenna pointing |
| CN108173591A (en) * | 2017-12-29 | 2018-06-15 | 中国人民解放军空军工程大学 | A kind of cognition tropospheric scatter communication system and method, wireless communication system |
| CN114499652A (en) * | 2022-02-28 | 2022-05-13 | 中国电子科技集团公司第五十四研究所 | Power control method and control terminal of beyond-line-of-sight wireless communication system |
| WO2022233169A1 (en) * | 2021-05-06 | 2022-11-10 | 山东省科学院海洋仪器仪表研究所 | Method for calculating forward acoustic scattering doppler frequency shift of underwater moving target under irregular trajectory |
-
2023
- 2023-01-03 CN CN202310001061.9A patent/CN116032823B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103559356A (en) * | 2013-11-11 | 2014-02-05 | 中国电子科技集团公司第五十四研究所 | GIS-based (Geographic Information System-based) scatter link planning method |
| US20170062900A1 (en) * | 2015-08-26 | 2017-03-02 | Telecommunication Systems, Inc. | Troposcatter antenna pointing |
| CN108173591A (en) * | 2017-12-29 | 2018-06-15 | 中国人民解放军空军工程大学 | A kind of cognition tropospheric scatter communication system and method, wireless communication system |
| WO2022233169A1 (en) * | 2021-05-06 | 2022-11-10 | 山东省科学院海洋仪器仪表研究所 | Method for calculating forward acoustic scattering doppler frequency shift of underwater moving target under irregular trajectory |
| CN114499652A (en) * | 2022-02-28 | 2022-05-13 | 中国电子科技集团公司第五十四研究所 | Power control method and control terminal of beyond-line-of-sight wireless communication system |
Non-Patent Citations (2)
| Title |
|---|
| 李红雨 等: "基于MapX 的散射通信链路分析系统的开发", 无线电通信技术, vol. 37, no. 1, pages 53 - 55 * |
| 李红雨 等: "基于MapX的散射通信链路分析系统的开发", 无线电通信技术, vol. 37, no. 1, pages 53 - 55 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116634474A (en) * | 2023-07-05 | 2023-08-22 | 中国人民解放军61905部队 | Method, system, equipment and medium for predicting communication coverage of large-area cluster |
| CN116634474B (en) * | 2023-07-05 | 2024-06-25 | 中国人民解放军61905部队 | Method, system, equipment and medium for predicting communication coverage of large-area cluster |
| CN117135639A (en) * | 2023-08-09 | 2023-11-28 | 中国人民解放军61905部队 | Wireless communication coverage prediction method, system and equipment based on neural network |
| CN117135639B (en) * | 2023-08-09 | 2024-03-08 | 中国人民解放军61905部队 | Wireless communication coverage prediction method, system and equipment based on neural network |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116032823B (en) | 2023-11-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN116032823B (en) | Method and system for rapidly calculating scattering communication link based on digital elevation map | |
| CN112652036B (en) | Road data processing method, device, equipment and storage medium | |
| CN104793224A (en) | GPS positioning rectification method and device | |
| CN102799763B (en) | A kind of based on a cloud attitude standardized some cloud line feature extraction method | |
| CN116559928B (en) | Pose information determining method, device and equipment of laser radar and storage medium | |
| CN107341284B (en) | Bidirectional parabolic equation method for predicting low-frequency radio wave propagation characteristics with high precision | |
| CN104869641A (en) | Wi-Fi indoor positioning method based on AP optimization | |
| CN109522385B (en) | Method for judging M-N matching mode of multi-scale road network | |
| CN112526470A (en) | Method and device for calibrating radar parameters, electronic equipment and storage medium | |
| Huang et al. | MAPS: Indoor localization algorithm based on multiple AP selection | |
| CN110096564B (en) | Route point positioning method, device and system based on BIM + GIS | |
| CN114283343B (en) | Map updating method, training method and device based on remote sensing satellite image | |
| CN116309026A (en) | Point cloud registration method and system based on statistical local feature description and matching | |
| CN113239136B (en) | Data processing method, device, equipment and medium | |
| US20190347835A1 (en) | Maze solving method based on line-surface spatial relation | |
| CN109067678B (en) | WFRFT signal cascade modulation identification method based on high-order cumulant and wireless communication system | |
| CN104635206B (en) | A kind of method and device of wireless location | |
| CN115952248B (en) | Pose processing method, device, equipment, medium and product of terminal equipment | |
| CN112130166A (en) | AGV positioning method and device based on reflector network | |
| Wu et al. | Research on WiFi Location Fingerprint Positioning Algorithm Based on DPC‐FCM Clustering | |
| CN116704024A (en) | Vehicle-mounted laser point cloud pose graph optimization method and system combining multiclass constraints | |
| US20220414929A1 (en) | Portable device positioning data processing method and apparatus, device, and storage medium | |
| CN104199937A (en) | Multi-website POI position mapping method and device | |
| CN104156938A (en) | Image communication area description method, and application method thereof in image registration | |
| CN112614162A (en) | Indoor vision rapid matching positioning method and system based on space optimization strategy |
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 |