CN108562899A - High-resolution polarimetric SAR target image rapid simulation method - Google Patents
High-resolution polarimetric SAR target image rapid simulation method Download PDFInfo
- Publication number
- CN108562899A CN108562899A CN201810362394.3A CN201810362394A CN108562899A CN 108562899 A CN108562899 A CN 108562899A CN 201810362394 A CN201810362394 A CN 201810362394A CN 108562899 A CN108562899 A CN 108562899A
- Authority
- CN
- China
- Prior art keywords
- light
- face element
- matrix
- radar
- scattering
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
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
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9076—Polarimetric features in SAR
-
- 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
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9004—SAR image acquisition techniques
- G01S13/9005—SAR image acquisition techniques with optical processing of the SAR signals
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)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a kind of high-resolution polarimetric SAR target image rapid simulation method, mainly solve the problem of that current polarimetric SAR image obtains that difficulty is big, of high cost, speed is slow and traditional ray tracing image simulation does not consider polarization characteristic.Implementation step is：Target scene model is constructed, radar parameter is set；The position for being irradiated to face element is determined according to the distance between each face element of target on opticpath and radar, calculates the scattering,single matrix and rescattering matrix of the face element；The scattering,single energy of face element is calculated according to incident ray and scattering,single matrix；The rescattering energy of face element is calculated according to incident ray and rescattering matrix；Gross energy is scattered according to scattering,single energy and rescattering energy balane；It completes to be imaged target scene according to scattering gross energy.The more true polarimetric SAR image of the easy quick obtaining of energy of the invention, reduces cost, can be used for target identification.
Description
Technical field
The invention belongs to Radar Technology field, more particularly to a kind of SAR target images emulation mode can be used for radar target
Discriminance analysis builds library.
Background technology
Synthetic aperture radar SAR is a kind of radar obtaining high-definition picture using microwave, it has round-the-clock, complete
The advantages of weather.It dexterously utilizes pulse compression technique, the method for synthetic aperture technique and some signal processings, with true
Small-bore antenna obtains orientation and apart from upward double high-definition pictures.In orientation, SAR utilizes small-bore day
Line to make beam angle become smaller, and then makes azimuth resolution become by the continuous movement of radar to synthesize big aperture
It is high；Upward, the linear FM signal of the big bandwidth of radar emission in distance carries out pulse after receiving target echo to echo
Compression so that distance to resolution ratio get higher.
With the continuous development of SAR technologies, the analogue technique of SAR also achieves good progress.SAR analogue techniques are profits
The course of work of SAR is simulated with computer, the final technology for realizing imaging.SAR imagings are main and the reflection characteristic of target with
And SAR system itself is related, and the reception due to transmitting, signal from signal and the processing to received signal are to the end
Imaging, many steps therein be required for aboard practical measure realize, by carrier of radar movement, radar system hardware
It is limited with the influence of signal processing technology, ultimately generating image, to generate error larger and uncontrollable, sees that [Li Jin are based on FPGA's
Carried SAR simulated radar echo studies the Chengdu [D]：University of Electronic Science and Technology .1-65].In addition, needed for the acquisition of SAR experimental datas
Cost it is relatively high, therefore, using computer combination target scattering characteristics, SAR image-forming principles and electromagnetism computer sim- ulation technology come
SAR image is generated, this design for SAR system, the researchs such as imaging verification and target characteristics analysis suffer from extremely important
Effect, analog simulation has become a kind of economical and effective of SAR imaging research and highly important method.This method can be with
Greatly reduce the cost of actual flying test, and image introducing physical condition is comprehensively controllable, can save search time and warp
Ji expense.
Traditional ray-tracing simulations SAR imagings, are the irradiations that the transmitting of electromagnetic wave is regarded as to a rule ray, pass through
The scattering and reflection of simulated light are imaged according to the size of its propagation path and scattering energy.But in imaging process, by
Polarization characteristic is propagated in the vector for the electromagnetic wave for not accounting for actual transmission, thus causes its image generated that cannot reflect mesh
Target Polarization scattering feature so that emulation gained image and the SAR image of practical not same polarization admission have larger difference, can not
Target Signature Analysis for polarization SAR and identification.
Invention content
It is an object of the invention in view of the above shortcomings of the prior art, propose a kind of high-resolution polarimetric SAR target image
Rapid simulation method keeps emulating image truer to reduce the difference of emulating image and practical polarimetric SAR image, and can be with
Target identification for polarization SAR.
The technical scheme is that：Small Triangular object model one by one forms target scene, with the light with polarization characteristic
Line irradiates target scene；The position that the small Triangular object model of target being irradiated to is determined using biggest advantage of light track algorithm, according to light
The posture of direction of illumination and small Triangular object model calculates the back scattering matrix at the Triangular object model, and then by each illuminated
The back scattering energy of face element is imaged, and implementation step includes as follows：
(1) parameter setting
Airborne radar height h, carrier aircraft flying speed v, light minimum irradiation angle, θ are setmin, azimuth resolution Δ A,
Range resolution Δ R, radar transmitting wave frequency ω, target surface permittivity ε, horizontal polarization directions electromagnetism wave amplitude
Spend EihWith vertical polarization directions electromagnetism wave amplitude Eiv；
(2) the target scene 3D models combined by Triangular object model are built, and are conducted into commercial matlab softwares, are extracted
The coordinates matrix T on the 3D model Triangular object models vertex, the orientation of radar is obtained most by apex coordinate matrix T and radar altitude h
Big value Amax, minimum value Amin, oblique distance maximum value RmaxWith oblique distance minimum value Rmin；
(3) parameter obtained according to (2), calculate radar bearing to sampling number K and distance to sampling number N；
(4) control radar does the linear uniform motion that speed is v, constantly emits beam during the motion, light passes through
In radar illumination to target scene, the light matrix of a K × N is formed；
(5) according to the parameter being arranged in (1), radar motion is in the process in each light direction of illumination r of position P { l, m }
{ l, m }, wherein l are the line number of light matrix, and 1≤l≤K, m are the columns of light matrix, 1≤m≤N；
(6) Triangular object model being irradiated to is determined：
(6a) is the light E of m according to the parameter in (1) and (5), Accounting Line Number l, columnsiTarget face on propagation path
The number ζ of member, and calculate the distance between radar and each target face element on propagation path Di；
(6b) from small to large sorts a distance, and the face element corresponding to minimum range is determined as being irradiated to by light
Face element, if the face element coordinate is in the q rows of apex coordinate matrix T, 1≤q≤K；
(7) following parameter is calculated according to Triangular object model：
(7a) calculates face element normal vector g according to line number q where Triangular object model；
(7b) calculates light and Triangular object model normal vector according to light direction of illumination r { l, m } and face element normal vector g
Between angle α；
(7c) according to the angle α in target surface permittivity ε in (1) and (7b) between light and face element normal vector,
The polarized Fresnel reflection coefficient R of calculated level1With the Fresnel reflection coefficient R of vertical polarization2；
Horizontal Fresnel reflection coefficient in angle α, (7c) in (7d) basis (7b) between light and face element normal vector
R1With vertical Fresnel reflection coefficient R2, calculate the back scattering matrix S of the scattering,single of face element1Backward with rescattering dissipates
Penetrate matrix S2；
(7e) is according to light E in (6)iWith the back scattering matrix S of middle scattering,single1, calculate the scattering,single energy of face element
E1s(l,m)；
(7f) is according to incident ray EiWith rescattering matrix S2, calculate the rescattering ENERGY E of face element2s(l,m)；
(8) step (6)-(7) are pressed, the scattering,single ENERGY E of all K × N light is calculated1sWith rescattering ENERGY E2s；
(9) by scattering,single ENERGY E1sWith rescattering ENERGY E2sThe sum of as pixel matrix Es, utilize business software
Matlab draws high-resolution polarimetric SAR image.
The invention has the advantages that
1) present invention is due in the polarization characteristic for using ray tracing add light when SAR image emulation, making to imitate
True process more has authenticity, and it is also more true to be formed by image.
2) polarization characteristic of the present invention since light is utilized, gained image have reacted the polarization characteristic of target, so as to
Data are provided and are supported for polarization SAR target identification and polarization SRA characteristics of image study.
Description of the drawings
Fig. 1 is the implementation flow chart of the present invention；
Fig. 2 is simulation result diagram of the target of the present invention under horizontal polarization；
Fig. 3 is simulation result diagram of the target of the present invention under vertical polarization.
Specific implementation mode
The embodiment of the present invention and effect are described further below in conjunction with the accompanying drawings：
The targeted airborne radar of the present invention is positive side view imaging pattern, and the operating mode of carrier aircraft is " step one is stopped ", false
If carrier aircraft flight environment of vehicle is ideal situation, without abnormal shake when flight, target scene model is small triangle combination one by one
It forms.
Referring to Fig.1, it is the high-resolution polarimetric SAR target image rapid simulation method of the present invention, includes the following steps：
Step 1：Setting radar parameter simultaneously imports object module.
1.1) setting airborne radar height h, carrier aircraft flying speed v, light minimum irradiation angle, θmin, azimuth resolution
Δ A, range resolution Δ R, radar transmitting wave frequency ω, target surface permittivity ε, horizontal polarization directions electromagnetism
Wave amplitude Eih, vertical polarization directions electromagnetism wave amplitude Eiv；
1.2) by among existing target max format models importing business software 3Dmax, the model of obj formats is exported, it should
The model of obj formats be by multiple triangle sets at；
1.3) object module of obj formats is imported in business software matlab, extracts object module Triangular object model vertex
Coordinates matrix T.
Step 2：Calculate the parameter of object module.
2.1) compare in coordinates matrix T indicate orientation numerical value size, obtain object module orientation maximum value
AmaxWith the minimum value A in orientationmin；
2.2) compare in coordinates matrix T expression size of the distance to numerical value, obtain object module distance to maximum value
dmaxWith distance to minimum value dmin, according to dmaxAnd dminCalculate storage radar slant-range maximum value Rmax, minimum value RminRespectively：
Step 3：The orientation and distance of calculating radar form the light matrix of a K × N to sampled point number.
3.1) according to object module obtained by (2.1) in the maximum value A of orientationmaxWith the minimum value A in orientationmin, meter
Calculating sampling number K is：
3.2) according to radar slant-range maximum value A obtained by (2.2)maxWith radar slant-range minimum value Amin, distance is calculated to sampling
Points N is：
3.3) control radar does the linear uniform motion that speed is v, constantly emits beam during the motion, and the light is logical
It crosses in radar illumination to target scene, forms the light matrix of a K × N.
Step 4：Calculate radar site and light direction of illumination.
The position coordinates P { l, m } of radar and each light during imaging simulation are calculated according to the parameter in step 1 to shine
The direction r { l, m } penetrated：
P { l, m }=[- htan (θmin),Amax(l-1) Δ A, h],
R { l, m }=[sin (θm),0,-cos(θm)],
Wherein l is the line number in light matrix, and m is the columns in light matrix, θmFor m row light irradiating angle,θm-1For the light irradiating angle of m-1 row, θ1For the light irradiating angle of first row,
This example takes θ1=θmin。
Step 5：Illuminated face element is determined using the parameter in step 1.
5.1) face element number ζ on opticpath is calculated：
Make a ray along direction r { l, m } by radar site P { l, m }, seeks light EiWith the friendship of plane where Triangular object model
Point (x, y, z)：
vp1、vp2、vp3Three numerical value of the normal vector of plane, v where indicating Triangular object model respectively1, v2, v3Respectively three
Three vertex of edged surface member, v1(1)、v1(2)、v1(3) v is indicated respectively1Three coordinate values, p (1, l, m), p (2, l, m), p (3,
L, m) respectively indicate p { l, m } three coordinate values, r (1, l, m), r (2, l, m), r (3, l, m) indicate three of r { l, m } respectively
Coordinate value；
5.2) judge intersection point whether in the inside of Triangular object model：
If intersection point, inside Triangular object model, which takes the ζ values to be on the propagation path of this light1；
If without intersection point or intersection point not inside Triangular object model, the face element is not in the propagation path of this light
On, the value of ζ is 0；
5.3) judge whether light intersects with all Triangular object models successively according to (5.1) and (5.2), ζ values obtained by counting
Number for 1 is the number ζ of required face element；
5.4) according to the number ζ of target face element, the distance between radar and each face element on propagation path D are calculatedi：
Di=| (p { l, m }-(x, y, z)) |, 1≤i≤ζ；
5.5) size for comparing each face element distance in (5.4), by minimum value assignment to apart from minimum value R, the minimum range
Corresponding face element is the face element being irradiated to, and obtains three apex coordinates of face element, respectively v1, v2, v3, these three coordinates exist
The q rows of coordinates matrix.
Step 6：Calculate the normal vector for being irradiated to face element.
According to three vertex vs of (5.5) intermediate cam face element1, v2, v3, calculate face element normal vector g：
Step 7：Calculate the angle α between light and face element normal vector.
According to face element normal vector in light direction of illumination r { l, m } in step 4 and (6.2), light and face element normal direction are calculated
The angle α of vector：
Since its angle is acute angle, ifThen α=π-α.
Step 8：Calculated level polarizes and the Fresnel reflection coefficient under vertical polarization.
According to the angle α in target surface permittivity ε in step 1 and step 7 between light and face element normal vector, meter
Calculate Fresnel reflection coefficient R under horizontal polarization1With Fresnel reflection coefficient R under vertical polarization2：
Step 9：Calculate the back scattering matrix S of scattering,single1With the back scattering matrix S of rescattering2。
According to angle α, horizontal Fresnel reflection coefficient R between light and face element normal vector1, vertical Fresnel reflection system
Number R2, calculate the back scattering matrix S of scattering,single1With the back scattering matrix S of rescattering2：
Wherein：
S1vh=S1hv,
S2hh=m1[-2R1cosγcos2β]+n1[sin2γsin2β+R2sin2β(1+cos2γ)],
S2hv=2n1R2cosγcos2β+m1[sin2γsin2β+R1sin2β(1+cos2γ)],
S2vh=S2hv,
S2vv=2n2R1cosγcos2β+m2[sin2γsin2β+R2sin2β(1+cos2γ)],
m1=R2cos2γcos2β-R1sin2β
m2=-(R1+R2)cosγcosβsinβ
n1=-(R1+R2)cosγcosβsinβ
n2=-R1sin2β+R2cos2γcos2β
β=α,
a1、a2、a3Indicate that three coordinate values of the normal vector g of face element, β indicate radar motion flight path and triangular facet respectively
The angle of member, γ indicate antenna to the visual angle at ground trace center.
Step 10：Calculate the scattering energy of scattering,single and rescattering.
According to light E in (5.1)i, scattering,single collision matrix S in step 91With rescattering matrix S2, calculate single and dissipate
Penetrate ENERGY E1s(l, m) and rescattering ENERGY E2s(l,m)：
E1s(l, m)=S1·Ei,
E2s(l, m)=S2·Ei,
R { l, m } is light direction of illumination, and R is apart from minimum value.
Step 11：Calculate the scattering,single energy and rescattering energy of all light.
The scattering,single ENERGY E of all K × N light is calculated according to step 5- steps 101s(l, m) and rescattering energy
E2s(l, m) obtains scattering,single energy matrix E1sWith rescattering energy matrix E2s。
Step 12：Draw polarimetric SAR image.
By obtained scattering,single energy matrix E in step 111sWith rescattering energy matrix E2sIt is added, obtains total
Scatter energy matrix Es, by total scattering energy matrix EsAs picture element matrix, high-resolution polarimetric is drawn using business software matlab
SAR image completes above-mentioned emulation.
Below by a simulation example, invention is further explained：
1) experiment parameter is set：
Airborne radar height h is 5000m, and carrier aircraft speed v is 100m/s, minimum irradiation angle, θminIt is 60 degree, emits electricity
Magnetic wave frequency is 10GHz, and azimuth resolution Δ A is 0.1 meter, and range resolution Δ R is 0.1 meter, and target surface dielectric is normal
Number ε is 50, when emission level line polarization wave, horizontal polarization directions amplitude Eih=100m, vertical polarization directions amplitude Eiv=0；
When emitting vertical line polarization wave, horizontal polarization directions amplitude Eih=0, vertical polarization directions amplitude Eiv=100m.
2) simulating scenes
Target scene be T95 tank models in a larger plane, model of place size be 17.957m ×
19.771m×3.857m。
3) emulation content
To above-mentioned scene, polarization SAR imaging is carried out to it with basic inventive method.When transmitted wave is horizontal polarized wave,
The results are shown in Figure 2；When transmitted wave is vertically polarized wave, the results are shown in Figure 3.
By Fig. 2 and Fig. 3 as it can be seen that the present invention can accurately be imaged target scene, it can really reflect object
Position, shape these information, wherein Fig. 3 can also reflect that the polarization characteristic of target, the shade in Fig. 3 also reflect its authenticity
And reliability.
Above description is only example of the present invention, does not constitute any limitation of the invention, it is clear that for
It, all may be without departing substantially from the principle of the invention, structure after having understood the content of present invention and principle for one of skill in the art
In the case of, carry out various modifications and variations in form and details, but these modifications and variations based on inventive concept
Still within the claims of the present invention.
Claims (10)
1. high-resolution polarimetric SAR target image rapid simulation method, including：
(1) parameter setting
Airborne radar height h, carrier aircraft flying speed v, light minimum irradiation angle, θ are setmin, azimuth resolution Δ A, distance
To resolution ax R, radar transmitting wave frequency ω, target surface permittivity ε, horizontal polarization directions electromagnetism wave amplitude Eih
With vertical polarization directions electromagnetism wave amplitude Eiv；
(2) the target scene 3D models combined by Triangular object model are built, and are conducted into business software matlab, the 3D is extracted
The coordinates matrix T on model Triangular object model vertex obtains the orientation maximum value of radar by apex coordinate matrix T and radar altitude h
Amax, minimum value Amin, oblique distance maximum value RmaxWith oblique distance minimum value Rmin；
(3) parameter obtained according to (2), calculate radar bearing to sampling number K and distance to sampling number N；
(4) control radar does the linear uniform motion that speed is v, constantly emits beam during the motion, which passes through thunder
Up to being irradiated in target scene, the light matrix of a K × N is formed；
(5) according to the parameter being arranged in (1), during radar motion each light direction of illumination r of position P { l, m } l,
M }, wherein l is the line number in matrix, and 1≤l≤K, m are matrix column number, 1≤m≤N；
(6) Triangular object model being irradiated to is determined：
(6a) according to the parameter in (1) and (5), calculating light matrix line number is l, the light E that columns is miTarget on propagation path
The number ζ of face element, and calculate the distance between radar and each target face element on propagation path Di；
(6b) from small to large sorts a distance, by minimum value assignment to apart from minimum value R, the face corresponding to the minimum range
Member is the face element being irradiated to, if the face element coordinate is in the q rows of apex coordinate matrix T, 1≤q≤K；
(7) following parameter is calculated according to Triangular object model：
(7a) calculates face element normal vector g according to line number q where Triangular object model；
(7b) is calculated according to light direction of illumination r { l, m } and face element normal vector g between light and Triangular object model normal vector
Angle α；
(7c) is calculated according to the angle α in target surface permittivity ε in (1) and (7b) between light and face element normal vector
The Fresnel reflection coefficient R of horizontal polarization1With the Fresnel reflection coefficient R of vertical polarization2；
Horizontal Fresnel reflection coefficient R in angle α, (7c) in (7d) basis (7b) between light and face element normal vector1With hang down
Straight Fresnel reflection coefficient R2, calculate the back scattering matrix S of the scattering,single of face element1With the back scattering matrix of rescattering
S2；
(7e) is according to light E in (6)iThe back scattering matrix S of scattering,single in (7d)1, calculate the scattering,single energy of face element
E1s(l,m)；
(7f) is according to incident ray EiWith rescattering matrix S2, calculate the rescattering ENERGY E of face element2s(l,m)；
(8) step (6)-(7) are pressed, the scattering,single ENERGY E of all K × N light is calculated1sWith rescattering ENERGY E2s；
(9) by scattering,single ENERGY E1sWith rescattering ENERGY E2sThe sum of as pixel matrix Es, utilize business software
Matlab draws high-resolution polarimetric SAR image.
2. the method as described in claim 1, it is characterised in that：Calculated in step (3) radar bearing to sampling number K and away from
The sampling number N of descriscent, carries out as follows：
3. the method as described in claim 1, it is characterised in that：Radar motion is calculated in step 5 in the process at position P { l, m }
Each light direction of illumination r { l, m }, be calculated as follows：
P { l, m }=[- htan (θmin),Amax(l-1) Δ A, h],
R { l, m }=[sin (θm),0,-cos(θm)],
Wherein, θmFor m row light irradiating angle,
θm-1For m-1 row light irradiating angle,
θ1For the light irradiating angle of first row, θ1=θmin。
4. the method as described in claim 1, it is characterised in that：Accounting Line Number is l in step (6a), and columns is the light E of miIt passes
The number ζ of target face element on path is broadcast, and calculates the distance between radar and each face element on propagation path Ri, as follows into
Row：
(6a1) makees a ray by radar site P { l, m } along direction r { l, m }, asks the friendship of light and plane where Triangular object model
Point (x, y, z)：
vp1、vp2、vp3Three numerical value of the normal vector of plane, v where indicating Triangular object model respectively1(1)、v1(2)、v1(3) respectively
Indicate v1Three coordinate values, p (1, l, m), p (2, l, m), p (3, l, m) indicate three coordinate values of p { l, m } respectively, r (1,
L, m), r (2, l, m), r (3, l, m) indicate three coordinate values of r { l, m } respectively；
Whether (6a2) judges intersection point in the inside of Triangular object model：
If intersection point, inside Triangular object model, for the face element on the propagation path of this light, the value of ζ is then 1,
If without intersection point or intersection point not inside Triangular object model, the face element not on the propagation path of this light, ζ's
Value is then 0；
(6a3) judges whether light intersects with all Triangular object models successively according to (6a1) and (6a2) described step, statistics gained
The number for being 1 to ζ values is the number ζ of required face element；
The distance between radar and each face element on propagation path D is calculated as follows according to the number ζ of target face element in (6a4)i：
Di=| (p { l, m }-(x, y, z)) |, 1≤i≤ζ.
5. the method as described in claim 1, it is characterised in that：According to line number q where Triangular object model in step (7a), face is calculated
First normal vector g, is calculated as follows：
Wherein v1, v2, v3Respectively three vertex of vertex matrix q row Triangular object models.
6. the method as described in claim 1, it is characterised in that：According to light direction of illumination r { l, m } and face element in step (7b)
The angle α between light and Triangular object model normal vector is calculated as follows in normal vector g：
r(:, l, m) indicate radiation direction r { l, m } coordinate value.
7. the method as described in claim 1, it is characterised in that：According to target permittivity ε and light and method in step (7c)
To the angle α between vector, calculated level polarization Fresnel reflection coefficient R1With vertical polarization Fresnel reflection coefficient R2, by such as
Lower formula calculates：
8. the method as described in claim 1, it is characterised in that：It is pressed from both sides according between light and face element normal vector in step (7d)
Angle α, horizontal Fresnel reflection coefficient R1With vertical Fresnel reflection coefficient R2, calculate the back scattering matrix of face element scattering,single
S1With the back scattering matrix S of rescattering2, it is calculated as follows：
Wherein
S1vh=S1hv
S2hh=m1[-2R1cosγcos2β]+n1[sin2γsin2β+R2sin2β(1+cos2γ)]
S2hv=2n1R2cosγcos2β+m1[sin2γsin2β+R1sin2β(1+cos2γ)]
S2vh=S2hv
S2vv=2n2R1cosγcos2β+m2[sin2γsin2β+R2sin2β(1+cos2γ)]
m1=R2cos2γcos2β-R1sin2β
m2=-(R1+R2)cosγcosβsinβ
n1=-(R1+R2)cosγcosβsinβ
n2=-R1sin2β+R2cos2γcos2β
β=α
(a1,a2,a3) indicate face element normal vector g coordinate value, β indicate radar motion flight path and Triangular object model angle, γ
Visual angle of the expression antenna to ground trace center.
9. the method as described in claim 1, it is characterised in that：According to light E in step (7e)iBackward with scattering,single dissipates
Penetrate matrix S1, calculate the scattering,single ENERGY E of face element1s(l, m) is calculated as follows：
E1s(l, m)=S1·Ei,
Wherein
R { l, m } is light direction of illumination, and R is apart from minimum value.
10. the method as described in claim 1, it is characterised in that：According to incident ray E in step (7f)iWith rescattering matrix
S2, calculate the rescattering ENERGY E of face element2s(l, m) is calculated as follows：
E2s(l, m)=S2·Ei,
Wherein
R { l, m } is light direction of illumination, and R is apart from minimum value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810362394.3A CN108562899B (en) | 2018-04-20 | 2018-04-20 | High-resolution polarization SAR target image rapid simulation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810362394.3A CN108562899B (en) | 2018-04-20 | 2018-04-20 | High-resolution polarization SAR target image rapid simulation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108562899A true CN108562899A (en) | 2018-09-21 |
| CN108562899B CN108562899B (en) | 2021-08-20 |
Family
ID=63536081
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810362394.3A Active CN108562899B (en) | 2018-04-20 | 2018-04-20 | High-resolution polarization SAR target image rapid simulation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108562899B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111123215A (en) * | 2019-11-11 | 2020-05-08 | 南京长峰航天电子科技有限公司 | Multi-channel-based full polarization target implementation system and method |
| CN112531924A (en) * | 2020-10-13 | 2021-03-19 | 电子科技大学 | Method for rapidly designing antenna array based on maximum wireless energy transmission efficiency |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101236247A (en) * | 2008-03-07 | 2008-08-06 | 北京航空航天大学 | Star-carrying multichannel antenna SAR data channel amplitude and phase error correction platform |
| CN101329402A (en) * | 2008-06-20 | 2008-12-24 | 西安电子科技大学 | Multi-dimension SAR image edge detection method based on improved Wedgelet |
| US20090102705A1 (en) * | 2006-11-30 | 2009-04-23 | Obermeyer Henry K | Spectrometric synthetic aperture radar |
| CN101894368A (en) * | 2010-07-06 | 2010-11-24 | 西安电子科技大学 | Method for dividing semi-monitoring SAR image water area based on geodesic distance |
| CN102901965A (en) * | 2012-09-26 | 2013-01-30 | 同济大学 | Radar imaging method of electrically large size target in ocean clutter environment |
| CN103983968A (en) * | 2014-03-20 | 2014-08-13 | 西安电子科技大学 | Complete polarization type SAR super-resolution imaging method based on distributed compressed sensing |
| CN104360317A (en) * | 2014-12-09 | 2015-02-18 | 中国科学院遥感与数字地球研究所 | Method for making linear target all-azimuth angle SAR (Synthetic Aperture Radar) scene matching reference diagram |
| CN103792521B (en) * | 2014-01-17 | 2016-05-04 | 西安空间无线电技术研究所 | A kind of radar target back scattering analogy method that is subject to Faraday rotation variable effect |
| CN106680812A (en) * | 2016-10-27 | 2017-05-17 | 西安空间无线电技术研究所 | Microwave correlated imaging simulation method based on analysis surface element |
| CN107356924A (en) * | 2017-07-14 | 2017-11-17 | 中国科学院电子学研究所 | For the SAR image emulation mode of cavity class formation |
| CN107561535A (en) * | 2017-09-07 | 2018-01-09 | 中国科学院电子学研究所 | A kind of synthetic aperture radar range ambiguity suppressing method and device |
| CN107636693A (en) * | 2015-03-20 | 2018-01-26 | 弗劳恩霍夫应用研究促进协会 | Assigned for the relevance scores of artificial neural network |
-
2018
- 2018-04-20 CN CN201810362394.3A patent/CN108562899B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090102705A1 (en) * | 2006-11-30 | 2009-04-23 | Obermeyer Henry K | Spectrometric synthetic aperture radar |
| CN101236247A (en) * | 2008-03-07 | 2008-08-06 | 北京航空航天大学 | Star-carrying multichannel antenna SAR data channel amplitude and phase error correction platform |
| CN101329402A (en) * | 2008-06-20 | 2008-12-24 | 西安电子科技大学 | Multi-dimension SAR image edge detection method based on improved Wedgelet |
| CN101894368A (en) * | 2010-07-06 | 2010-11-24 | 西安电子科技大学 | Method for dividing semi-monitoring SAR image water area based on geodesic distance |
| CN102901965A (en) * | 2012-09-26 | 2013-01-30 | 同济大学 | Radar imaging method of electrically large size target in ocean clutter environment |
| CN103792521B (en) * | 2014-01-17 | 2016-05-04 | 西安空间无线电技术研究所 | A kind of radar target back scattering analogy method that is subject to Faraday rotation variable effect |
| CN103983968A (en) * | 2014-03-20 | 2014-08-13 | 西安电子科技大学 | Complete polarization type SAR super-resolution imaging method based on distributed compressed sensing |
| CN104360317A (en) * | 2014-12-09 | 2015-02-18 | 中国科学院遥感与数字地球研究所 | Method for making linear target all-azimuth angle SAR (Synthetic Aperture Radar) scene matching reference diagram |
| CN107636693A (en) * | 2015-03-20 | 2018-01-26 | 弗劳恩霍夫应用研究促进协会 | Assigned for the relevance scores of artificial neural network |
| CN106680812A (en) * | 2016-10-27 | 2017-05-17 | 西安空间无线电技术研究所 | Microwave correlated imaging simulation method based on analysis surface element |
| CN107356924A (en) * | 2017-07-14 | 2017-11-17 | 中国科学院电子学研究所 | For the SAR image emulation mode of cavity class formation |
| CN107561535A (en) * | 2017-09-07 | 2018-01-09 | 中国科学院电子学研究所 | A kind of synthetic aperture radar range ambiguity suppressing method and device |
Non-Patent Citations (2)
| Title |
|---|
| FENG XU ET.AL: "《Imaging Simulation of polarimetric SAR for a Comprehensive Terrain Scene Using the Mapping and Projection Algorithm》", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 * |
| 王国军 等: "《建筑物目标高分辨率SAR图像模拟方法》", 《遥感技术与应用》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111123215A (en) * | 2019-11-11 | 2020-05-08 | 南京长峰航天电子科技有限公司 | Multi-channel-based full polarization target implementation system and method |
| CN111123215B (en) * | 2019-11-11 | 2021-11-05 | 南京长峰航天电子科技有限公司 | Multi-channel-based full polarization target implementation system and method |
| CN112531924A (en) * | 2020-10-13 | 2021-03-19 | 电子科技大学 | Method for rapidly designing antenna array based on maximum wireless energy transmission efficiency |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108562899B (en) | 2021-08-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104865562B (en) | Identification method for radar disoperative target based on mixed model | |
| Schlager et al. | State-of-the-art sensor models for virtual testing of advanced driver assistance systems/autonomous driving functions | |
| KR101809067B1 (en) | Determination of mobile display position and orientation using micropower impulse radar | |
| CN101923166B (en) | Method for simulating composite hole diameter radar image of three-dimensional complex object | |
| CN104050716B (en) | A kind of marine multiple target SAR image visual modeling method | |
| CN109471080A (en) | High speed platform radar simulated radar echo system based on simulink | |
| CN101281249A (en) | Method for modeling target dispersion characteristic of high resolution synthetic aperture radar | |
| CN106324571A (en) | Fast Realization method for simulation 3D scene SAR radar echo based on forward method | |
| CN106872978A (en) | A kind of Electromagnetic Modeling emulation mode of complex scene | |
| Kleinkort et al. | Visual hull method for realistic 3D particle shape reconstruction based on high-resolution photographs of snowflakes in free fall from multiple views | |
| CN106770045A (en) | The simple calculating method of target laser scattering properties under a kind of local irradiation | |
| CN202735513U (en) | Holographic active microwave imaging device | |
| CN107607951A (en) | A kind of SAR image rescattering characteristic simulation method | |
| CN108761397A (en) | Polarization SAR model decomposition evaluation method based on electromagnetic scattering simulation | |
| CN108562899A (en) | High-resolution polarimetric SAR target image rapid simulation method | |
| CN110471061A (en) | A kind of emulation mode and its system for realizing airborne synthetic aperture radar imaging | |
| Kusk et al. | Synthetic SAR image generation using sensor, terrain and target models | |
| CN109188384A (en) | The Electromagnetic Simulation method of extraterrestrial target echo dynamic observation | |
| CN107024693A (en) | A kind of radar relevance imaging method of single-shot beam | |
| CN117269940A (en) | Point cloud data generation method and perception capability verification method of laser radar | |
| Zherdev et al. | Object recognition using real and modelled SAR images | |
| Chen et al. | SAR imaging simulation for urban structures based on analytical models | |
| Li et al. | A novel DEM extraction method based on chain correlation of CSAR subaperture images | |
| Cheng et al. | M 3 SC: A generic dataset for mixed multi-modal (MMM) sensing and communication integration | |
| CN111522007A (en) | SAR imaging simulation method and system with real scene and target simulation fused |
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 |