CN109581349B - Radar terrain shading analysis and display method - Google Patents

Radar terrain shading analysis and display method Download PDF

Info

Publication number
CN109581349B
CN109581349B CN201811325201.3A CN201811325201A CN109581349B CN 109581349 B CN109581349 B CN 109581349B CN 201811325201 A CN201811325201 A CN 201811325201A CN 109581349 B CN109581349 B CN 109581349B
Authority
CN
China
Prior art keywords
radar
shielding
distance
irng
angle
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811325201.3A
Other languages
Chinese (zh)
Other versions
CN109581349A (en
Inventor
周慧德
李帅
钮俊清
王顺利
安良
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CETC 38 Research Institute
Original Assignee
CETC 38 Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CETC 38 Research Institute filed Critical CETC 38 Research Institute
Priority to CN201811325201.3A priority Critical patent/CN109581349B/en
Publication of CN109581349A publication Critical patent/CN109581349A/en
Application granted granted Critical
Publication of CN109581349B publication Critical patent/CN109581349B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06—Systems determining position data of a target
    • G01S13/42—Simultaneous measurement of distance and other co-ordinates
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing

Abstract

The invention discloses a radar terrain shielding analysis and display method, which loads elevation data by taking a radar site as a center and taking a longitude and latitude range as a radius according to the longitude, the latitude and the height of the radar site set by a user and the longitude and latitude range to be calculated; initializing radar shielding data according to a radar detection minimum elevation angle, a maximum detection distance, azimuth precision and distance precision set by a user; analyzing the radar terrain shielding condition according to the loaded elevation data; and constructing a multi-dimensional display mode by using the calculated radar terrain shading data. The space is divided in two dimensions of distance and direction, and the terrain shielding information of each unit is respectively analyzed; the calculation method of the shielding angle adopts an empirical formula which takes various factors such as atmospheric refraction and the like into consideration; calculating the visibility range of the radar to a target with fixed height on the basis of calculating the shielding angle; the radar terrain shielding condition is observed and analyzed from multiple directions and angles by the aid of perspective range display and shielding angle curve display.

Description

Radar terrain shading analysis and display method
Technical Field
The invention relates to a display method of radar terrain, in particular to a radar terrain shielding analysis and display method.
Background
In the analytical research of radar, the selection of a place is an important part. The detection power of a radar is not only determined by the performance of the radar itself, but also by the terrain surrounding the radar site.
The traditional radar position selection is generally to survey according to the field, directly combine the landform and the landform, and then combine the technical indexes to carry out position screening. The traditional method has the defects of spending a large amount of time for surveying, occupying a large amount of manpower and material resources, having low precision and the like. Therefore, the problem to be solved is urgent when a proper radar position is selected quickly.
The radar terrain shading analysis and display method based on the Digital Elevation Model (DEM) overcomes the defects of the traditional method, does not limit the geographical position, and can obtain the part which can be seen by the radar in a certain area in an off-line manner. The specific application of the radar terrain shielding analysis and display method comprises multiple links, parameter setting, terrain shielding analysis and calculation and terrain shielding display. How to enable users to analyze and observe terrain shielding through various angles meets the requirements of users on continuous climbing, and the method is used for analyzing the current radar terrain shielding.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: how to provide a plurality of parameter setting interfaces and realize multi-dimensional display analysis results provides a radar terrain shading analysis and display method.
The invention solves the technical problems through the following technical scheme, and the invention comprises the following steps:
(1) according to the longitude, latitude and height of a radar site set by a user and a longitude and latitude range to be calculated, loading elevation data by taking the radar site as a center and taking the longitude and latitude range as a radius;
(2) initializing radar shielding data according to a radar detection minimum elevation angle, a maximum detection distance, azimuth precision and distance precision set by a user;
(3) analyzing the radar terrain shielding condition according to the loaded elevation data;
(4) and constructing a multi-dimensional display mode by using the calculated radar terrain shading data.
In the step (2), the mask angle target data includes:
elevation array dataEle [ N ]Azi][NRng]，
Height array dataAlt [ NAzi][NRng]，
Wherein
Figure BDA0001858552360000011
NAziAnd NRngRounding off, RmaxFor maximum detection distance, PAziOrientation accuracy, PRngIs distance accuracy.
And (3) performing coordinate transformation according to the longitude and latitude information of the loaded elevation data and the radar station address, and not performing shading angle calculation at the position exceeding the maximum detection distance.
The calculation process of the shielding angle in the step (3) is as follows:
(31) constructing shielding angle target data according to the azimuth precision and the distance precision set by a user, wherein the shielding angle target data comprises a shielding angle array and a height array;
(32) calculating the distance, the azimuth and the shielding angle of each grid of the elevation data relative to the radar site, and recalculating the shielding angle by using an empirical formula to obtain the numerical value of the shielding angle DemEle:
Figure BDA0001858552360000021
dataEle[iAzi][iRng]＝DemEle；
wherein the content of the first and second substances,
Figure BDA0001858552360000022
the iAzi and iRng are rounded, CorrectAll is the corrected height,
DemAzi is the azimuth, DemRng is the distance, and RadarAlt is the height;
(33) and calculating the visibility range of the radar for the target with the given height according to the shielding angle data.
The step (33) comprises the steps of:
(331) using the mask angle array, calculate the height array:
in the same direction, the distance is searched from small to large according to the current distance sequence iRng (2 is more than or equal to iRng is less than or equal to NRng) Maximum value of inner shielding angle, Elemax＝max{dataEle[iAzi][1],dataEle[iAzi][2]…dataEle[iAzi][iRng-1]}，
If dataEle [ iAzi][iRng]>ElemaxThen dataEle [ iAzi ]][iRng]＝Elemax；
Calculating the height by using an empirical formula to obtain a height array:
Alt＝Rng*sin(dataEle[iAzi][iRng])+Rng2/(2*8.5*e6)+RadarAlt+CorrectAlt，
wherein Rng ═ (iRng +0.5) × PRng，dataAlt[iAzi][iRng]＝Alt；
(332) Acquiring a panoramic array of preset target point heights:
the following equation is calculated at each orientation iAzi,
Figure BDA0001858552360000023
obtaining the iRng of the reaction solution,
selection of communication arrayRng[iAzi]＝iRng*PRng。
In the step (4), the display mode is an optional display mode: polar coordinate perspective range display and azimuth-shielding angle two-dimensional rectangular coordinate shielding angle curve display.
In the step (4), the display mode is a through-view range display: according to the target height selected by a user, the visible range of the radar to the target with the designated target height is displayed, and the selection of a plurality of designated target heights is supported.
In the step (4), the display mode is a shading angle curve display: according to the target distance selected by a user, the shielding condition of the radar on the designated distance is displayed, data of a mouse pickup curve are supported, and an azimuth corresponding to the position of the mouse, the shielding angle of the azimuth and the shielding distance causing the shielding are displayed.
Compared with the prior art, the invention has the following advantages: according to the invention, a more flexible calculation parameter setting interface is adopted, so that a user can set the radar station address, the latitude and longitude calculation range, the radar detection minimum elevation angle, the maximum detection distance, the azimuth precision and the distance precision according to actual requirements, and meanwhile, the loading and calculation of redundant data are avoided, and the calculation efficiency is improved; the space is divided in two dimensions of distance and direction, and the terrain shielding information of each unit is respectively analyzed; the calculation method of the shielding angle adopts an empirical formula considering various factors such as atmospheric refraction and the like, so that the result is more accurate; calculating the visibility range of the radar to a target with fixed height on the basis of calculating the shielding angle; two display methods are adopted, the through-view range display and the shielding angle curve display provide a way for a user to observe and analyze the radar terrain shielding condition from multiple directions and angles.
Drawings
FIG. 1 is a flow chart of the operation of the present invention;
FIG. 2 is a flow chart of the present invention calculation;
FIG. 3 is a flowchart of the operation of the through view range display;
FIG. 4 is an effect diagram of a through view range display;
FIG. 5 is a flowchart of the operation of the shading angle curve display;
fig. 6 is an effect diagram of a shading angle curve display.
Detailed Description
The following examples are given for the detailed implementation and specific operation of the present invention, but the scope of the present invention is not limited to the following examples.
As shown in fig. 1, the working flow of the method for calculating and displaying the shielding angle of the intelligence radar of the present embodiment includes the following processes:
1. the user sets calculation parameters including longitude (degree) RadarLat, latitude (degree) RadarLon, altitude (meter) RadarAlt of radar site, longitude and latitude range to be calculated and maximum detection distance R according to the requirementmaxAzimuth accuracy PAziDistance accuracy (m) PRng；
2. Initializing the shielding angle target data according to the set calculation parameters, and initializing coordinate transformation:
the shading angle target data includes:
elevation array dataEle [ N ]Azi][NRng]，
Height array dataAlt [ NAzi][NRng]，
Wherein
Figure BDA0001858552360000041
NAziAnd NRngGetting the whole;
using longitude, latitude and height of radar site to initialize coordinate transformation matrix;
3. loading required elevation data according to calculation parameters set by a user: longitude (degrees) DemLon, longitude (degrees) DemLat, altitude (meters) DemAlt;
4. and calculating a shielding angle by combining radar site and elevation data, wherein the specific flow is as shown in fig. 2:
4-1: calculating the distance, the azimuth and the shielding angle of each grid of the elevation data relative to the radar site:
distance (m) DemRng, azimuth (degree) DemAzi, shielding angle (degree) DemEle;
4-2: recalculating the shielding angle by using an empirical formula to obtain a shielding angle array:
Figure BDA0001858552360000042
dataEle[iAzi][iRng]＝DemEle；
wherein:
Figure BDA0001858552360000043
rounding iAzi and iRng, and CorrectAll is the corrected height;
4-3: and calculating a height array by using the shading angle array obtained in the step 4-2:
in the same direction, the distance is searched from small to large according to the current distance sequence iRng (2 is more than or equal to iRng is less than or equal to NRng) Maximum value of inner shielding angle, Elemax＝max{dataEle[iAzi][1],dataEle[iAzi][2]…dataEle[iAzi][iRng-1]If dataEle [ iAzi ]][iRng]>ElemaxThen dataEle [ iAzi ]][iRng]＝Elemax；
Calculating the height by using an empirical formula to obtain a height array:
Alt＝Rng*sin(dataEle[iAzi][iRng])+Rng2/(2*8.5*e6) + RadarAlt + corractalt, where Rng ═ iRng +0.5 × PRng，
dataAlt[iAzi][iRng]＝Alt；
5. Acquiring a panoramic array of preset target point heights:
the following equation is calculated at each orientation iAzi,
Figure BDA0001858552360000044
obtaining iRng, SelectRng[iAzi]＝iRng*PRng；
6. And (3) shielding angle through display: the shielding angle is used for displaying the visual range of the target with the designated height in a perspective way, the polar coordinate P display mode is adopted for displaying, and the azimuth-distance (array Select) is drawnRng[iAzi]The graph of (2) can be superposed with graphs corresponding to a plurality of targets with specified heights, the rendering effect of the display mode is shown in fig. 4, the specified heights can be switched in real time, the curve is updated, and the operation flow is shown in fig. 3.
7. Acquiring a specified distance shielding condition:
acquiring a shielding angle array of a specified distance: finding a designated distance Select from an elevation arrayRngShielding angle data Select ofele[iAzi]＝dataEle[iAzi][iRng]Wherein iRng ═ SelectRng/PRng，iRng≤NRngAnd rounding;
obtaining a distance array that results in a specified distance mask angle:
RngEle[iAzi]＝iRng*PRngwherein iRng is {0, SelectrngMax [ dataEle [ iAzi ] at a distance within the range from small to large][iRng]iRng corresponding to the obtained frequency spectrum;
8. shading angle curve shows: shading Angle Curve shows the azimuth-shading Angle (array Select) at a given distanceele[iAzi]) The graph of (2) can support the mouse to pick up the data of the curve, and display the azimuth corresponding to the position of the mouse, the shielding angle of the azimuth and the shielding distance causing the shielding, the rendering effect of the display mode is shown in fig. 6, and the specified distance can be switched in real time to update the curve, and the operation flow is shown in fig. 5.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents and improvements made within the spirit and principle of the present invention are intended to be included within the scope of the present invention.

Claims (6)

1. A radar terrain shading analysis and display method is characterized by comprising the following steps:
(1) according to the longitude, latitude and height of a radar site set by a user and a longitude and latitude range to be calculated, loading elevation data by taking the radar site as a center and taking the longitude and latitude range as a radius;
(2) initializing radar shielding data according to a radar detection minimum elevation angle, a maximum detection distance, azimuth precision and distance precision set by a user;
(3) analyzing the radar terrain shielding condition according to the loaded elevation data;
(4) constructing a multi-dimensional display mode by using the calculated radar terrain shading data;
in the step (2), the mask angle target data includes:
elevation array dataEle [ N ]Azi][NRng]，
Height array dataAlt [ NAzi][NRng]，
Wherein
Figure FDA0002409289540000011
NAziAnd NRngRounding off, RmaxFor maximum detection distance, PAziOrientation accuracy, PRngDistance accuracy;
the calculation process of the shielding angle in the step (3) is as follows:
(31) constructing shielding angle target data according to the azimuth precision and the distance precision set by a user, wherein the shielding angle target data comprises a shielding angle array and a height array;
(32) calculating the distance, the azimuth and the shielding angle of each grid of the elevation data relative to the radar site, and recalculating the shielding angle by using an empirical formula to obtain the numerical value of the shielding angle DemEle:
Figure FDA0002409289540000012
dataEle[iAzi][iRng]＝DemEle；
wherein the content of the first and second substances,
Figure FDA0002409289540000013
the iAzi and iRng are rounded, CorrectAll is the corrected height,
DemAzi is the azimuth, DemRng is the distance, and RadarAlt is the height;
(33) and calculating the visibility range of the radar for the target with the given height according to the shielding angle data.
2. The radar terrain shading analysis and display method according to claim 1, wherein in the step (3), coordinate transformation is performed according to longitude and latitude information of the loaded elevation data and a radar site, and shading angle calculation is not performed for a position beyond a maximum detection distance.
3. A radar terrain occlusion analysis and display method according to claim 1, characterized in that said step (33) comprises the steps of:
(331) using the mask angle array, calculate the height array:
in the same direction, the distance is searched from small to large according to the current distance sequence iRng (2 is more than or equal to iRng is less than or equal to NRng) Maximum value of inner shielding angle, Elemax＝max{dataEle[iAzi][1]，dataEle[iAzi][2]...dataEle[iAzi][iRng-1]}，
If dataEle [ iAzi][iRng]＞ElemaxThen dataEle [ iAzi ]][iRng]＝Elemax；
Calculating the height by using an empirical formula to obtain a height array:
Alt＝Rng*sin(dataEle[iAzi][iRng])+Rng2/(2*8.5*e6)+RadarAlt+CorrectAlt，
wherein Rng ═ (iRng +0.5) × PRng，dataAlt[iAzi][iRng]＝Alt；
(332) Acquiring a panoramic array of preset target point heights:
the following equation is calculated at each orientation iAzi,
Figure FDA0002409289540000021
obtaining the iRng of the reaction solution,
selection of communication arrayRng[iAzi]＝iRng*PRng。
4. The radar terrain shading analysis and display method of claim 1, wherein in the step (4), the display modes are selectable display modes: polar coordinate perspective range display and azimuth-shielding angle two-dimensional rectangular coordinate shielding angle curve display.
5. The radar terrain shading analysis and display method according to claim 1, wherein in the step (4), the display mode is a through-view range display: according to the target height selected by a user, the visible range of the radar to the target with the designated target height is displayed, and the selection of a plurality of designated target heights is supported.
6. The radar topography shading analyzing and displaying method according to claim 1, wherein in the step (4), the displaying manner is shading angle curve display: according to the target distance selected by a user, the shielding condition of the radar on the designated distance is displayed, data of a mouse pickup curve are supported, and an azimuth corresponding to the position of the mouse, the shielding angle of the azimuth and the shielding distance causing the shielding are displayed.
CN201811325201.3A 2018-11-08 2018-11-08 Radar terrain shading analysis and display method Active CN109581349B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811325201.3A CN109581349B (en) 2018-11-08 2018-11-08 Radar terrain shading analysis and display method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811325201.3A CN109581349B (en) 2018-11-08 2018-11-08 Radar terrain shading analysis and display method

Publications (2)

Publication Number Publication Date
CN109581349A CN109581349A (en) 2019-04-05
CN109581349B true CN109581349B (en) 2020-06-26

Family

ID=65921908

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811325201.3A Active CN109581349B (en) 2018-11-08 2018-11-08 Radar terrain shading analysis and display method

Country Status (1)

Country Link
CN (1) CN109581349B (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110456319A (en) * 2019-08-29 2019-11-15 西安电子工程研究所 A kind of radar intervisibility calculation method based on SRTM
CN111060882B (en) * 2020-03-06 2021-03-16 成都蓉奥科技有限公司 Rapid calculation method for radar terrain shielding blind area under multiple heights
CN111475916B (en) * 2020-03-09 2023-08-01 北京仿真中心 DEM-based radar shielding angle calculation method and system
CN111523240B (en) * 2020-04-28 2023-05-05 鲁东大学 View analysis method based on DEM data
CN114509042A (en) * 2020-11-17 2022-05-17 易图通科技（北京）有限公司 Shielding detection method, shielding detection method of observation route and electronic equipment
CN113093112B (en) * 2021-03-23 2023-05-23 中国电子科技集团公司第二十九研究所 Ground equipment-oriented reconnaissance envelope calculation method, device and storage medium under obstacle shielding
CN113156441B (en) * 2021-03-25 2022-05-27 广东省气象探测数据中心 Effective three-dimensional airspace subdivision approximation calculation method for meteorological radar detection
CN115032605A (en) * 2022-08-10 2022-09-09 中国航天科工集团八五一一研究所 Radar target system analysis method based on target capability portrait

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TW358247B (en) * 1997-11-06 1999-05-11 Chung Shan Inst Of Science Radar mask direction tolerance automatic testing method and device
CN103778477B (en) * 2013-11-20 2017-07-18 中国民用航空总局第二研究所 Monitor station cloth station method and apparatus
CN108537361A (en) * 2018-03-06 2018-09-14 西安大衡天成信息科技有限公司 A kind of addressing appraisal procedure of mobile radar and communication station

Also Published As

Publication number Publication date
CN109581349A (en) 2019-04-05

Similar Documents

Publication Publication Date Title
CN109581349B (en) Radar terrain shading analysis and display method
US8300920B2 (en) System and method for measuring form and position tolerances of an object
EP2343659B1 (en) Techniques for density mapping
US6173067B1 (en) System and method for rapid determination of visibility-based terrain properties over broad regions
Tabik et al. Simultaneous computation of total viewshed on large high resolution grids
Fresno et al. Survey on the performance of source localization algorithms
US8817093B2 (en) Photogrammetric networks for positional accuracy
CN107833279B (en) DEM-based terrain slope analysis method
US10380454B2 (en) Identifying a visual center of a polygon
WO2019100188A1 (en) Method for planning operation route of unmanned aerial vehicle and ground end device
CN110686704A (en) Pose calibration method, system and medium for laser radar and combined inertial navigation
KR102042773B1 (en) Method to select subset of tracking satellites for receiver of satellite navigational system, and receiver, system, and apparatus therewith
CN108763811A (en) Dynamic data drives forest fire appealing prediction technique
US7778808B2 (en) Geospatial modeling system providing data thinning of geospatial data points and related methods
CN104331927A (en) Method and system for simulating reservoir digital outcrop surface
Fawcett et al. Two decades of digital photogrammetry: Revisiting Chandler’s 1999 paper on “Effective application of automated digital photogrammetry for geomorphological research”–a synthesis
CN111709988B (en) Method and device for determining characteristic information of object, electronic equipment and storage medium
CN109033181B (en) Wind field geographic numerical simulation method for complex terrain area
KR101425425B1 (en) Apparatus for gridding radar data and method thereof
JP2006300602A (en) Acquiring method, acquiring device, and acquiring system for positioning data
US20230186594A1 (en) Information processing apparatus, information processing method, and non-transitory computer readable medium
Lohani et al. Generating LiDAR data in laboratory: LiDAR simulator
CN115690339A (en) High-precision airspace perspective analysis method in cylindrical coordinate space
JP2021149063A (en) Topographic map output device, topographic map output method, and program
JP2022537647A (en) Method and satellite access system for determining the likelihood of satellite-borne sensors accessing a target area

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