CN111475916A - DEM-based radar shielding angle calculation method and system - Google Patents
DEM-based radar shielding angle calculation method and system Download PDFInfo
- Publication number
- CN111475916A CN111475916A CN202010156881.1A CN202010156881A CN111475916A CN 111475916 A CN111475916 A CN 111475916A CN 202010156881 A CN202010156881 A CN 202010156881A CN 111475916 A CN111475916 A CN 111475916A
- Authority
- CN
- China
- Prior art keywords
- radar
- dem
- target
- distance
- azimuth
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F9/00—Arrangements for program control, e.g. control units
- G06F9/06—Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
- G06F9/46—Multiprogramming arrangements
- G06F9/50—Allocation of resources, e.g. of the central processing unit [CPU]
- G06F9/5005—Allocation of resources, e.g. of the central processing unit [CPU] to service a request
- G06F9/5027—Allocation of resources, e.g. of the central processing unit [CPU] to service a request the resource being a machine, e.g. CPUs, Servers, Terminals
-
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The embodiment of the application discloses a method and a system for calculating a radar shielding angle based on DEM, wherein the method comprises the following steps: preprocessing a target DEM file; calculating a longitude and latitude range of a target radar detection range according to the target radar deployment position and the detection distance, and acquiring a target DEM file covered by the longitude and latitude range; determining a limit distance according to the detection distance and a line-of-sight distance, wherein the line-of-sight distance is determined according to the height of a target radar and the highest point of the terrain; and sampling and taking points in the azimuth of 0-360 degrees at a set angle interval by taking the deployment position of the target radar as an origin, the limit distance as a radius and the true north as an azimuth, and calculating the radar shielding angle in each azimuth. The accuracy is guaranteed, meanwhile, the calculation amount is reduced, meanwhile, the GPU is used for performing parallel calculation, high-performance hardware resources are fully utilized, and calculation time is reduced.
Description
Technical Field
The embodiment of the application relates to the technical field of radars, in particular to a method and a system for calculating a radar shielding angle based on a DEM.
Background
The military radar is used for aiming at finding out an object of an oncoming enemy in the process of battle, and the target finding capability of the military radar is not only related to the performance of the radar and the characteristics of the oncoming object, but also closely related to the surrounding environment of a radar deployment position. The shielding angle is an included angle between a connecting line of the sight line and the top end of the barrier and a horizon line, and when the elevation angle of the radar antenna is smaller than the included angle, the radar cannot find a target due to the fact that the radar is shielded by the terrain. And a detection blind area formed by the landform blocking radar signal propagation can be obtained through the calculation of the shielding angle.
In a military simulation system, in order to reasonably deploy a radar to improve the operational capability of the radar, a Digital Elevation Model (DEM) is usually used for acquiring terrain data to calculate a radar shielding angle, so that simulation personnel can visually know the shielding condition of a position environment around a radar deployment area. The current shading angle calculation has problems: firstly, because DEM data volume is huge, there are the loading consuming time big, take up many scheduling problems of memory when using. Secondly, the calculation of the shielding angle needs a large amount of computation of terrain sampling points, and the existing method has large calculation time consumption.
The problems of large DEM loading time consumption, large memory occupation and large shading angle calculation time consumption are solved.
Disclosure of Invention
Therefore, the embodiment of the application provides a method and a system for calculating a radar shielding angle based on a DEM (digital elevation model), so as to solve the problems that the DEM in the prior art is large in loading time consumption, large in occupied memory and large in shielding angle calculation time consumption.
In order to achieve the above object, the embodiments of the present application provide the following technical solutions:
according to a first aspect of embodiments of the present application, there is provided a method for calculating a radar shielding angle based on DEM, the method including:
preprocessing a target DEM file, wherein the preprocessing comprises binary conversion of a file format and addition of a two-dimensional index identifier;
calculating a longitude and latitude range of a target radar detection range according to the target radar deployment position and the detection distance, and acquiring a target DEM file covered by the longitude and latitude range;
determining a limit distance according to the detection distance and a line-of-sight distance, wherein the line-of-sight distance is determined according to the height of a target radar and the highest point of the terrain;
sampling and taking points in the azimuth of 0-360 degrees at a set angle interval by taking the deployment position of the target radar as an origin, the limit distance as a radius and the included angle with the true north as an azimuth, and calculating the radar shielding angle in each azimuth;
the calculation of the radar shielding angles in all directions is divided into two layers by adopting a GPU (graphics processing unit) to be executed in parallel, the first layer is arranged among a plurality of DEM terrain blocks in parallel, and the second layer is arranged among a plurality of sampling points on each DEM terrain block in parallel.
Optionally, the calculating a radar shielding angle in each direction by taking the target radar deployment position as an origin, the limit distance as a radius, and an included angle with the true north as an azimuth, sampling and taking points in the directions of 0 degree to 360 degrees at a set angle interval, includes:
obtaining an azimuth line segment according to the origin, the radius and the azimuth angle;
sampling and taking points on the azimuth line segment at set distance intervals;
calculating an included angle between a connecting line of the surface position of each sampling point and the position of the target radar and the ground plane;
and determining the maximum value of the included angle of each sampling point as a radar shielding angle.
Optionally, the set distance interval of the sampling points is half of the accuracy of the topographic data, the set distance interval of the distant sampling points is twice of the accuracy of the topographic data, and the topographic data is a longitude and latitude range.
Optionally, the method further comprises: and sending the radar shielding angles in each direction from the GPU to the CPU, drawing the radar shielding angles in a radar chart mode, and releasing the memory occupied by the target DEM file.
Optionally, the information of the target DEM file includes a DEM file storage location, a start longitude, a start latitude, an end longitude, an end latitude, a highest point elevation, and a azimuth line segment sampling point.
According to a second aspect of embodiments of the present application, there is provided a DEM-based radar shading angle calculation system, the system including:
the DEM file management module is used for preprocessing a target DEM file, wherein the preprocessing comprises binary conversion of a file format and addition of a two-dimensional index identifier;
the topographic data determining module is used for calculating a longitude and latitude range of a target radar detection range according to the target radar deployment position and the detection distance, and acquiring a target DEM file covered by the longitude and latitude range;
the distance determining module is used for determining a limiting distance according to the detection distance and a line-of-sight distance, wherein the line-of-sight distance is determined according to the height of a target radar and the highest point of a terrain;
the radar shielding angle determining module is used for sampling and taking points in the directions of 0-360 degrees at set angle intervals by taking the deployment position of the target radar as an origin, the limit distance as a radius and a true north as an azimuth angle, and calculating the radar shielding angle in each direction; the calculation of the radar shielding angles in all directions is divided into two layers by adopting a GPU (graphics processing unit) to be executed in parallel, the first layer is arranged among a plurality of DEM terrain blocks in parallel, and the second layer is arranged among a plurality of sampling points on each DEM terrain block in parallel.
Optionally, the radar shielding angle determining module is further configured to: obtaining an azimuth line segment according to the origin, the radius and the azimuth angle; sampling and taking points on the azimuth line segment at set distance intervals; calculating an included angle between a connecting line of the surface position of each sampling point and the position of the target radar and the ground plane; and determining the maximum value of the included angle of each sampling point as a radar shielding angle.
Optionally, the set distance interval of the sampling points is half of the accuracy of the topographic data, the set distance interval of the distant sampling points is twice of the accuracy of the topographic data, and the topographic data is a longitude and latitude range.
Optionally, the system further comprises: and the sending module is used for sending the radar shielding angles on each direction from the GPU to the CPU, drawing the radar shielding angles in a radar map mode, and releasing the memory occupied by the target DEM file.
Optionally, the information of the target DEM file includes a DEM file storage location, a start longitude, a start latitude, an end longitude, an end latitude, a highest point elevation, and a azimuth line segment sampling point.
In summary, the embodiment of the application provides a method and a system for calculating a radar shielding angle based on a DEM, which perform preprocessing aiming at a target DEM file, wherein the preprocessing comprises binary conversion of a file format and addition of a two-dimensional index identifier; calculating a longitude and latitude range of a target radar detection range according to the target radar deployment position and the detection distance, and acquiring a target DEM file covered by the longitude and latitude range; determining a limit distance according to the detection distance and a line-of-sight distance, wherein the line-of-sight distance is determined according to the height of a target radar and the highest point of the terrain; sampling and taking points in the azimuth of 0-360 degrees at a set angle interval by taking the deployment position of the target radar as an origin, the limit distance as a radius and the included angle with the true north as an azimuth, and calculating the radar shielding angle in each azimuth; the calculation of the radar shielding angles in all directions is divided into two layers by adopting a GPU (graphics processing unit) to be executed in parallel, the first layer is arranged among a plurality of DEM terrain blocks in parallel, and the second layer is arranged among a plurality of sampling points on each DEM terrain block in parallel. The time consumption for reading the file is effectively reduced through DEM file conversion, the memory occupation is reduced through DEM file management, the calculation amount is reduced while the precision is ensured through the shielding angle calculation, meanwhile, the GPU is used for parallel calculation, the high-performance hardware resources are fully utilized, the time consumption for calculation is reduced, and the purpose of rapid calculation is achieved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below. It should be apparent that the drawings in the following description are merely exemplary, and that other embodiments can be derived from the drawings provided by those of ordinary skill in the art without inventive effort.
The structures, ratios, sizes, and the like shown in the present specification are only used for matching with the contents disclosed in the specification, so that those skilled in the art can understand and read the present invention, and do not limit the conditions for implementing the present invention, so that the present invention has no technical significance, and any structural modifications, changes in the ratio relationship, or adjustments of the sizes, without affecting the functions and purposes of the present invention, should still fall within the scope of the present invention.
Fig. 1 is a schematic flowchart of a method for calculating a radar shielding angle based on a DEM according to an embodiment of the present disclosure;
FIG. 2 is a schematic view of a shielding angle provided in an embodiment of the present application;
fig. 3 is a schematic diagram of a calculation result of a shading angle according to an embodiment of the present application;
fig. 4 is a block diagram of a system for calculating a radar shielding angle based on DEM according to an embodiment of the present disclosure.
Detailed Description
The present invention is described in terms of particular embodiments, other advantages and features of the invention will become apparent to those skilled in the art from the following disclosure, and it is to be understood that the described embodiments are merely exemplary of the invention and that it is not intended to limit the invention to the particular embodiments disclosed. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Fig. 1 shows a method for calculating a radar shielding angle based on a DEM according to an embodiment of the present application, which solves the problems of long DEM loading time, large memory occupation, and long shielding angle calculation time, and as shown in fig. 1, the method includes the following steps:
step 101: and preprocessing the target DEM file, wherein the preprocessing comprises binary conversion of a file format and addition of a two-dimensional index identifier.
Step 102: and calculating the longitude and latitude range of the detection range of the target radar according to the deployment position and the detection distance of the target radar, and acquiring a target DEM file covered by the longitude and latitude range.
Step 103: and determining a limit distance according to the detection distance and a line-of-sight distance, wherein the line-of-sight distance is determined according to the height of the target radar and the highest point of the terrain.
Step 104: sampling and taking points in the azimuth of 0-360 degrees at a set angle interval by taking the deployment position of the target radar as an origin, the limit distance as a radius and the included angle with the true north as an azimuth, and calculating the radar shielding angle in each azimuth; the calculation of the radar shielding angles in all directions is divided into two layers by adopting a GPU (graphics processing unit) to be executed in parallel, the first layer is arranged among a plurality of DEM terrain blocks in parallel, and the second layer is arranged among a plurality of sampling points on each DEM terrain block in parallel.
In a possible implementation, in step 101, the DEM file is subjected to a binary stream conversion process, and the binary stream file is directly read during calculation to reduce time consumption. Taking a DEM file with an asc format with a grid size of 6000 × 6000 as an example, defining each data in the DEM file as a double length write-in binary stream file, and sequentially writing the data as follows: row number (6000), column number (6000), starting longitude (70.0), starting latitude (40.0), grid cell size (0.000833), invalid value (-9999), 6000 x 6000 height values. For a DEM file with the same content, it takes 89s to load an asc file, and it takes 0.016s to load a binary stream file.
Each DEM file comprises elevation data of a specific area, and if all DEM files in the global range are loaded at the same time, the problem of overlarge memory cost is caused. The method sequentially traverses the topographic information of all DEM files during initialization, establishes a two-dimensional DEM file index according to the longitude and latitude sequence, then loads a plurality of related DEM files as required during calculation, and releases the files in time after use, thereby effectively reducing the memory occupation.
Each DEM file information includes DEM file storage location, start longitude, start latitude, end longitude, end latitude, highest point elevation, azimuth line segment sampling point, and the like. The size of the two-dimensional matrix is dynamically set according to the longitude and latitude span of the DEM file, the number of rows is 360/longitude span, and the number of columns is 180/latitude span.
Taking the example that each DEM file covers a longitude and latitude range area of 5 × 5, a two-dimensional matrix of 72 × 36 is defined to store global DEM information (longitude range 0-360, latitude range-90- +90), and the index of each DEM file in two-dimensional evidence is as follows: row number i ceil (start longitude/5); column number j equals ceil (starting latitude/5 + 18).
In a possible implementation manner, in step 102, for the shielding angle of each target radar, a longitude and latitude range of a detection range of the target radar is calculated according to the deployment position and the detection distance of the target radar, a covered DEM file of the longitude and latitude range is determined, and the covered DEM file is sent to a memory to calculate the shielding angle. Because the DEM data is usually large, the time for loading the data file is large. According to the method, each DEM file is converted into a binary stream file before calculation, and the binary stream file is directly read during calculation, so that time consumption can be effectively reduced. Defining each data in the DEM file as a double length and writing the double length into a binary stream file, wherein the written data sequentially comprises the following data: number of rows, number of columns, starting longitude, starting latitude, grid cell size, invalid value, number of rows x number of columns, number of height values.
And initializing, sequentially traversing all DEM file information, and establishing a two-dimensional DEM file index according to the topographic information, so that the related DEM file is loaded during calculation and released when the DEM file is used up. The DEM file information comprises a DEM file storage position, topographic information, a highest point elevation and a azimuth line section sampling point. The terrain information includes a start longitude, a start latitude, an end longitude, and an end latitude.
In a possible implementation manner, in step 102, for each radar for which the shielding angle needs to be calculated, a longitude and latitude range of a radar detection range is calculated according to a radar deployment position and a detection distance, then DEM file information covering the longitude and latitude range is determined, and the DEM file information is read into an internal memory for use in the calculation of the shielding angle.
In one possible implementation, in step 103, the line-of-sight distance is calculated from the radar height h1 and the highest point h2 of the terrain, taking into account the influence of the earth curvature and the influence of the atmosphere on the propagation of electromagnetic waves, and the magnitudes of both the detection distance and the line-of-sight distance are compared, taking the smaller data as the limit distance R. The line of sight distance
In a possible implementation, in step 104, as shown in fig. 2, the sampling points at 0 to 360 degrees of azimuth at set angle intervals with the target radar deployment position as an origin, the limiting distance as a radius, and the north as an azimuth, and calculating the radar shielding angle in each azimuth includes:
obtaining an azimuth line segment according to the origin, the radius and the azimuth angle; sampling and taking points on the azimuth line segment at set distance intervals; calculating an included angle between a connecting line of the surface position of each sampling point and the position of the target radar and the ground plane; and determining the maximum value of the included angle of each sampling point as a radar shielding angle.
In one possible embodiment, the set distance interval of the sampling points is half of the accuracy of the topographic data, the set distance interval of the distant sampling points is twice the accuracy of the topographic data, and the topographic data is the latitude and longitude range. The smaller the sampling distance interval setting, the higher the calculation accuracy, and therefore the larger the amount of calculation generated. Because the calculation accuracy of the shielding angle is related to the accuracy of the topographic data, the sampling distance interval delta d is set to be half of the accuracy of the topographic data, so that the full coverage of the topographic data points is ensured, and the calculated amount is controlled to be minimum. Meanwhile, the influence of terrain shielding is smaller when the distance is farther due to the curvature of the earth, and the distance sampling distance interval is set to be 2 delta d, so that the unnecessary calculation amount is reduced under the condition of ensuring the accuracy.
Taking DEM data with the precision of 90 meters as an example, the short-distance sampling interval is set to be 45 meters, and the long-distance sampling interval is set to be 90 meters. The terrain data uses geographic coordinates and the shading angle calculation uses a geocentric rectangular coordinate system.
The conversion from the geographic coordinate system to the geocentric coordinate system is X ═ N + H ═ cos (b) · cos (L), Y ═ N + H) · cos (b) · sin (L), and Z ═ N · (1-e _2_ C) + H) · sin (b).
The conversion relationship from the geocentric coordinate system to the geographic coordinate system is converted according to the following formulas (1), (2) and (3):
wherein, a is 6378137.0, b is 6356752.3142,
e2C＝0.00669437999013，
in one possible embodiment, the method further comprises: and sending the radar shielding angles in each direction from the GPU to the CPU, drawing the radar shielding angles in a radar chart mode as shown in figure 3, and releasing the memory occupied by the target DEM file. Due to the fact that the shading angles are large in calculation amount and independent from each other, the method further reduces calculation time by using GPU parallel calculation. The method is divided into two layers to be executed in parallel, a first layer of parallel operation is carried out among a plurality of DEM terrain blocks, and a second layer of parallel operation is carried out among a plurality of sampling points on each terrain block.
With the radar deployment position at (88.5, 39.1, 940), the detection distance is 210 kilometers, 120 shielding angles at intervals of 3 degrees are calculated, the time consumption of the original serial calculation based on a CPU is about 4s, and the time consumption of the method is about 0.7 s.
The embodiment of the application provides a radar shielding angle calculation method, time consumed by reading files is effectively reduced through DEM file conversion, memory occupation is reduced through DEM file management, the calculation amount is reduced while the precision is guaranteed through shielding angle calculation, GPU parallel calculation is used, high-performance hardware resources are fully utilized, the time consumed by calculation is reduced, and the purpose of rapid calculation is achieved.
In summary, the embodiment of the application provides a method for calculating a radar shielding angle based on a DEM, which includes preprocessing a target DEM file, wherein the preprocessing includes binary conversion of a file format and addition of a two-dimensional index identifier; calculating a longitude and latitude range of a target radar detection range according to the target radar deployment position and the detection distance, and acquiring a target DEM file covered by the longitude and latitude range; determining a limit distance according to the detection distance and a line-of-sight distance, wherein the line-of-sight distance is determined according to the height of a target radar and the highest point of the terrain; sampling and taking points in the azimuth of 0-360 degrees at a set angle interval by taking the deployment position of the target radar as an origin, the limit distance as a radius and the included angle with the true north as an azimuth, and calculating the radar shielding angle in each azimuth; the calculation of the radar shielding angles in all directions is divided into two layers by adopting a GPU (graphics processing unit) to be executed in parallel, the first layer is arranged among a plurality of DEM terrain blocks in parallel, and the second layer is arranged among a plurality of sampling points on each DEM terrain block in parallel. The time consumption for reading the file is effectively reduced through DEM file conversion, the memory occupation is reduced through DEM file management, the calculation amount is reduced while the precision is ensured through the shielding angle calculation, meanwhile, the GPU is used for parallel calculation, the high-performance hardware resources are fully utilized, the time consumption for calculation is reduced, and the purpose of rapid calculation is achieved.
Based on the same technical concept, an embodiment of the present application further provides a radar shading angle calculation system based on DEM, as shown in fig. 4, the system includes:
the DEM file management module 401 is configured to perform preprocessing on the target DEM file, where the preprocessing includes binary conversion of a file format and addition of a two-dimensional index identifier.
And the topographic data determining module 402 is used for calculating the longitude and latitude range of the detection range of the target radar according to the deployment position and the detection distance of the target radar, and acquiring a target DEM file covered by the longitude and latitude range.
A distance determining module 403, configured to determine a limiting distance according to the detection distance and a line-of-sight distance, where the line-of-sight distance is determined according to a height of the target radar and a highest point of the terrain.
A radar shielding angle determining module 404, configured to sample and take points in an azimuth of 0 degree to 360 degrees at a set angle interval with the target radar deployment position as an origin, the limited distance as a radius, and a true north included angle as an azimuth, and calculate a radar shielding angle in each azimuth; the calculation of the radar shielding angles in all directions is divided into two layers by adopting a GPU (graphics processing unit) to be executed in parallel, the first layer is arranged among a plurality of DEM terrain blocks in parallel, and the second layer is arranged among a plurality of sampling points on each DEM terrain block in parallel.
In a possible implementation, the radar shielding angle determining module 404 is further configured to: obtaining an azimuth line segment according to the origin, the radius and the azimuth angle; sampling and taking points on the azimuth line segment at set distance intervals; calculating an included angle between a connecting line of the surface position of each sampling point and the position of the target radar and the ground plane; and determining the maximum value of the included angle of each sampling point as a radar shielding angle.
In one possible embodiment, the set distance interval of the sampling points is half of the accuracy of the topographic data, the set distance interval of the distant sampling points is twice the accuracy of the topographic data, and the topographic data is the latitude and longitude range.
In one possible embodiment, the system further comprises: and the sending module is used for sending the radar shielding angles on each direction from the GPU to the CPU, drawing the radar shielding angles in a radar map mode, and releasing the memory occupied by the target DEM file.
In one possible implementation, the information of the target DEM file includes DEM file storage location, start longitude, start latitude, end longitude, end latitude, highest point elevation, and azimuth line segment sampling point.
In the present specification, each embodiment of the method is described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. Reference is made to the description of the method embodiments.
It is noted that while the operations of the methods of the present invention are depicted in the drawings in a particular order, this is not a requirement or suggestion that the operations must be performed in this particular order or that all of the illustrated operations must be performed to achieve desirable results. Additionally or alternatively, certain steps may be omitted, multiple steps combined into one step execution, and/or one step broken down into multiple step executions.
Although the present application provides method steps as in embodiments or flowcharts, additional or fewer steps may be included based on conventional or non-inventive approaches. The order of steps recited in the embodiments is merely one manner of performing the steps in a multitude of orders and does not represent the only order of execution. When an apparatus or client product in practice executes, it may execute sequentially or in parallel (e.g., in a parallel processor or multithreaded processing environment, or even in a distributed data processing environment) according to the embodiments or methods shown in the figures. The terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, the presence of additional identical or equivalent elements in a process, method, article, or apparatus that comprises the recited elements is not excluded.
The units, devices, modules, etc. set forth in the above embodiments may be implemented by a computer chip or an entity, or by a product with certain functions. For convenience of description, the above devices are described as being divided into various modules by functions, and are described separately. Of course, in implementing the present application, the functions of each module may be implemented in one or more software and/or hardware, or a module implementing the same function may be implemented by a combination of a plurality of sub-modules or sub-units, and the like. The above-described embodiments of the apparatus are merely illustrative, and for example, the division of the units is only one logical division, and other divisions may be realized in practice, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.
Those skilled in the art will also appreciate that, in addition to implementing the controller as pure computer readable program code, the same functionality can be implemented by logically programming method steps such that the controller is in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers and the like. Such a controller may therefore be considered as a hardware component, and the means included therein for performing the various functions may also be considered as a structure within the hardware component. Or even means for performing the functions may be regarded as being both a software module for performing the method and a structure within a hardware component.
The application may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, classes, etc. that perform particular tasks or implement particular abstract data types. The application may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
From the above description of the embodiments, it is clear to those skilled in the art that the present application can be implemented by software plus necessary general hardware platform. Based on such understanding, the technical solutions of the present application may be embodied in the form of a software product, which may be stored in a storage medium, such as a ROM/RAM, a magnetic disk, an optical disk, or the like, and includes several instructions for enabling a computer device (which may be a personal computer, a mobile terminal, a server, or a network device) to execute the method according to the embodiments or some parts of the embodiments of the present application.
The embodiments in the present specification are described in a progressive manner, and the same or similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. The application is operational with numerous general purpose or special purpose computing system environments or configurations. For example: personal computers, server computers, hand-held or portable devices, tablet-type devices, multiprocessor systems, microprocessor-based systems, set top boxes, programmable electronic devices, network PCs, minicomputers, mainframe computers, distributed computing environments that include any of the above systems or devices, and the like.
The above-mentioned embodiments are further described in detail for the purpose of illustrating the invention, and it should be understood that the above-mentioned embodiments are only illustrative of the present invention and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements, etc. made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (10)
1. A radar shielding angle calculation method based on DEM is characterized by comprising the following steps:
preprocessing a target DEM file, wherein the preprocessing comprises binary conversion of a file format and addition of a two-dimensional index identifier;
calculating a longitude and latitude range of a target radar detection range according to the target radar deployment position and the detection distance, and acquiring a target DEM file covered by the longitude and latitude range;
determining a limit distance according to the detection distance and a line-of-sight distance, wherein the line-of-sight distance is determined according to the height of a target radar and the highest point of the terrain;
sampling and taking points in the azimuth of 0-360 degrees at a set angle interval by taking the deployment position of the target radar as an origin, the limit distance as a radius and the included angle with the true north as an azimuth, and calculating the radar shielding angle in each azimuth;
the calculation of the radar shielding angles in all directions is divided into two layers by adopting a GPU (graphics processing unit) to be executed in parallel, the first layer is arranged among a plurality of DEM terrain blocks in parallel, and the second layer is arranged among a plurality of sampling points on each DEM terrain block in parallel.
2. The method of claim 1, wherein the calculating the radar shading angle in each azimuth by taking the target radar deployment location as an origin, the limit distance as a radius, and the true north as an azimuth, and sampling points at 0 to 360 degrees of azimuth at set angular intervals comprises:
obtaining an azimuth line segment according to the origin, the radius and the azimuth angle;
sampling and taking points on the azimuth line segment at set distance intervals;
calculating an included angle between a connecting line of the surface position of each sampling point and the position of the target radar and the ground plane;
and determining the maximum value of the included angle of each sampling point as a radar shielding angle.
3. The method of claim 2, wherein the set distance intervals of the sampling points are half the accuracy of the topographical data, the set distance intervals of the distant sampling points are twice the accuracy of the topographical data, and the topographical data is latitude and longitude.
4. The method of claim 1, wherein the method further comprises:
and sending the radar shielding angles in each direction from the GPU to the CPU, drawing the radar shielding angles in a radar chart mode, and releasing the memory occupied by the target DEM file.
5. The method of claim 1, wherein the information for the target DEM file includes DEM file storage location, start longitude, start latitude, end longitude, end latitude, highest point elevation, and azimuth line segment sampling points.
6. A DEM-based radar shading angle calculation system, the system comprising:
the DEM file management module is used for preprocessing a target DEM file, wherein the preprocessing comprises binary conversion of a file format and addition of a two-dimensional index identifier;
the topographic data determining module is used for calculating a longitude and latitude range of a target radar detection range according to the target radar deployment position and the detection distance, and acquiring a target DEM file covered by the longitude and latitude range;
the distance determining module is used for determining a limiting distance according to the detection distance and a line-of-sight distance, wherein the line-of-sight distance is determined according to the height of a target radar and the highest point of a terrain;
the radar shielding angle determining module is used for sampling and taking points in the directions of 0-360 degrees at set angle intervals by taking the deployment position of the target radar as an origin, the limit distance as a radius and a true north as an azimuth angle, and calculating the radar shielding angle in each direction; the calculation of the radar shielding angles in all directions is divided into two layers by adopting a GPU (graphics processing unit) to be executed in parallel, the first layer is arranged among a plurality of DEM terrain blocks in parallel, and the second layer is arranged among a plurality of sampling points on each DEM terrain block in parallel.
7. The system of claim 6, wherein the radar-obscuring angle determination module is further to:
obtaining an azimuth line segment according to the origin, the radius and the azimuth angle;
sampling and taking points on the azimuth line segment at set distance intervals;
calculating an included angle between a connecting line of the surface position of each sampling point and the position of the target radar and the ground plane;
and determining the maximum value of the included angle of each sampling point as a radar shielding angle.
8. The system of claim 7, wherein the set distance intervals for the sampling points are half the accuracy of the topographical data, and the set distance intervals for the distant sampling points are twice the accuracy of the topographical data, the topographical data being latitude and longitude.
9. The system of claim 6, wherein the system further comprises:
and the sending module is used for sending the radar shielding angles on each direction from the GPU to the CPU, drawing the radar shielding angles in a radar map mode, and releasing the memory occupied by the target DEM file.
10. The system of claim 6, wherein the information for the target DEM file comprises a DEM file storage location, a start longitude, a start latitude, a stop longitude, a stop latitude, a highest point elevation, and a azimuth line segment sampling point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010156881.1A CN111475916B (en) | 2020-03-09 | 2020-03-09 | DEM-based radar shielding angle calculation method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010156881.1A CN111475916B (en) | 2020-03-09 | 2020-03-09 | DEM-based radar shielding angle calculation method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111475916A true CN111475916A (en) | 2020-07-31 |
| CN111475916B CN111475916B (en) | 2023-08-01 |
Family
ID=71748177
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010156881.1A Active CN111475916B (en) | 2020-03-09 | 2020-03-09 | DEM-based radar shielding angle calculation method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111475916B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113093112A (en) * | 2021-03-23 | 2021-07-09 | 中国电子科技集团公司第二十九研究所 | Reconnaissance envelope calculation method and device for ground-oriented equipment under obstacle shielding and storage medium |
| CN113569446A (en) * | 2021-07-02 | 2021-10-29 | 上海机电工程研究所 | Air defense radar optimal deployment method and system based on three-dimensional space grid discretization |
| CN113589246A (en) * | 2021-05-24 | 2021-11-02 | 中国南方电网有限责任公司超高压输电公司广州局 | Target detection method and device based on radar and computer equipment |
| CN114509042A (en) * | 2020-11-17 | 2022-05-17 | 易图通科技（北京）有限公司 | Shielding detection method, shielding detection method of observation route and electronic equipment |
| CN115032605A (en) * | 2022-08-10 | 2022-09-09 | 中国航天科工集团八五一一研究所 | Radar target system analysis method based on target capability portrait |
| CN118112523A (en) * | 2023-12-18 | 2024-05-31 | 北京四象爱数科技有限公司 | Visualization method, device and storage medium for radar detection capability prediction |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109118583A (en) * | 2018-08-23 | 2019-01-01 | 中国科学院电子学研究所苏州研究院 | A kind of high-speed parallel terrain masking calculation method mixed based on CPU and GPU |
| CN109581349A (en) * | 2018-11-08 | 2019-04-05 | 中国电子科技集团公司第三十八研究所 | A kind of analysis of radar terrain masking and display methods |
-
2020
- 2020-03-09 CN CN202010156881.1A patent/CN111475916B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109118583A (en) * | 2018-08-23 | 2019-01-01 | 中国科学院电子学研究所苏州研究院 | A kind of high-speed parallel terrain masking calculation method mixed based on CPU and GPU |
| CN109581349A (en) * | 2018-11-08 | 2019-04-05 | 中国电子科技集团公司第三十八研究所 | A kind of analysis of radar terrain masking and display methods |
Non-Patent Citations (2)
| Title |
|---|
| 栾厚斌 等: ""多GPU并行技术在雷达信号处理中的应用"", 《电子技术与软件工程》 * |
| 陈超凡等: ""一种用于雷达遮蔽分析的DEM数据处理方法"", 《空军预警学院学报》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114509042A (en) * | 2020-11-17 | 2022-05-17 | 易图通科技（北京）有限公司 | Shielding detection method, shielding detection method of observation route and electronic equipment |
| CN114509042B (en) * | 2020-11-17 | 2024-05-24 | 易图通科技（北京）有限公司 | Shading detection method, shading detection method of observation route and electronic equipment |
| CN113093112A (en) * | 2021-03-23 | 2021-07-09 | 中国电子科技集团公司第二十九研究所 | Reconnaissance envelope calculation method and device for ground-oriented equipment under obstacle shielding and storage medium |
| CN113589246A (en) * | 2021-05-24 | 2021-11-02 | 中国南方电网有限责任公司超高压输电公司广州局 | Target detection method and device based on radar and computer equipment |
| CN113589246B (en) * | 2021-05-24 | 2023-09-01 | 中国南方电网有限责任公司超高压输电公司广州局海口分局 | Target detection method and device based on radar and computer equipment |
| CN113569446A (en) * | 2021-07-02 | 2021-10-29 | 上海机电工程研究所 | Air defense radar optimal deployment method and system based on three-dimensional space grid discretization |
| CN113569446B (en) * | 2021-07-02 | 2023-08-18 | 上海机电工程研究所 | Air defense radar optimal deployment method and system based on three-dimensional space grid discretization |
| CN115032605A (en) * | 2022-08-10 | 2022-09-09 | 中国航天科工集团八五一一研究所 | Radar target system analysis method based on target capability portrait |
| CN118112523A (en) * | 2023-12-18 | 2024-05-31 | 北京四象爱数科技有限公司 | Visualization method, device and storage medium for radar detection capability prediction |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111475916B (en) | 2023-08-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111475916A (en) | DEM-based radar shielding angle calculation method and system | |
| US11526536B2 (en) | Systems and methods for identifying grids of geographical region in map | |
| US11226431B2 (en) | Method and device for filling invalid regions of terrain elevation model data | |
| CN111398918B (en) | Radar detection capability analysis method under complex mountain environment | |
| Andrade et al. | Efficient viewshed computation on terrain in external memory | |
| US20170314947A1 (en) | Computing method for ridesharing paths, computing apparatus and recording medium using the same | |
| CN113190515A (en) | Heterogeneous parallel computing-based urban mass point cloud coordinate transformation method | |
| Ferreira et al. | A parallel algorithm for viewshed computation on grid terrains | |
| CN114549671A (en) | Grid coding method and computer system | |
| CN113936046A (en) | Object positioning method and device, electronic equipment and computer readable medium | |
| CN114155391A (en) | Trajectory processing method and device and computer equipment | |
| CN116521676A (en) | Map data management method, device, terminal equipment and readable storage medium | |
| US11002515B2 (en) | Intelligent artillery fire supporting device and operation method thereof | |
| Guan et al. | Fast approximate viewshed analysis based on the regular-grid digital elevation model: X-type partition proximity-direction-elevation spatial reference line algorithm | |
| CN115511701A (en) | Method and device for converting geographic information | |
| CN113421190B (en) | Super-resolution method, device, computer equipment and storage medium | |
| CN115205475A (en) | Method and device for determining terrain complexity, terminal equipment and storage medium | |
| Vo et al. | A parallel algorithm for local point density index computation of large point clouds | |
| CN113869715A (en) | Site grouping method, home site determining method, device and computer equipment | |
| CN113724229A (en) | Height difference determination method and device and electronic equipment | |
| CN109803268B (en) | Method and device for determining associated sector and computer readable storage medium | |
| CN118153915B (en) | Satellite maximum coverage calculating method, equipment and medium based on ellipsoid model | |
| CN117870651B (en) | Map high-precision acquisition method, memory and storage medium based on RTK-SLAM technology | |
| CN115880276B (en) | Tower pole operation state assessment method and system based on multi-period point cloud comparison | |
| CN115529601A (en) | Data processing method and device based on wireless propagation model and electronic equipment |
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 |