CN109633583A - The target shape inversion method of Surface scan Polarization scattering data is surrounded based near field - Google Patents
The target shape inversion method of Surface scan Polarization scattering data is surrounded based near field Download PDFInfo
- Publication number
- CN109633583A CN109633583A CN201811230301.8A CN201811230301A CN109633583A CN 109633583 A CN109633583 A CN 109633583A CN 201811230301 A CN201811230301 A CN 201811230301A CN 109633583 A CN109633583 A CN 109633583A
- Authority
- CN
- China
- Prior art keywords
- target
- scattering
- polarization
- parameter
- near field
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The present invention relates to a kind of target shape inversion methods that Surface scan Polarization scattering data are surrounded based near field comprising the steps of: S1, according to target length direction, building Target near field scan elliptic cylindrical shape enveloping surface；S2, transmitting antenna and receiving antenna using two cross-dipoles as near-field scan obtain the distributed data of the near-field scattering function on Target ellipse column enveloping surface；S3, calculate Target near field scattering function distributed data Huynen parameter；S4, it is based on Huynen parameter, inverting is carried out to target shape contour feature.The present invention can be realized the diagnosis and positioning of target strong scattering source or scattering center under different polarized states, gross data reference is provided for the low detectable geometric shape design of radar target, and changes using polarization parameter caused by target geometric shape and provides important reference frame for radar target acquisition and identification.
Description
Technical field
The present invention relates to the methods of radar target Polarization scattering device modeling and characterization, in particular to a kind of near field that is based on to wrap
Envelope surface scans the target shape inversion method of Polarization scattering data, belongs to radar target signature analysis and the technology of target identification is led
Domain.
Background technique
Radar target polarization scattering characteristics not only reflect target radar scattering cross-section under different polarized states (RCS,
Radar Cross Section) dimensional information variation, the amplitude and phase difference between different POLARIZATION CHANNELs can also be passed through
The geometric shape information of inverting target.
Target Polarization scattering data can pass through Inverse Synthetic Aperture Radar (ISAR, the Inverse of turntable target
Synthetic Aperture Radar) imaging is to obtain, and range resolution is realized by radar incidence wave bandwidth, orientation
It is realized to resolution ratio by turntable incidence angle varying width, thus obtains target in distance to the two-dimensional imaging with orientation,
Target geometry feature and strong scattering source or scattering center distribution characteristics are obtained, it is " Huang Peikang, blackish red that for details, reference can be made to documents
At, Xu little Jian, radar target signature [M], Beijing, Electronic Industry Press, 2005 ".But any target is all three-dimensional space
In solid geometry, on the one hand cause the target geometry of shadow region special blocking for incident electromagnetic wave due to space multistory geometry
Sign can not be imaged, and there are multiple target geometrical characteristics in the same imaging point in ISAR two-dimensional imaging for another aspect objective
Folded in resolution unit covers effect, these factors are unfavorable for determining the inverting of target shape geometrical characteristic and target strong scattering center
Position and diagnosis.
In order to carry out inverting to target geometric profile feature, in document, " Xing Shuguang, Lv Xiaode, Ding Chibiao are swept based on cylinder
Retouch the rcs measurement technique study [J] of Near-Field Radar Imaging, radar journal, 2015,4 (2): 172-177 " in, it is close to propose a kind of target
Field cylinder scans three-D imaging method, and incidence wave bandwidth realizes range resolution, the vertical and horizontal observation point along cylindrical surface
Not Shi Xian pitching to the resolution ratio with orientation.In order to meet complicated TV university scaled target imaging demand, the incident wave frequency of one side
Rate, angle and cylinder longitudinal sampling interval are very small, guarantee that target area, another aspect incidence wave is completely covered in imaging region
Bandwidth wants bigger, and long enough is wanted to guarantee higher imaging resolution in scanning cylindrical surface, but whether being to pass through electromagnetic scattering
Emulation, or by darkroom measurement experiment, all there is very big workload in the three-dimensional imaging of target cylinder scanning.
Several Kinds of Target Polar theory shows that the variation of target geometry and posture can be decomposed by target scattering matrix
Polarization parameter characterize.In document " Huynen, J.R.Phenomenological theory of radar targets
[D] .Delft, The Netherlands:University of Technology 1970 " in, Huynen divides collision matrix
Solution is the polarization parameters with physical meaning such as characterization target scattering amplitude, orientation, symmetry.Document " Cameron, W.L.,
Youssef,N.N.,Leung,L.K.Simulated polarimetric signatures of primitive
geometrical shapes[J].Geoscience and Remote Sensing,IEEE Transactions on,
1996,34 (3): 793-803 " in, target scattering is classified as conveyor screw, dihedral angle, narrow dihedral angle, trihedral angle, circle by Cameron
A variety of scatterer types such as column, dipole, a quarter wave device, never ipsilateral reflects the knot of target geometry for they
Structure variation.In document " Paladini, R., Ferro Famil, L., Pottier, E., etc.Point Target
Classification via Fast Lossless and Sufficient Omega-Psi-Phi Invariant
Decomposition of High-Resolution and Fully Polarimetric SAR/ISAR Data[J]
.Proceedings of the IEEE, 2013,101 (3): 798-830 " in, Palandini in actual measurement then to being unsatisfactory for
The asymmetric collision matrix of reciprocity carries out that energy is lossless, polarization invariance parameter decomposition, and obtaining one group, to describe target several
The polarization parameter of what feature.In addition, there are bigger differences for the polarization parameter of different target, same target is in different postures
Also there can be bigger polarization parameter variation under angle, therefore target polarization parameter can be used as one group of target scattering feature and be used for
Radar target recognition research, for details, reference can be made to document " Baird, C., Kersey, W., Giles, R.,
etc.Classification of targets using optimized ISAR Euler imagery[C].In
Defense and Security Symposium,2006；62100A-62100A-11".
Based on above-mentioned, the present invention proposes a kind of target shape inverting side that Surface scan Polarization scattering data are surrounded based near field
Method has solved the problems, such as of the existing technology and defect.
Summary of the invention
The object of the present invention is to provide a kind of target shape inverting sides that Surface scan Polarization scattering data are surrounded based near field
Method realizes the diagnosis and positioning of target strong scattering source or scattering center under different polarized states, is that radar target is low detectable
Geometric shape design provides gross data reference, and is radar target using the variation of polarization parameter caused by target geometric shape
Detection provides important reference frame with identification.
To achieve the above object, the present invention provides a kind of target shape that Surface scan Polarization scattering data are surrounded based near field
Inversion method comprising the steps of:
S1, according to target length direction, building Target near field scan elliptic cylindrical shape enveloping surface；
S2, transmitting antenna and receiving antenna using two cross-dipoles as near-field scan obtain Target ellipse column
The distributed data of near-field scattering function on shape enveloping surface；
S3, calculate Target near field scattering function distributed data Huynen parameter；
S4, it is based on Huynen parameter, inverting is carried out to target shape contour feature.
In the S1, using the reference axis of the length direction of target as center axis, formed by target be fully enclosed in
Elliptic cylindrical shape enveloping surface, and the cylinder length of the elliptic cylindrical shape enveloping surface is greater than target length, the major semiaxis of oval cross section and short
Maximum sized half of the length of semiaxis respectively than two orthogonal directions in object cross section is 1m long.
In the S2, comprising the following steps:
S21, on elliptic cylindrical shape enveloping surface, be parallel to the tangential of the direction x along enveloping surface and place two orthogonal electricity respectively
Type or magnetic type dipole antenna, respectively as transmitting dipole antenna and receiving dipole antenna；
S22, the back scattering electric field or magnetic field data that target is obtained by measurement, and then calculate corresponding scattering and write a letter in reply
Number, is denoted as:
Wherein, p, q=V, H respectively indicate the polarization of transmitting dipole antenna and receiving dipole antenna；V indicates vertical
Polarization, direction and the oval cross section of elliptic cylindrical shape enveloping surface are tangent；H indicates horizontal polarization, is oriented parallel to elliptic cylindrical shape encirclement
The cylinder in face is axial；
S23, it is scattered according to near-field targetAnd transmitting dipole antenna and receiving dipole is respectively set
The different polarized states of antenna obtain the distributed data of the near field Complete polarimetry function on Target ellipse column enveloping surface, i.e.,
Target near field collision matrix are as follows:
Wherein,Indicate the azimuthal variation of the oval cross section along elliptic cylindrical shape enveloping surface；X indicates elliptic cylindrical shape enveloping surface
On be parallel to the distance change of cylinder central axis.
In the S22, the cylinder central axis along elliptic cylindrical shape enveloping surface carries out equidistant intervals measurement, obtains target
Back scattering electric field or the longitudinal data in magnetic field；Oval cross section along elliptic cylindrical shape enveloping surface carries out isoazimuth interval measurement,
Obtain the back scattering electric field of target or the lateral data in magnetic field.
In the S3, Huynen parameter includes:
M indicates maximum same polarization response amplitude；
ψ,τmThe rotation angle of elliptic polarization base and oval angle when respectively indicating maximum same polarization response, and -90 °≤ψ≤
90 °, -45 °≤τm≤45°；
ν indicates the phase difference between Con-eigenvalue, and -45 °≤ν≤45 °；
γ is target signature angle or polarizing angle, indicate target to incident polarized sensibility, and 0 °≤γ≤45 °.
It is K matrix by Target near field collision matrix S-transformation in the S3, to realize the rapid solving of Huynen parameter,
Comprising the following steps:
S31, to collision matrix S-transformation, obtain the expression-form of K matrix are as follows:
Wherein,Indicate Kronecker product；And:
S32, each component for setting K matrix, and indicate are as follows:
Wherein, A0,B0It is angle of orientation invariance parameter；
K matrix is expressed as Huynen parameter ψ, τm, Euler's rotation of ν, specifically:
S33, by the substitution of variable to the formula in S32, release and Huynen parameter quickly calculated by K matrix component
It is expressed as follows:
Wherein, the expression way of parameter L and M are as follows:
In the S32, SD, U (ψ), U (τm), U (ν) is respectively indicated are as follows:
In the S4, by the Huynen parameter of obtained collision matrix, obtains target and polarize constant polarization characteristic
Parameter establishes its being associated between target geometric shape by the physical meaning and variation range of polarization characteristic parameter, thus
Pass through the shape geometrical characteristic of the variation inverting target of polarization characteristic parameter.
In the S4, comprising the following steps:
S41, the polarization characteristic parameter for polarizing constant according to Huynen parameter, acquisition target are as follows: indicate target maximum homopolarity
Change the m of response amplitude, and indicates the γ of target signature angle or polarizing angle；
S42, target scattering response lag is arranged to parameter m, more than or equal to the region quilt of the target scattering response lag
Regard as target scattering region；
S43, the change threshold that parameter γ is arranged position object edge feature by the γ value lower than the change threshold.
In conclusion the target shape inverting side provided by the present invention for surrounding Surface scan Polarization scattering data based near field
Method, can by target polarization parameter inverting target shape geometrical characteristic, realize under different polarized states target strong scattering source or
The diagnosis and positioning of scattering center, the only low detectable geometric shape design of radar target does not provide gross data reference, instead
Come over to change using polarization parameter caused by target geometric shape and provides important reference frame for radar target acquisition and identification.
Detailed description of the invention
Fig. 1 is the schematic diagram of the Target near field scanning elliptic cylindrical shape enveloping surface in the present invention；
Fig. 2 is the calculation flow chart of the Huynen parameter of the distributed data of the Target near field scattering function in the present invention；
Fig. 3 is that the near field of the metal plate in the present invention surrounds the schematic diagram of Surface scan Polarization scattering data result；
Fig. 4 is the schematic diagram of the near field polarization scan data Huynen parameter decomposition result of the metal plate in the present invention；
Fig. 5 is the schematic diagram of the metal plate Edge Gradient Feature based on Huynen parameter in the present invention；
Fig. 6 is the process of the target shape inversion method that Surface scan Polarization scattering data are surrounded based near field in the present invention
Figure.
Specific embodiment
Below in conjunction with FIG. 1 to FIG. 6, by preferred embodiment to technology contents of the invention, construction feature, reached purpose
And effect is described in detail.
As shown in fig. 6, for the target shape inverting provided by the invention for surrounding Surface scan Polarization scattering data based near field
Method comprising the steps of:
S1, according to target length direction, building Target near field scan elliptic cylindrical shape enveloping surface；
S2, transmitting antenna and receiving antenna using two cross-dipoles as near-field scan obtain Target ellipse column
The distributed data of near-field scattering function on shape enveloping surface；
S3, calculate Target near field scattering function distributed data Huynen parameter；
S4, it is based on Huynen parameter, inverting is carried out to target shape contour feature.
Under Near Field, target scattering is related to radar antenna, target observation distance, i.e. near-field target scattering needs comprehensive
Close the influence for considering radar directional pattern and observed range.For this purpose, proposing to need near field observed range and radar antenna class
Type progress is unified and standard, to form the canonical representation form of near-field target scattering.
The electromagnetic radiation characteristic for being uniformly based primarily upon antenna of radar antenna type.Variation charge and variation electricity on antenna
Stream is as the radiation source of excitation electromagnetic wave, themselves constitutes dipole radiation.Actual antenna is then considered as serving as reasons
Many dipoles are composed, and the electromagnetic field that antenna is excited is considered as the folded of the electromagnetic field excited for these dipoles
Add.According to electromagnetic field principle of stacking, any incident or scattering field can be analyzed to infinitesimal desired electrical or magnetic dipole
The combination of radiation field.
Using ideal dipole as radar emission and receiving antenna, then near-field target scattering can indicate the present invention are as follows:
Wherein, u, v=e, m respectively indicate the electric or magnetic type of transmitting and receiving dipole, and wherein e indicates electric type, m
Indicate magnetic type；RiTo emit the distance between dipole and target surface；RsBetween receiving dipole and target surface away from
From；L is scattering complex function, and the scattering complex function is related with target scattering characteristics, dipole-type and orientation.
Assuming that transmitting and receiving antenna are the dipole of electric type, incident electromagnetic wave is ideal electric type dipole spoke
The non-homogeneous spherical wave penetrated, then scattering complex function l may be expressed as:
Wherein, λ is electromagnetic wavelength；η is the wave impedance of free space, and Idl is the intensity for emitting electric type dipole, Es
Expression propagates to the electromagnetic wave electric field intensity of receiving dipole antenna.In order to express easily, subscript is omitted in the l in above formula
" ee ", indicate transmitting and it is received be all electric type dipole antenna.
Further, near field observed range is subjected to unified standard.It is required that emitting dipole and receiving dipole to object table
The distance in face is used uniformly unit distance, i.e. Ri=Rs=1m, as a result, near-field target scattering then indicate are as follows:
In general, the geometry of target surface be not rule variation, observation antenna to target surface unit away from
It will be complex from also.Therefore, need to carry out certain approximate agreement in practical applications.And in the present invention, then arrange to use
The concept of Target ellipse column enveloping surface carrys out specification antenna range-to-go.That is in the S1, with the length direction of target
Reference axis be center axis, formed target be fully enclosed in interior elliptic cylindrical shape enveloping surface, and the elliptic cylindrical shape enveloping surface
Cylinder length is greater than target length, and the length of the major semiaxis of oval cross section and semi-minor axis is more orthogonal than in object cross section two respectively
The also long 1m of the maximum sized half in direction.
As shown in Figure 1, using the reference axis in length of aircraft direction as center axis, being formed aircraft so that target is aircraft as an example
It is fully enclosed in interior elliptic cylindrical shape enveloping surface, and requires the cylinder length of this elliptic cylindrical shape enveloping surface greater than the length of aircraft
Degree, while the length of the major semiaxis of oval cross section and semi-minor axis is respectively than the extreme length of two orthogonal directions in aircraft cross section
The also long 1m of half.
In the S2, comprising the following steps:
S21, on elliptic cylindrical shape enveloping surface, be parallel to the tangential of the direction x along enveloping surface and place two orthogonal electricity respectively
Type or magnetic type dipole antenna, respectively as transmitting dipole antenna and receiving dipole antenna；As shown in fig. 1, H, V
Respectively indicate the placement form of two dipoles；
S22, the back scattering electric field or magnetic field data that target is obtained by measurement, and then calculate corresponding scattering and write a letter in reply
Number, is denoted as:
Wherein, p, q=V, H respectively indicate the polarization of transmitting dipole antenna and receiving dipole antenna；V indicates vertical
Polarization, direction and the oval cross section of elliptic cylindrical shape enveloping surface are tangent；H indicates horizontal polarization, is oriented parallel to elliptic cylindrical shape encirclement
The cylinder in face is axial；
S23, it is scattered according to near-field targetAnd transmitting dipole antenna and receiving dipole day is respectively set
The different polarized states of line (namely by changing the polarization of dipole antenna, realize the not same polarization of transmitting with receiving antenna
State), dummy level polarization is scattered with the near-field target under vertical polarization, and the near field obtained on Target ellipse column enveloping surface is complete
The distributed data of Polarization scattering function, i.e. Target near field collision matrix are as follows:
Wherein,Indicate the azimuthal variation of the oval cross section along elliptic cylindrical shape enveloping surface；X indicates elliptic cylindrical shape enveloping surface
On be parallel to the distance change of cylinder central axis.
In the S22, the cylinder central axis along elliptic cylindrical shape enveloping surface carries out equidistant intervals measurement, obtains target
Back scattering electric field or the longitudinal data in magnetic field；Oval cross section along elliptic cylindrical shape enveloping surface carries out isoazimuth interval measurement,
Obtain the back scattering electric field of target or the lateral data in magnetic field.
For Target near field collision matrix S, it is assumed that target scattering meets reciprocity, then collision matrix is symmetrical matrix, i.e.,
SHV=SVH.So, Huynen, which is decomposed, decomposes the symmetrical collision matrix S for meeting scattering reciprocity are as follows:
S=U*(ψ,τm,ν)SDUH(ψ,τm,ν)；
Wherein, subscript * indicates conjugation, and subscript H indicates conjugate transposition；And:
That is, for U (ψ), U (τm), U (ν) can be respectively indicated are as follows:
Above-mentioned decomposed form is also referred to as congruent transformation, U (ψ, τm, ν) and indicate that congruent transformation matrix, m indicate maximum same polarization
Response amplitude；The absolute phase of ρ expression target scattering；ψ,τmThe rotation of elliptic polarization base when respectively indicating maximum same polarization response
Angle and oval angle, and -90 °≤ψ≤90 °, -45 °≤τm≤45°；Meanwhile τmValue indicate target perpendicular to radar line of sight
Symmetry in plane, τm=0 ° is symmetry scatterer, τm≠ 0 ° is asymmetric scatterer；ν is indicated between Con-eigenvalue
Phase difference, and -45 °≤ν≤45 °；ν is related to the scattering imaging of echo simultaneously, ν=0 ° when odd times scatter, ν when even scatters
=± 45 °, therefore also referred to as target skip angular；γ is target signature angle or polarizing angle, indicates target to incident polarized sensitivity
Property, and 0 °≤γ≤45 °.
Parameter ψ, τ mentioned abovem, m, ν, γ is Huynen parameter.Due to directly by the decomposition of collision matrix S come
It calculates each Huynen parameter to be not easy to realize, therefore in the present invention, considers collision matrix S-transformation to be Kennaugh square
Battle array (K matrix), to realize the rapid solving of Huynen parameter.
As shown in Fig. 2, in the S3, comprising the following steps:
S31, the incident wave vector and scattering Stokes vector of K matrix expression target similar with the definition of collision matrix
Transformation relation；And using the transformation relation between incident wave vector and Stokes vector, K can be obtained to collision matrix S-transformation
The form of matrix are as follows:
Wherein,Indicate Kronecker product；And:
S32, each component for setting K matrix, and indicate are as follows:
Wherein, A0,B0It is angle of orientation invariance parameter；Other each components are to be calculated according to the mode in S31；
By matrixing it can be proved that K matrix can be expressed as Huynen parameter ψ, τm, Euler's rotation of ν, expression are as follows:
S33, it can be released logical by the substitution of variable to the formula in S32 in addition to target scattering absolute phase parameter ρ
It crosses K matrix component and quickly calculates being expressed as follows for Huynen parameter:
Wherein, the expression way of parameter L and M are as follows:
In a preferred embodiment of the invention, it is surveyed using side length for the metal plate of 1m × 0.6m as target
The angle of examination, the metal plate and x-axis is 15 °.The length that elliptic cylindrical shape enveloping surface is arranged is 2m, is carried out along cylinder central axis
The equidistant intervals of 0.005m measure and sampling, and 0.5 ° of isoazimuth interval measurement and sampling are carried out along oval cross section, incident
Wave frequency rate is 15GHz, obtains the near-field scan of complete polarization elliptic cylindrical shape enveloping surface as a result, as shown in Figure 3.Wherein, image dynamic
Range unit is dB, and Fig. 3 (a) left-half is HH polarization, and right half part is HV polarization, and Fig. 3 (b) left-half is the pole VH
Change, right half part is VV polarization.
Huynen decomposition then carried out to the near field Polarization scattering data of metal plate, obtained Huynen parameter is respectively such as
Shown in Fig. 4.Wherein, Fig. 4 (a) display -35dB≤10log10M≤15dB, and width is responded by the maximum same polarization of Fig. 4 (a) target
Value can orient target strong scattering region, in target scattering region.Fig. 4 (b) shows γ≤1 0≤tan, wherein low γ value
Region shows object edge contoured profile.Fig. 4 (c) shows -90 °≤ψ≤90 °, and Fig. 4 (d) shows -10 °≤τm≤ 10 °, Fig. 4
(e) -45 °≤ν≤45 °, and ψ, τ are shownm, ν goes to zero, and shows metal plate under the conditions of current pose and antenna polarization
Target property.
In the S4, by the Huynen parameter of obtained collision matrix, obtains target and polarize constant polarization characteristic
Parameter establishes its being associated between target geometric shape by the physical meaning and variation range of polarization characteristic parameter, thus
Pass through the shape geometrical characteristic of the variation inverting target of polarization characteristic parameter.
In the S4, comprising the following steps:
S41, the polarization characteristic parameter for polarizing constant according to Huynen parameter, acquisition target are as follows: indicate target maximum homopolarity
Change the m of response amplitude, and indicates the γ of target signature angle or polarizing angle；
S42, target scattering response lag is arranged to parameter m, more than or equal to the region quilt of the target scattering response lag
Target scattering region is regarded as, the region less than the target scattering response lag is identified as ambient noise, to orientTarget area in plane；
S43, according to the correlation between target geometric shape and Huynen parameter, determine that object edge scattering shows as idol
Extremely sub- scattering properties, parameter γ is close to zero at this time；Thus the change threshold of parameter γ is set, by lower than the change threshold
γ value positions object edge feature.
By taking the Huynen parameter of the metal plate scattering in Fig. 4 as an example, target area is determined by m setup parameter threshold value first
It is the target area bianry image of m >=1 shown in such as Fig. 5 (a) of domain.Secondly the γ of target area is judged, obtains low γ ginseng
Number scattering point is tan as shown in Fig. 5 (b)2The plate edge feature (expression of grey dot) of γ≤1.By low in Fig. 5 (b)
γ parameter scattering point is compared with object edge as can be seen that object edge contour feature not only may be implemented in Huynen parameter m, γ
Inverting, while m also located the distribution in target strong scattering source.
In conclusion the target shape inverting side provided by the present invention for surrounding Surface scan Polarization scattering data based near field
Method is established under target-based coordinate system and target is enclosed in interior elliptic cylindrical shape enveloping surface, is carried out using dipole antenna to target
Complete polarization scanning near field changes dipole antenna orientation to obtain the target scattering data under different polarized states；Pass through encirclement
Face laterally and longitudinally scanning result can directly obtain target two-dimensional scattering image, and then obtain target strong scattering source or scattering
Central cross and genesis analysis position, and relevant point is carried out to the target polarization scattering matrix of target strong scattering center position
Solution, using Coherent Targets polarization decomposing method obtain target polarize invariance parameter, finally by target polarization parameter with it is several
Correlation (physical meaning and variation range) between what shape carrys out inverting target geometric shape.
The present invention can realize target under different polarized states by target polarization parameter inverting target shape geometrical characteristic
The diagnosis and positioning of strong scattering source or scattering center, the only low detectable geometric shape design of radar target does not provide theoretical value
According to reference, provided in turn for radar target acquisition with identification using the variation of polarization parameter caused by target geometric shape important
Reference frame.
It is discussed in detail although the contents of the present invention have passed through above preferred embodiment, but it should be appreciated that above-mentioned
Description is not considered as limitation of the present invention.After those skilled in the art have read above content, for of the invention
A variety of modifications and substitutions all will be apparent.Therefore, protection scope of the present invention should be limited to the appended claims.
Claims (9)
1. it is a kind of based near field surround Surface scan Polarization scattering data target shape inversion method, which is characterized in that comprising with
Lower step:
S1, according to target length direction, building Target near field scan elliptic cylindrical shape enveloping surface；
S2, transmitting antenna and receiving antenna using two cross-dipoles as near-field scan obtain Target ellipse column packet
The distributed data of near-field scattering function in envelope surface；
S3, calculate Target near field scattering function distributed data Huynen parameter；
S4, it is based on Huynen parameter, inverting is carried out to target shape contour feature.
2. the target shape inversion method of Surface scan Polarization scattering data is surrounded based near field as described in claim 1, it is special
Sign is, in the S1, using the reference axis of the length direction of target as center axis, formed by target be fully enclosed in it is ellipse
Cylindrical enveloping surface, and the cylinder length of the elliptic cylindrical shape enveloping surface is greater than target length, the major semiaxis of oval cross section and short by half
Maximum sized half of the length of axis respectively than two orthogonal directions in object cross section is 1m long.
3. the target shape inversion method of Surface scan Polarization scattering data is surrounded based near field as described in claim 1, it is special
Sign is, in the S2, comprising the following steps:
S21, on elliptic cylindrical shape enveloping surface, be parallel to the tangential of the direction x along enveloping surface and place two orthogonal electric types respectively
Or magnetic type dipole antenna, respectively as transmitting dipole antenna and receiving dipole antenna；
S22, the back scattering electric field or magnetic field data that target is obtained by measurement, and then corresponding scattering complex function is calculated, note
Are as follows:
Wherein, p, q=V, H respectively indicate the polarization of transmitting dipole antenna and receiving dipole antenna；
V indicates vertical polarization, and direction and the oval cross section of elliptic cylindrical shape enveloping surface are tangent；H indicates horizontal polarization, is oriented parallel to
The cylinder of elliptic cylindrical shape enveloping surface is axial；
S23, it is scattered according to near-field targetAnd transmitting dipole antenna and receiving dipole antenna is respectively set
Different polarized states, obtain the distributed data of the near field Complete polarimetry function on Target ellipse column enveloping surface, i.e. target is close
Field collision matrix are as follows:
Wherein,Indicate the azimuthal variation of the oval cross section along elliptic cylindrical shape enveloping surface；X indicates to put down on elliptic cylindrical shape enveloping surface
Row is in the distance change of cylinder central axis.
4. the target shape inversion method of Surface scan Polarization scattering data is surrounded based near field as claimed in claim 3, it is special
Sign is, in the S22, the cylinder central axis along elliptic cylindrical shape enveloping surface carries out equidistant intervals measurement, obtains target
Back scattering electric field or the longitudinal data in magnetic field；Oval cross section along elliptic cylindrical shape enveloping surface carries out isoazimuth interval measurement,
Obtain the back scattering electric field of target or the lateral data in magnetic field.
5. the target shape inversion method of Surface scan Polarization scattering data is surrounded based near field as claimed in claim 3, it is special
Sign is, in the S3, Huynen parameter includes:
M indicates maximum same polarization response amplitude；
ψ,τmThe rotation angle of elliptic polarization base and oval angle when respectively indicating maximum same polarization response, and -90 °≤ψ≤90 °, -45 °
≤τm≤45°；
ν indicates the phase difference between Con-eigenvalue, and -45 °≤ν≤45 °；
γ is target signature angle or polarizing angle, indicate target to incident polarized sensibility, and 0 °≤γ≤45 °.
6. the target shape inversion method of Surface scan Polarization scattering data is surrounded based near field as claimed in claim 5, it is special
Sign is, is K matrix by Target near field collision matrix S-transformation in the S3, to realize the rapid solving of Huynen parameter,
Comprising the following steps:
S31, to collision matrix S-transformation, obtain the expression-form of K matrix are as follows:
Wherein,Indicate Kronecker product；And:
S32, each component for setting K matrix, and indicate are as follows:
Wherein, A0,B0It is angle of orientation invariance parameter；
K matrix is expressed as Huynen parameter ψ, τm, Euler's rotation of ν, specifically:
S33, by the substitution of variable to the formula in S32, release the expression that Huynen parameter is quickly calculated by K matrix component
It is as follows:
Wherein, the expression way of parameter L and M are as follows:
7. the target shape inversion method of Surface scan Polarization scattering data is surrounded based near field as claimed in claim 6, it is special
Sign is, in the S32, SD, U (ψ), U (τm), U (ν) is respectively indicated are as follows:
8. the target shape inversion method of Surface scan Polarization scattering data is surrounded based near field as claimed in claim 6, it is special
Sign is, in the S4, by the Huynen parameter of obtained collision matrix, obtains the polarization characteristic ginseng that target polarizes constant
Number, by the physical meaning and variation range of polarization characteristic parameter, establishes its being associated between target geometric shape, thus logical
Cross the shape geometrical characteristic of the variation inverting target of polarization characteristic parameter.
9. the target shape inversion method of Surface scan Polarization scattering data is surrounded based near field as claimed in claim 8, it is special
Sign is, in the S4, comprising the following steps:
S41, the polarization characteristic parameter for polarizing constant according to Huynen parameter, acquisition target are as follows: indicate that the same polarization of target maximum is rung
The m of amplitude is answered, and indicates the γ of target signature angle or polarizing angle；
S42, target scattering response lag is arranged to parameter m, the region more than or equal to the target scattering response lag is identified
For target scattering region；
S43, the change threshold that parameter γ is arranged position object edge feature by the γ value lower than the change threshold.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811230301.8A CN109633583A (en) | 2018-10-22 | 2018-10-22 | The target shape inversion method of Surface scan Polarization scattering data is surrounded based near field |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811230301.8A CN109633583A (en) | 2018-10-22 | 2018-10-22 | The target shape inversion method of Surface scan Polarization scattering data is surrounded based near field |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109633583A true CN109633583A (en) | 2019-04-16 |
Family
ID=66066525
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811230301.8A Pending CN109633583A (en) | 2018-10-22 | 2018-10-22 | The target shape inversion method of Surface scan Polarization scattering data is surrounded based near field |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109633583A (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110826643A (en) * | 2019-11-20 | 2020-02-21 | 上海无线电设备研究所 | Offshore target identification method based on polarized Euler feature fusion deep learning |
| CN111505627A (en) * | 2020-04-28 | 2020-08-07 | 湖南正申科技有限公司 | Life detection and positioning method for full-polarization radar |
| CN111880168A (en) * | 2020-08-04 | 2020-11-03 | 上海无线电设备研究所 | Target positioning method based on passive digital array radar |
| CN112347913A (en) * | 2020-11-05 | 2021-02-09 | 中国科学院国家空间科学中心 | Disaster-suffered building damage level estimation method and system based on rapid Huynen-Euler decomposition |
| CN112597434A (en) * | 2020-11-05 | 2021-04-02 | 中国科学院国家空间科学中心 | Rapid implementation method and system for polarized KHT decomposition |
| CN113671269A (en) * | 2021-07-05 | 2021-11-19 | 南京航空航天大学 | Spherical near-far field method based on scattering matrix principle |
| CN114814772A (en) * | 2022-05-09 | 2022-07-29 | 北京航空航天大学 | Single-target polarization scattering matrix-based RCS (radar cross section) calculation method for electrically large cluster targets |
| CN115144856A (en) * | 2022-09-05 | 2022-10-04 | 中国人民解放军国防科技大学 | Method and device for inverting space target parabolic antenna parameters of polarized radar |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103761741A (en) * | 2014-01-23 | 2014-04-30 | 电子科技大学 | Variation polarized SAR image partition method based on polarization target decomposition feature |
| CN104636752A (en) * | 2015-02-03 | 2015-05-20 | 中国科学院空间科学与应用研究中心 | Method for classifying scattering pyramids facing polarimetric synthetic aperture radar images |
| CN104931950A (en) * | 2015-06-23 | 2015-09-23 | 中国科学院空间科学与应用研究中心 | Target decomposition method based on model for fully-polarized synthetic aperture radar |
| CN108061883A (en) * | 2017-11-28 | 2018-05-22 | 上海无线电设备研究所 | The method of the near-field scattering function conversion missile target encounter echo of locally scattered source inverting |
-
2018
- 2018-10-22 CN CN201811230301.8A patent/CN109633583A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103761741A (en) * | 2014-01-23 | 2014-04-30 | 电子科技大学 | Variation polarized SAR image partition method based on polarization target decomposition feature |
| CN104636752A (en) * | 2015-02-03 | 2015-05-20 | 中国科学院空间科学与应用研究中心 | Method for classifying scattering pyramids facing polarimetric synthetic aperture radar images |
| CN104931950A (en) * | 2015-06-23 | 2015-09-23 | 中国科学院空间科学与应用研究中心 | Target decomposition method based on model for fully-polarized synthetic aperture radar |
| CN108061883A (en) * | 2017-11-28 | 2018-05-22 | 上海无线电设备研究所 | The method of the near-field scattering function conversion missile target encounter echo of locally scattered source inverting |
Non-Patent Citations (3)
| Title |
|---|
| DALLMANN T: "Discrimination of scattering mechanisms via polarimetric rcs imaging", 《IEEE ANTENNAS AND PROPAGATION MAGAZINE》 * |
| YONGCHEN LI: "Imagmg and structural feature decomposition of complex target using multi-aspect polarimetric scattering", 《SCIENCE CHINA INFORMATION SCIENCE》 * |
| 梁子长 等: "基于单位距离包围面的近场RCS定义及应用", 《制导与引信》 * |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110826643A (en) * | 2019-11-20 | 2020-02-21 | 上海无线电设备研究所 | Offshore target identification method based on polarized Euler feature fusion deep learning |
| CN110826643B (en) * | 2019-11-20 | 2022-08-19 | 上海无线电设备研究所 | Offshore target identification method based on polarized Euler feature fusion deep learning |
| CN111505627A (en) * | 2020-04-28 | 2020-08-07 | 湖南正申科技有限公司 | Life detection and positioning method for full-polarization radar |
| CN111880168A (en) * | 2020-08-04 | 2020-11-03 | 上海无线电设备研究所 | Target positioning method based on passive digital array radar |
| CN112347913A (en) * | 2020-11-05 | 2021-02-09 | 中国科学院国家空间科学中心 | Disaster-suffered building damage level estimation method and system based on rapid Huynen-Euler decomposition |
| CN112597434A (en) * | 2020-11-05 | 2021-04-02 | 中国科学院国家空间科学中心 | Rapid implementation method and system for polarized KHT decomposition |
| CN112597434B (en) * | 2020-11-05 | 2024-03-22 | 中国科学院国家空间科学中心 | Method and system for rapidly realizing polarization KHT decomposition |
| CN112347913B (en) * | 2020-11-05 | 2024-03-29 | 中国科学院国家空间科学中心 | Disaster-stricken building damage level estimation method and system based on fast Huynen-Euler decomposition |
| CN113671269A (en) * | 2021-07-05 | 2021-11-19 | 南京航空航天大学 | Spherical near-far field method based on scattering matrix principle |
| CN114814772A (en) * | 2022-05-09 | 2022-07-29 | 北京航空航天大学 | Single-target polarization scattering matrix-based RCS (radar cross section) calculation method for electrically large cluster targets |
| CN114814772B (en) * | 2022-05-09 | 2024-04-09 | 北京航空航天大学 | Electric large cluster target RCS calculation method based on single target polarization scattering matrix |
| CN115144856A (en) * | 2022-09-05 | 2022-10-04 | 中国人民解放军国防科技大学 | Method and device for inverting space target parabolic antenna parameters of polarized radar |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109633583A (en) | The target shape inversion method of Surface scan Polarization scattering data is surrounded based near field | |
| Liu et al. | Hybrid polarimetric GPR calibration and elongated object orientation estimation | |
| Zhang et al. | MIMO ground penetrating radar imaging through multilayered subsurface using total variation minimization | |
| Xu et al. | Polarimetric–anisotropic decomposition and anisotropic entropies of high-resolution SAR images | |
| Dallmann et al. | Discrimination of scattering mechanisms via polarimetric rcs imaging [measurements corner] | |
| Demirci et al. | Interpretation and analysis of target scattering from fully-polarized ISAR images using Pauli decomposition scheme for target recognition | |
| Sun et al. | The orientation estimation of elongated underground objects via multipolarization aggregation and selection neural network | |
| Dvorsky et al. | Polarimetric synthetic aperture radar imaging with radially polarized antennas | |
| Dvorsky et al. | Synthetic aperture radar 3-D polarimetry | |
| Nassib et al. | A dyadic target model for multistatic SAR/ISAR imaging | |
| Frezza et al. | Short-pulse electromagnetic scattering by buried perfectly conducting cylinders | |
| Wang et al. | A novel approach to full-polarimetric short-range imaging with copolarized data | |
| Zhang et al. | Characteristics analysis and image processing for full-polarization synthetic aperture radar based on electromagnetic scattering from flat horizontal perfect electric conducting reflector | |
| Duquenoy et al. | Study of dispersive and anisotropic scatterers behavior in radar imaging using time-frequency analysis and polarimetric coherent decomposition | |
| Theofanopoulos et al. | A novel THz radar imaging system using the radon transform | |
| Tao et al. | Comparison between monostatic and bistatic SAR image based on spatial spectrum analysis | |
| Chauveau et al. | A new algorithm of 3D image reconstruction of radar targets from ramp responses in low frequency | |
| del Rio et al. | Statistics of the degree of polarization | |
| Popov et al. | Detection of metal objects against mountain ranges using polarization transfer functions | |
| Baird et al. | Exploitation of ISAR imagery in Euler parameter space | |
| O'Reilly et al. | Target decomposition and polarimetric radar applied to concealed threat detection | |
| Ali et al. | Compact Measurement Apparatus for 3D Spherical Read Range Characterization of Chipless RFID | |
| Zhang et al. | Similarity-enhanced target detection algorithm using polarimetric SAR images | |
| Salman et al. | Scattering matrix decomposition of quad-polarized short-range ultrawideband Radar data | |
| Phruksahiran et al. | Improving SAR data processing with polarimetric reference functions in the range Doppler algorithm |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190416 |