CN106932772A - A kind of radar coverage display methods by the influence of topography towards digital earth - Google Patents
A kind of radar coverage display methods by the influence of topography towards digital earth Download PDFInfo
- Publication number
- CN106932772A CN106932772A CN201710153349.2A CN201710153349A CN106932772A CN 106932772 A CN106932772 A CN 106932772A CN 201710153349 A CN201710153349 A CN 201710153349A CN 106932772 A CN106932772 A CN 106932772A
- Authority
- CN
- China
- Prior art keywords
- radar
- point
- points
- topography
- influence
- 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.)
- Pending
Links
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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a kind of radar coverage display methods by the influence of topography towards digital earth, the present invention had both simplified problem, had improve speed from the angle solve problem of geometric optics, and theoretically ensured the correctness for solving；And, it is separate between each sectional side elevation of radar coverage, therefore the present invention also uses parallel computing, concomitantly calculate each sectional side elevation is influenceed by landform, further increases calculating speed；Finally, the present invention copes with the situation that radar coverage is influenceed by many mountain peaks suitable for various complicated terrain environments.
Description
Technical field
The invention belongs to technical field, and in particular to a kind of radar coverage by the influence of topography towards digital earth
Display methods.
Background technology
Radar is that radar is used as the main unit for influenceing Battle Field Electromagnetic using the electronics of electromagnetic wave detection target
Element, is widely used in the aspects such as target detection, target following, target positioning, tachometric survey, target identification and precise guidance.
But, the electromagnetic information of radar is not seen, be can not touch, and how to show the direction that it makes joint efforts as scientific research personnel.And thunder
Can be influenceed by many environmental factors up to environment, such as landform, atural object, air, electromagnetic interference etc., present invention mainly solves receiving
The radar coverage three-dimensional visualization problem of the influence of topography.
With virtual reality technology extensive use militarily, extensive three-dimensional battlefield surroundings generation, based on ground
Manage the technologies such as the Virtual Battlefield battle state display and multidimensional information general performance of information system and be obtained for significant progress.And it is empty
It is the indispensable part in future digital battlefield to intend battlefield surroundings, and it is to landform, atural object, landforms, ocean, meteorology etc.
Relative maturity is showed, but in terms of the three-dimensional representation of electromagnetic environment, is studied also very weak.How by the effect of radar with shape
The mode of elephant is presented to commander and trainer, and how to represent radar is influenceed by environment such as landform, ocean, meteorologies, is all turned into
The unavoidable problem of future digital battlefield surroundings.
Traditional radar electromagnetic information representation, mainly after the related data for obtaining electromagnetic information, by manual plotting
Or area of computer aided is showed radar coverage-diagram with two dimensional form, to realize the performance of Electromagnetic Situation.The above method is deposited
Substantially not enough, not only limited accuracy, and also plane curve is not directly perceived enough, and when target component changes, radar coverage-diagram can not
Change therewith.To obtain accurate radar coverage, it is necessary to solve the three-dimensional vector ripple as obtained by Maxwell equation is derived
Dynamic equation, due to complexity of the landform on its three component directions and the polytropism of air, therefore solving speed is very slow.
The content of the invention
In view of this, it is an object of the invention to provide a kind of radar coverage by the influence of topography towards digital earth
Display methods, have studied sphere of action of the radar under the influence of by terrain environment, according to scientific algorithm from the angle of geometric optics
Visual thought intuitively shows in virtual battlefield environment the three-dimensional radar volume.
A kind of radar coverage display methods by the influence of topography towards digital earth, comprises the following steps：
Step 1, radar coverage is sampled, specially：
In the range of radar directional pattern, by azimuth direction by setting step-lengthPieces is divided into cut
Piece；Then on the center of radar directional pattern, by each section of azimuth direction on pitching angular direction by setting
Fixed step size Δ θ is divided into Layers section, section one radar lobe of correspondence of azimuth direction；
Step 2, for each radar lobe, on lobe border by setting step-length sampling, obtain sampled point；Judge radar
Each sampled point and the position relationship of corresponding topographic(al) point, determine whether sampled point will be adapted to ground according to position relationship on lobe border
On form point, detailed process is：
S20, first from radar center point, along radar lobe the latter half, counterclockwise, find first
By the sampled point of terrain shading；Calculate the height Ht of the corresponding topographic(al) point of radar center point R, and central point height value Hr；
If Hr>Ht, continues more next sampled point and the height relationships of corresponding topographic(al) point, is less than until finding first height value
The accordingly sampled point of form point height value, then the sampled point is the sampled point by the influence of topography, is defined as P points, correspondence topographic(al) point
Q points are defined as, this is adapted at topographic(al) point Q by the sampled point of the influence of topography；
S21, the sampled point for continuing to judge behind on lobe border according to the method for S20, it is next by landform until finding
The sampled point of influence, is defined as current sampled point P1 to be modified；Connection RQ simultaneously extends, the friendship on extended line and radar lobe border
Point is defined as Q '；Use the length ‖ QQ ' ‖ of line segment QQ ' divided by the number of sampled point between radar lobe border PQ ', obtain step delta
r1；From Q points, along QQ ' directions advance Δ r1, the point for being reached is defined as K1, and the corresponding topographic(al) points of K1 are Q1；Connection
RQ1, and extend, extended line is defined as sampled point Q1 ' with the intersection point on radar lobe border, then judges：
If 1., ‖ RQ1 ‖<‖ RQ1 ' ‖, then topographic(al) point Q1 P1 is adapted at Q1 still in radar maximum detection range；
If 2., ‖ RQ1 ‖ >=‖ RQ1 ' ‖, perform S22；
1. and 2. judged according to situation for all sampled points between PQ ', until all sampled points between completion PQ '
Judgement, finally perform S22；
S22, hypothesis P1 are previous processed sampled point, and are corrected at Q1；Use the length ‖ of line segment Q1Q1 '
Q1Q1 ' ‖ obtain step delta r divided by the number of sampled point between radar lobe border Q1Q1 '；From Q1 points, along Q1Q1 '
The point that direction advance Δ r is reached is defined as K points, judges：
If A, K point are higher than corresponding topographic(al) point：Next sampled point P2 of sampled point P1 is adjusted at K points, along
Q1Q1 ' the directions Δ r that readvances reaches new K points, subsequent sampling point is adjusted at new K points, then judged：
If new K points reach radar lobe boundary point Q1 ', analysis terminates, and performs step 3；
If new K points do not reach radar lobe boundary point Q1 ', execution S22 is returned to；
If B, K point are lower than corresponding topographic(al) point, S23 is performed；
S23, hypothesis P1 are previous processed sampled point, and are corrected at Q1；From Q1 points, along Q1Q1 '
Direction advance Δ r reaches K points；Judge：If K points are lower than corresponding topographic(al) point Q2, but extended line and the radar border of RQ2 friendship
In point Q2 ', and ‖ RQ2 ‖>‖ RQ2 ' ‖, then be adapted to Q2 ' places by P2 points, and analysis terminates, and performs step 3；If ‖ RQ2 ‖≤‖
RQ2 ' ‖, P2 is adapted at Q2, and the Δ r that moves on reaches another new point, is returned and is performed S23；
Step 3, according to the revised sampled point of step 2, new radar lobe border is drawn, so as to obtain under the influence of topography
Radar coverage.
Preferably, each radar lobe is processed using the method for step 2 simultaneously, after obtaining result, then using step 3
Method draw radar bounds.
Preferably, expression way of the borderline sampled point of radar lobe using spherical coordinates, asθ is
The angle of pitch,It is azimuthRepresent at pitching angle theta and azimuthThe corresponding detection range of upper radar.
Preferably, detection range
In formula, PtIt is the power of transmission signal；τ is pulse bandwidth；GtIt is transmitting antenna power gain；GrIt is reception antenna
Power gain；σ is cross section of radar targets；λ is wavelength；FtIt is the directional diagram propagation factor from transmitting antenna to target；FrBe from
Directional diagram propagation factor of the target to reception antenna；K is Boltzmann constant；TsIt is reception system noise temperature；D0For detection because
Son；CBIt is bandwidth correction factor；L is the system loss factor.
Preferably, detection rangeIt is reduced to following form：
Wherein, RmaxIt is maximum detectable range of the radar on antenna gain direction in free space.
Preferably, in the step 2, by the sample point coordinate in radar lobe from asking coordinate system to be converted to geodetic coordinates
System, then judgement is compared with the coordinate under earth coordinates of topographic(al) point, radar side is finally drawn under earth coordinates
Boundary's scope.
Preferably, the step-lengthSpan be 0.1 degree to 10 degree.
Preferably, the span of the step delta θ is 0.1 degree to 10 degree.
The present invention has the advantages that：
The method of traditional calculating influence of topography radar coverage needs to solve the ripple derived by Maxwell equation
Dynamic equation, not only process is complicated, and speed is slow.Theoretical according to duality principle, electromagnetic wave can be regarded as the expansion of light wave
Exhibition.Therefore, the present invention had both simplified problem, had improve speed, and theoretically ensured from the angle solve problem of geometric optics
The correctness of solution.And, it is separate between each sectional side elevation of radar coverage, therefore the present invention is also using simultaneously
Row computing technique, concomitantly calculate each sectional side elevation is influenceed by landform, further increases calculating speed.Finally, the present invention is suitable
For various complicated terrain environments, the situation that radar coverage is influenceed by many mountain peaks is coped with.
Brief description of the drawings
Fig. 1 represents investigative range schematic diagram of the radar in free space；
Fig. 2 represents the sectional side elevation of radar coverage under terrain shading；
Fig. 3 represents the schematic diagram that radar beam is blocked by landform part；
Fig. 4 represents radar beam by the completely cut off schematic diagram of landform；
Fig. 5 represents the first schematic diagram by influence of topography sampled point of amendment；
Fig. 6 represents amendment subsequent sampling point schematic diagram；
Fig. 7 represents that subsequent sampling point schematic diagram is corrected in continuation；
Fig. 8 represents amendment subsequent sampling point schematic diagram；
Fig. 9 represents the flow chart of influence of topography radar coverage algorithm；
Figure 10 is represented by the radar coverage schematic diagram under the influence of topography.
Specific embodiment
Develop simultaneously embodiment below in conjunction with the accompanying drawings, and the present invention will be described in detail.
The investigative range of radar is a solid, irregular occluding surface body, its spatial form by radar antenna
Pattern function is determined.Thought based on discrete sampling, can carry out discretization to radar space exploration range boundary.Radar is visited
The discretization for surveying scope is carried out under spherical coordinates, at its pitching angle theta and azimuthIt is upper to carry out discretization respectively, will be whole
Continuous area of space is divided into many space vectorsThen these vector institutes are calculated by pattern function right
The detection range answeredThe boundary point for then detecting is
Because the antenna radiation pattern function of radar is given under spherical coordinates, the angle of pitch of θ correspondence radars,Correspondence thunder
The azimuth for reaching, R correspondence pole spans.
Assuming that azimuthal scope of the antenna radiation pattern of radar is 0 °~360 °, the scope of the angle of pitch is 0 °~90 °, is
By radar coverage boundary discrete method, azimuthal and the angle of pitch are sampled respectively.
Sample at azimuth:By the way of similar graticules, in the range of radar directional pattern, by azimuth direction
Pieces section is divided into, each one azimuth of section correspondenceThe angle of cut between contiguous slices is azimuth sampling
Step-lengthStep-lengthSpan be 0.1 degree to 10 degree.
The angle of pitch is sampled:On the center of radar directional pattern, each is cut into slices in the form of direction vector again
Secondary subdivision, pitching angular direction is divided into Layers section, each direction vector one azimuth angle theta of correspondencem, between adjacent vector
Pitching angular difference is angle of pitch sampling step length Δ θ.The span of step delta θ is 0.1 degree to 10 degree.
One of azimuth direction section one lobe of correspondence, all lobes be stitched together be exactly radar space exploration model
Enclose.
Because the resolution ratio of terrain data is not quite similar, adopting for radar coverage can be dynamically adjusted according to actual needs
Sample rate, makes itself and terrain match.Compared to traditional two dimensional surface display mode, the present invention display radar coverage more
Intuitively, more preferably, the three-dimensional visible effect of investigative range of the radar under free space is as shown in Figure 1 for display effect.
Electromagnetic wave can be produced reflex in spatial by terrain shading, and the back side of landform occurs a range of
Detection blind area, has various methods to calculate the influence that terrain shading brings, such as geometrical optics approach, electromagnetic field diffraction method.
Electromagnetic wave propagation is simulated using geometrical optics approach, new border will be adjusted to by the radar boundary point of the influence of topography
On, draw out by the actual radar coverage under terrain shading.Fig. 2 gives radar detection model in the case of consideration terrain shading
The sectional side elevation for enclosing.
Assuming that radar is at R, elliptical arc section is by one sphere of action in vertical section of radar before the influence of topography.Thunder
Between terrain shading is subject to since the A points, it is necessary to impacted radar boundary point is adjusted into hillside AB；After crossing mountain peak B,
Due to being blocked by AB, electromagnetic wave cannot be propagated in region below BC, and radar can only detect the region of more than BC, it is necessary to will
Affected radar boundary point is adjusted on BC；After reaching C points, electromagnetic wave runs into hillside CD higher, impacted radar
Between boundary point is also required to be adjusted to CD；After crossing mountain peak D again, the detection border E of radar is reached.The closure that RABCDER is constituted
Region is the investigative range after current lobe is subject to the influence of topography.
There are two kinds of possible situations under the influence of terrain shading in radar.
Situation 1：Radar center is located at landform top, and radar beam is partly blocked, as shown in figure 3, shadow region is
Actual radar coverage.
Situation 2：Radar center is located at landform top, and radar beam is completely interrupted, as shown in figure 4, shadow region is
Actual radar coverage.
Situation of the radar lobe under unimodal blocking is discussed above, the influence of multimodal can be sentenced by combinations thereof
It is disconnected, it is possible thereby to realize radar coverage repeatedly blocked by landform after display effect.Thunder under the influence of topography is first illustrated below
Up to the specific method of investigative range amendment, the flow chart of whole algorithm is then provided, as shown in Figure 9.
First from radar center point, along radar lobe the latter half, counterclockwise, first is found by ground
The sampled point that shape is blocked.Wherein, topographic(al) point refers to the sampled point on actual landform.With each sample in radar lobe in
Topographic(al) point on same vertical curve is the corresponding topographic(al) point of sampled point.Calculate the height of the corresponding topographic(al) points of radar center point R
Ht, and central point height value Hr.If Hr>Ht, continues the height relationships of more next sampled point and corresponding topographic(al) point,
Until finding first sampled point height value less than accordingly form point height value, then the sampled point is to be sampled by the influence of topography
Point, is defined as P points, and correspondence topographic(al) point is defined as Q points, as shown in figure 5, this is adapted into topographic(al) point by the sampled point of the influence of topography
At Q, i.e., replace the coordinate value of correspondence sampled point with the coordinate value of topographic(al) point Q.
To follow-up sampled point, it is modified according to different situations：
(1) continue to judge sampled point below according to the method described above, until finding next sampled point by the influence of topography,
Current sampled point P1 to be modified is defined as, as shown in fig. 6, connect RQ and extend, the intersection point on extended line and radar lobe border
It is defined as Q '.Use the length ‖ QQ ' ‖ of line segment QQ ' divided by the number of sampled point between radar lobe border PQ ', obtain step delta
r1.From Q points, along QQ ' directions advance Δ r1, the point for being reached is defined as K1, and the corresponding topographic(al) points of K1 are Q1.Connection
RQ1, and extend, extended line is defined as sampled point Q1 ' with the intersection point on radar lobe border, then judges：
If 1., ‖ RQ1 ‖<‖ RQ1 ' ‖, then topographic(al) point Q1 P1 is adapted at Q1 still in radar maximum detection range；
If 2., ‖ RQ1 ‖ >=‖ RQ1 ' ‖, perform (2)；
1. and 2. judged according to situation for all sampled points between PQ ', until all sampled points between completion PQ '
Judgement, finally perform (2)；
(2) as shown in Figure 7, it is assumed that P1 is previous processed sampled point, and it is corrected at Q1.Use line segment Q1Q1 '
Length ‖ Q1Q1 ' ‖ divided by the number of sampled point between radar lobe border Q1Q1 ', obtain step delta r；From Q1 points, edge
That Q1Q1 ' directions advance Δ r arrival is defined as K points.If K points are higher than corresponding topographic(al) point, illustrate to have passed over previous
Individual mountain peak, next sampled point P2 of sampled point P1 is adjusted at K points.New K is reached along Q1Q1 ' the directions Δ r that readvances
Point, subsequent sampling point is adjusted at new K points, is judged：If new K points reach radar lobe boundary point Q1 ', analysis knot
Beam；If new K points do not reach radar lobe boundary point Q1 ', execution (2) is returned to；If K points are than corresponding topographic(al) point
It is low, then perform (3).
(3) as shown in Figure 8, it is assumed that P1 is previous processed sampled point, and it is corrected at Q1.From Q1 points,
K points are reached along Q1Q1 ' directions advance Δ r.Judge：If although K points are lower than corresponding topographic(al) point Q2, the extension of RQ2
Line meets at point Q2 ', and ‖ RQ2 ‖ with radar border>‖ RQ2 ' ‖, this explanation Q2 alreadys exceed the maximum range of radar,
P2 points are adapted to Q2 ' places, analysis terminates；If ‖ RQ2 ‖≤‖ RQ2 ' ‖, P2 is adapted at Q2, the Δ r that moves on reaches another
One K point, returns and performs (3).
Be given above be radar single lobe up-sampling point modification method, according to azimuth coverage and sampling step length according to
After sampled point in the secondary all lobes for the treatment of, new radar border is drawn, so as to obtain the radar coverage under the influence of topography.
Because separate between each lobe of azimuth direction sampling, using parallel computing,
Concurrently calculate each lobe is influenceed by landform.After using parallel computation, successfully solve that conventional method calculating speed is slow asks
Topic, hence it is evident that shorten SRT.Using the three-dimensional of the radar coverage under the influence of topography that the method for the present invention is realized
Shown in design sketch 10.
Radar it is determined that pitching angle theta, azimuthOn detection range computing formula be：
In formula, PtIt is the power (antenna end) of transmission signal；τ is pulse bandwidth；GtIt is transmitting antenna power gain；GrFor
Reception antenna power gain；σ is cross section of radar targets；λ is wavelength；FtBe directional diagram from transmitting antenna to target propagate because
Son；FrIt is the directional diagram propagation factor from target to reception antenna；K is Boltzmann constant；Ts is reception system noise temperature；
D0 is detecting factor；CB is bandwidth correction factor；L is the system loss factor.
In formula (1), antenna radiation pattern propagation factor FrAnd FtThe variable relevant with space coordinates, they be θ andLetter
Number.When radar is drawn, the Partial Variable need to be in real time only calculated, so as to reach the purpose of simplified operation.Therefore, abbreviation formula (1)
Afterwards on any directionDetection range can be expressed as：
In formula, RmaxIt is maximum detectable range of the radar on antenna gain direction in free space, it is clear that now Fr=Ft
=1.Think that antenna completes transmitting for monostatic radar, i.e. a frame antenna and receives task simultaneously in radar coverage-diagram modeling.
Thus it is possible to think Fr=Ft.Final formula is:
Therefore it is determined that after the maximum detectable range of radar, want accurately to determine shape of the radar in space exploration region,
The exact expression of the antenna radiation pattern function of radar must be just given.The present invention is directed to typical F (θ)=sin θ/θ type radars
Carry out visual presentation.
In sum, presently preferred embodiments of the present invention is these are only, is not intended to limit the scope of the present invention.
All any modification, equivalent substitution and improvements within the spirit and principles in the present invention, made etc., should be included in of the invention
Within protection domain.
Claims (8)
1. a kind of radar coverage display methods by the influence of topography towards digital earth, it is characterised in that including as follows
Step：
Step 1, radar coverage is sampled, specially：
In the range of radar directional pattern, by azimuth direction by setting step-lengthIt is divided into Pieces section；Then
On the center of radar directional pattern, setting step-length is pressed into each section of azimuth direction on pitching angular directionIt is divided into Layers section, section one radar lobe of correspondence of azimuth direction；
Step 2, for each radar lobe, on lobe border by setting step-length sampling, obtain sampled point；Judge radar lobe
Each sampled point and the position relationship of corresponding topographic(al) point, determine whether sampled point will be adapted to topographic(al) point according to position relationship on border
On, detailed process is：
S20, first from radar center point, along radar lobe the latter half, counterclockwise, find first by ground
The sampled point that shape is blocked；Calculate the height Ht of the corresponding topographic(al) point of radar center point R, and central point height value Hr；If
Hr>Ht, continues the height relationships of more next sampled point and corresponding topographic(al) point, until finding first height value less than correspondence
The sampled point of topographic(al) point height value, then the sampled point is the sampled point by the influence of topography, is defined as P points, correspondence topographic(al) point definition
It is Q points, this is adapted at topographic(al) point Q by the sampled point of the influence of topography；
S21, the sampled point for continuing to judge behind on lobe border according to the method for S20, it is next by the influence of topography until finding
Sampled point, be defined as current sampled point P1 to be modified；Connection RQ simultaneously extends, and extended line is determined with the intersection point on radar lobe border
Justice is Q '；Use the length ‖ QQ ' ‖ of line segment QQ ' divided by the number of sampled point between radar lobe border PQ ', obtain step delta r1；
From Q points, along QQ ' directions advance Δ r1, the point for being reached is defined as K1, and the corresponding topographic(al) points of K1 are Q1；Connection RQ1,
And extend, extended line is defined as sampled point Q1 ' with the intersection point on radar lobe border, then judges：
If 1., ‖ RQ1 ‖<‖ RQ1 ' ‖, then topographic(al) point Q1 P1 is adapted at Q1 still in radar maximum detection range；
If 2., ‖ RQ1 ‖ >=‖ RQ1 ' ‖, perform S22；
1. and 2. judged according to situation for all sampled points between PQ ', all sampled points are sentenced between complete PQ '
It is disconnected, finally perform S22；
S22, hypothesis P1 are previous processed sampled point, and are corrected at Q1；Use the length ‖ Q1Q1 ' ‖ of line segment Q1Q1 '
Divided by the number of sampled point between radar lobe border Q1Q1 ', step delta r is obtained；From Q1 points, before Q1Q1 ' directions
The point for entering Δ r arrival is defined as K points, judges：
If A, K point are higher than corresponding topographic(al) point：Next sampled point P2 of sampled point P1 is adjusted at K points, along Q1Q1 '
The direction Δ r that readvances reaches new K points, subsequent sampling point is adjusted at new K points, then judged：
If new K points reach radar lobe boundary point Q1 ', analysis terminates, and performs step 3；
If new K points do not reach radar lobe boundary point Q1 ', execution S22 is returned to；
If B, K point are lower than corresponding topographic(al) point, S23 is performed；
S23, hypothesis P1 are previous processed sampled point, and are corrected at Q1；From Q1 points, along Q1Q1 ' directions
Advance Δ r reaches K points；Judge：If K points are lower than corresponding topographic(al) point Q2, but the extended line of RQ2 and meeting at a little for radar border
Q2 ', and ‖ RQ2 ‖>‖ RQ2 ' ‖, then be adapted to Q2 ' places by P2 points, and analysis terminates, and performs step 3；If ‖ RQ2 ‖≤‖
RQ2 ' ‖, P2 is adapted at Q2, and the Δ r that moves on reaches another new point, is returned and is performed S23；
Step 3, according to the revised sampled point of step 2, new radar lobe border is drawn, so as to obtain the thunder under the influence of topography
Up to investigative range.
2. a kind of radar coverage display methods by the influence of topography towards digital earth as claimed in claim 1, its
It is characterised by, each radar lobe is processed using the method for step 2 simultaneously, after obtaining result, then using the method for step 3
Draw radar bounds.
3. a kind of radar coverage display methods by the influence of topography towards digital earth as claimed in claim 1, its
It is characterised by, the borderline sampled point of radar lobe uses the expression way of spherical coordinates, asθ is the angle of pitch,It is azimuthRepresent at pitching angle theta and azimuthThe corresponding detection range of upper radar.
4. a kind of radar coverage display methods by the influence of topography towards digital earth as claimed in claim 3, its
It is characterised by, detection range
In formula, PtIt is the power of transmission signal；τ is pulse bandwidth；GtIt is transmitting antenna power gain；GrIt is reception antenna power
Gain；σ is cross section of radar targets；λ is wavelength；FtIt is the directional diagram propagation factor from transmitting antenna to target；FrIt is from target
To the directional diagram propagation factor of reception antenna；K is Boltzmann constant；TsIt is reception system noise temperature；D0It is detecting factor；
CBIt is bandwidth correction factor；L is the system loss factor.
5. if claim 4 is by a kind of described radar coverage display methods by the influence of topography towards digital earth,
Characterized in that, detection rangeIt is reduced to following form：
Wherein, RmaxIt is maximum detectable range of the radar on antenna gain direction in free space.
6. if claim 3 is by a kind of described radar coverage display methods by the influence of topography towards digital earth,
Characterized in that, in the step 2, by the sample point coordinate in radar lobe from asking coordinate system to be converted to earth coordinates, then
The coordinate under earth coordinates with topographic(al) point is compared judgement, and radar border model is finally drawn under earth coordinates
Enclose.
7. if claim 3 is by a kind of described radar coverage display methods by the influence of topography towards digital earth,
Characterized in that, the step-lengthSpan be 0.1 degree to 10 degree.
8. if claim 3 is by a kind of described radar coverage display methods by the influence of topography towards digital earth,
Characterized in that, the span of the step delta θ is 0.1 degree to 10 degree.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710153349.2A CN106932772A (en) | 2017-03-15 | 2017-03-15 | A kind of radar coverage display methods by the influence of topography towards digital earth |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710153349.2A CN106932772A (en) | 2017-03-15 | 2017-03-15 | A kind of radar coverage display methods by the influence of topography towards digital earth |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106932772A true CN106932772A (en) | 2017-07-07 |
Family
ID=59433254
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710153349.2A Pending CN106932772A (en) | 2017-03-15 | 2017-03-15 | A kind of radar coverage display methods by the influence of topography towards digital earth |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106932772A (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109283499A (en) * | 2018-09-10 | 2019-01-29 | 西北工业大学 | Investigative range three-dimensional visualization method under a kind of active interference based on radar equation |
| CN110661582A (en) * | 2019-09-29 | 2020-01-07 | 北京电子工程总体研究所 | Method for performing radio wave communication in mountains and mountains by computer simulation |
| CN110708124A (en) * | 2019-09-29 | 2020-01-17 | 北京电子工程总体研究所 | Method for performing radio wave communication in bimodal terrain through computer simulation |
| CN111398918A (en) * | 2019-12-20 | 2020-07-10 | 中国科学院电子学研究所苏州研究院 | Radar detection capability analysis method under complex mountain environment |
| CN111859255A (en) * | 2020-07-17 | 2020-10-30 | 中国人民解放军78125部队 | Radar detection range calculation method under influence of terrain shielding |
| CN112859108A (en) * | 2021-01-28 | 2021-05-28 | 中国科学院南京土壤研究所 | Method for extracting under-forest vegetation coverage under complex terrain condition by using ground laser radar data |
| CN113093112A (en) * | 2021-03-23 | 2021-07-09 | 中国电子科技集团公司第二十九研究所 | Reconnaissance envelope calculation method and device for ground-oriented equipment under obstacle shielding and storage medium |
| CN117368869A (en) * | 2023-12-06 | 2024-01-09 | 航天宏图信息技术股份有限公司 | Visualization method, device, equipment and medium for radar three-dimensional power range |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102012501A (en) * | 2010-08-02 | 2011-04-13 | 电子科技大学 | Three-dimensional visualization correcting method for radar detecting range influenced by terrain |
| CN102360070A (en) * | 2011-06-10 | 2012-02-22 | 中国科学技术大学 | Receiving apparatus for ultra wideband impulse signal and ultra wideband impulse radar system |
| CN103913776A (en) * | 2014-04-23 | 2014-07-09 | 深圳市市政工程总公司 | Method for detecting spherical boulders in shield construction |
| CN104267396A (en) * | 2014-10-17 | 2015-01-07 | 中国电子科技集团公司第二十九研究所 | Imaging method for three-dimensional diagram of radar detection range |
-
2017
- 2017-03-15 CN CN201710153349.2A patent/CN106932772A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102012501A (en) * | 2010-08-02 | 2011-04-13 | 电子科技大学 | Three-dimensional visualization correcting method for radar detecting range influenced by terrain |
| CN102360070A (en) * | 2011-06-10 | 2012-02-22 | 中国科学技术大学 | Receiving apparatus for ultra wideband impulse signal and ultra wideband impulse radar system |
| CN103913776A (en) * | 2014-04-23 | 2014-07-09 | 深圳市市政工程总公司 | Method for detecting spherical boulders in shield construction |
| CN104267396A (en) * | 2014-10-17 | 2015-01-07 | 中国电子科技集团公司第二十九研究所 | Imaging method for three-dimensional diagram of radar detection range |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109283499A (en) * | 2018-09-10 | 2019-01-29 | 西北工业大学 | Investigative range three-dimensional visualization method under a kind of active interference based on radar equation |
| CN109283499B (en) * | 2018-09-10 | 2022-09-13 | 西北工业大学 | Radar equation-based three-dimensional visualization method for detection range under active interference |
| CN110661582A (en) * | 2019-09-29 | 2020-01-07 | 北京电子工程总体研究所 | Method for performing radio wave communication in mountains and mountains by computer simulation |
| CN110708124A (en) * | 2019-09-29 | 2020-01-17 | 北京电子工程总体研究所 | Method for performing radio wave communication in bimodal terrain through computer simulation |
| CN111398918A (en) * | 2019-12-20 | 2020-07-10 | 中国科学院电子学研究所苏州研究院 | Radar detection capability analysis method under complex mountain environment |
| CN111398918B (en) * | 2019-12-20 | 2020-12-04 | 中国科学院电子学研究所苏州研究院 | Radar detection capability analysis method under complex mountain environment |
| CN111859255A (en) * | 2020-07-17 | 2020-10-30 | 中国人民解放军78125部队 | Radar detection range calculation method under influence of terrain shielding |
| CN112859108A (en) * | 2021-01-28 | 2021-05-28 | 中国科学院南京土壤研究所 | Method for extracting under-forest vegetation coverage under complex terrain condition by using ground laser radar data |
| CN112859108B (en) * | 2021-01-28 | 2024-03-22 | 中国科学院南京土壤研究所 | Method for extracting vegetation coverage under forests under complex terrain condition by using ground laser radar data |
| CN113093112A (en) * | 2021-03-23 | 2021-07-09 | 中国电子科技集团公司第二十九研究所 | Reconnaissance envelope calculation method and device for ground-oriented equipment under obstacle shielding and storage medium |
| CN117368869A (en) * | 2023-12-06 | 2024-01-09 | 航天宏图信息技术股份有限公司 | Visualization method, device, equipment and medium for radar three-dimensional power range |
| CN117368869B (en) * | 2023-12-06 | 2024-03-19 | 航天宏图信息技术股份有限公司 | Visualization method, device, equipment and medium for radar three-dimensional power range |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106932772A (en) | A kind of radar coverage display methods by the influence of topography towards digital earth | |
| US9430871B2 (en) | Method of generating three-dimensional (3D) models using ground based oblique imagery | |
| CN106249237B (en) | Big Squint SAR frequency domain imaging method under a kind of curvilinear path | |
| JP6328789B2 (en) | Method and apparatus for determining angle of arrival (AOA) in a radar warning receiver | |
| US9330504B2 (en) | 3D building model construction tools | |
| CN102226840B (en) | Radar cross-section layered calculation method of ship target within atmospheric duct range | |
| CN112070894B (en) | Real environment navigation multipath real-time simulation method, device, medium and electronic equipment | |
| CN109900298A (en) | A kind of vehicle location calibration method and system | |
| CN105137430B (en) | The sparse acquisition of echo of forward sight array SAR a kind of and its three-D imaging method | |
| CN110146089A (en) | A kind of crusing robot localization method based on Embedded A I computing platform | |
| CN107918115B (en) | Radar target positioning method based on multipath utilization | |
| CN111398918B (en) | Radar detection capability analysis method under complex mountain environment | |
| CN110412571A (en) | Synthetic aperture radar three-dimensional imaging method based on electromagnetism vortex wave | |
| CN106291496B (en) | Moving target and ground combination scattering modeling method based on four path models | |
| CN109283499A (en) | Investigative range three-dimensional visualization method under a kind of active interference based on radar equation | |
| CN107678022A (en) | The radar Direct Three-dimensional space target positioning method utilized based on multipath | |
| CN105160698B (en) | A kind of trigonometric ratio ray trace method for searching path | |
| CN106646421B (en) | MIMO radar waveform co-design method based on three-dimensional nonuniform noise | |
| CN110456304A (en) | Airborne DF and location method | |
| CN112379603B (en) | Compensating system and method for installation eccentricity of strapdown seeker in radio frequency guidance simulation | |
| CN109738890A (en) | A method of distance figure is generated based on missile-borne Bistatic SAR range Doppler image | |
| JP7220246B2 (en) | Position detection method, device, equipment and readable storage medium | |
| CN115575907A (en) | Method for tracing surface crawling wave of complex target model based on planar grid | |
| CN115480493B (en) | Maneuvering control method for concealing approaching target in air combat simulation environment | |
| Cacciamano et al. | A coordinate registration technique for OTH sky-wave radars based on 3D ray-tracing and sea-land transitions |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170707 |