CN103487808B - A kind of track Simulation method of variable element locking mode missile-borne Spotlight SAR Imaging - Google Patents
A kind of track Simulation method of variable element locking mode missile-borne Spotlight SAR Imaging Download PDFInfo
- Publication number
- CN103487808B CN103487808B CN201310413667.XA CN201310413667A CN103487808B CN 103487808 B CN103487808 B CN 103487808B CN 201310413667 A CN201310413667 A CN 201310413667A CN 103487808 B CN103487808 B CN 103487808B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- flight path
- sar
- plane
- coordinate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
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/9052—Spotlight mode
-
- 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
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 variable element locking mode missile-borne Spotlight SAR Imaging track Simulation method, comprise step one, determine four coordinate systems of employing in emulating, step 2, obtain the transformed matrix that coordinate system asks, step 3, setting radar system fundamental simulation parameter, step 4, utilize runge kutta method to emulate to obtain body at flight path plane coordinate system E
min thick emulation flight path coordinates matrix, step 5, to obtain in the synthetic aperture time SAR at earth rotation coordinate system E
gin thick emulation flight path coordinates matrix, step 6, obtain accurately body emulation flight path.The present invention proposes a kind of locking mode of operation of variable element missile-borne Spotlight SAR Imaging, the distinctive fixed strabismus angle of this pattern is that body guidance is provided convenience.
Description
Technical field
The present invention relates to synthetic-aperture radar (Synthetic Aperture Radar is called for short SAR) technical field, specifically say, refer to a kind of track Simulation method of variable element locking mode missile-borne Spotlight SAR Imaging.
Background technology
Synthetic-aperture radar is that a kind of signal analysis technology that relies on builds equivalent vast of heaven line, to obtain beam elevation upwards high-resolution round-the-clock two-dimensional imaging radar.The beam bunching mode synthetic-aperture radar of variable element system is for the high precision mapping of large stravismus zonule, the application demand of investigation and proposing, its by adjust carrier frequency of SAR system, frequency modulation rate, pulse width, sampling rate and pulse repeat to ask every, be provided with following advantage: compared to the SAR system of the constant parameter of routine, variable element SAR has the imaging processing flow process more simplified, higher image processing speed, the space availability ratio of data in higher SAR data storer, thus total amount of data is reduced, the volume of transmitted data that reduction SAR and ground base station are asked.
Variable element system SAR is particularly applicable to being applied in Missile-borne SAR guidance.In Missile-borne SAR guidance task, SAR often needs to work in large stravismus or super large strabismus mode.Different and airborne or satellite-borne SAR, for evading interception, the movement locus of Missile-borne SAR is not straight line usually, under missile-borne locking beam bunching mode, the direction of motion of body with play angle that order oblique distance asks in imaging process, need to be maintained a fixing angle, the flight path plane of body also no longer with ground level keeping parallelism.This demand is the difficulties in the track Simulation of variable element locking mode missile-borne Spotlight SAR Imaging.
Summary of the invention
The object of the invention is to solve the problem to solve, for variable element locking mode missile-borne Spotlight SAR Imaging, proposing a kind of track Simulation method.The flight path of setting body is restrained in three-dimensional single plane, angle is being locked according to known, fourth-order Runge-Kutta method is adopted to generate the flight path of thick emulation, and then according to the sampling location of the accurate flight path of variable element SAR system design Missile-borne SAR, utilize the method for interpolation to generate the flight path of body emulation accurately.
A kind of variable element locking mode missile-borne Spotlight SAR Imaging track Simulation method of the present invention, comprises following step:
Step one, determine emulate in adopt four coordinate systems: earth rotation coordinate system E
g, imaging scene coordinate system E
s, direction of visual lines polar coordinate system E
lwith body flight path plane coordinate system E
m.
In step 2, determining step one, four coordinate systems ask the method for transformation of coordinate, obtain the transformed matrix that coordinate system is asked.
Step 3, setting radar system fundamental simulation parameter, determine the accumulation corner Δ θ of SAR in the synthetic aperture time, determine the quantity N of the sampling pulse launched along flight path direction SAR
a.
The time step t of step 4, the emulation of setting runge kutta method
h.Utilize runge kutta method to emulate and obtain body at flight path plane coordinate system E
min thick emulation flight path coordinates matrix.
Step 5, integrating step four and the result in step 2, to obtain in the synthetic aperture time SAR at earth rotation coordinate system E after coordinate transformation
gin thick emulation flight path coordinates matrix.
Step 6, in conjunction with variable element SAR system design interpolation method, interpolation is carried out to the thick emulation flight path obtained in step 5, obtain the flight path of body emulation accurately.
The invention has the advantages that:
(1) propose a kind of locking mode of operation of variable element missile-borne Spotlight SAR Imaging, the distinctive fixed strabismus angle of this pattern is that body guidance is provided convenience;
(2) utilize runge kutta method to carry out the recursion calculation of body flight path, carry out interpolation processing in conjunction with variable element system, obtain the accurate flight path of locking mode missile-borne Spotlight SAR Imaging;
(3) this method is applicable to the body track Simulation of single flight path plane, and imaging scene center can not be arctic point.
Accompanying drawing explanation
Fig. 1 is method flow diagram of the present invention;
Fig. 2 is the schematic diagram of four coordinate systems of the present invention;
Fig. 3 is along orientation to the schematic diagram data be spacedly distributed in wavenumber domain of the present invention;
Fig. 4 is the flight path schematic diagram of gained body of the present invention under earth rotation coordinate system.
Embodiment
Below in conjunction with drawings and Examples, the present invention is described in further detail.
The present invention is a kind of track Simulation method of variable element locking mode missile-borne Spotlight SAR Imaging, and process flow diagram as shown in Figure 1, comprises following step:
Step one, determine emulate in adopt four coordinate systems: earth rotation coordinate system E
g, imaging scene coordinate system E
s, direction of visual lines polar coordinate system E
lwith body flight path plane coordinate system E
m, four coordinate systems as shown in Figure 2.Be specially:
Set with the earth tangent plane of scene center formation as the flight path plane of L, SAR is for P.
(1) earth rotation coordinate system E
g
True origin: the earth's core is designated as O
g;
X-axis is designated as X
g: under the line in plane, point to zero degree warp direction;
Y-axis is designated as Y
g: under the line in plane, point to east longitude 90 degree of warp directions;
Z axis is designated as Z
g: along earth's axis, point to the positive arctic (north latitude 90 degree of directions).
(2) imaging scene coordinate system E
s
True origin: scene center is designated as O
s;
X-axis is designated as X
s: in plane L, point to the east of scene;
Y-axis is designated as Y
s: in plane L, point to the north of scene;
Z axis is designated as Z
s: cross scene central point, deviate from direction, the earth's core perpendicular to the earth tangent plane formed with scene center.
(3) the polar coordinate system E of direction of visual lines
l
True origin is designated as O
l: scene center;
X-axis is designated as X
l: in plane L, with Y
laxle, Z
laxle forms right-handed coordinate system;
Y-axis is designated as Y
l: in plane L, point to SAR sensor along the projection of synthetic aperture central instant oblique distance in this plane;
Z axis is designated as Z
l: cross scene central point, deviate from direction, the earth's core perpendicular to the earth tangent plane formed with scene center.
(4) body flight path plane coordinate system E
m
True origin is designated as O
m: scene center;
X-axis is designated as X
m: in body flight path plane, with Y
maxle, Z
maxle forms right-handed coordinate system;
Y-axis is designated as Y
m: in body flight path plane, point to the direction of target and SAR sensor;
Z axis is designated as Z
m: deviate from direction, the earth's core, directed in orthogonal is in the direction of body flight path plane P.
In step 2, determining step one, four coordinate systems ask the method for transformation of coordinate, obtain the transformed matrix that coordinate system is asked.
If SAR is at earth rotation coordinate system E
gin synthetic aperture central instant (be designated as t
0) coordinate is (x
g0, y
g0, z
g0), at imaging scene coordinate system E
smiddle t
0the coordinate in moment is (x
s0, y
s0, z
s0), scene center is at E
gcoordinate (x in coordinate system
p, y
p, z
p).Certain moment SAR is at earth rotation coordinate system E
gin coordinate be (x
g, y
g, z
g), SAR is at imaging scene coordinate system E in the same time
sin coordinate be (x
s, y
s, z
s), at the polar coordinate system E of direction of visual lines
lin coordinate be (x
l, y
l, z
l), at body flight path plane coordinate system E
min coordinate be (x
m, y
m, z
m).
(1) imaging scene coordinate system E is obtained
sto earth rotation coordinate system E
gtransition matrix A
sg, earth rotation coordinate system E
gto imaging scene coordinate system E
stransition matrix A
gs.
The longitude φ of SAR at synthetic aperture central instant is obtained according to formula (1) and formula (2)
lowith latitude φ
la:
A
sgmeet the relation that formula (3) represents, can A be obtained according to formula (4)
sg
A
sga is obtained through inversion operation
gs.
(2) imaging scene coordinate system E is obtained
sto the polar coordinate system E of direction of visual lines
ltransition matrix A
sl, the polar coordinate system E of direction of visual lines
lto imaging scene coordinate system E
stransition matrix A
ls.
Obtain synthetic aperture central instant SAR according to formula (5) and depart from imaging scene coordinate system E
sthe angle of Y direction is α, then
A
slmeet the relation that formula (6) represents, can A be obtained according to formula (7)
sl
[x
ly
lz
l]=[x
sy
sz
s]A
sl(6)
A
sla is obtained through inversion operation
ls.
(3) the polar coordinate system E of direction of visual lines is obtained
lto body flight path plane coordinate system E
mtransition matrix A
lm, body flight path plane coordinate system E
mto the polar coordinate system E of direction of visual lines
ltransition matrix A
ml.
Obtaining synthetic aperture central instant SAR to the incident angle of scene center according to formula (8) is β, then
Setting flight path plane P and plane L angle are γ, obtain flight path plane relative rotation angle Ω to be according to formula (9)
A
lmmeet the relation that formula (10) represents, can A be obtained according to formula (11)
sl
A
lma is obtained through inversion operation
ml.
Step 3, setting radar system fundamental simulation parameter, determine the accumulation corner Δ θ of SAR in the synthetic aperture time, determine the quantity N of the sampling pulse launched along flight path direction SAR
a.
(1) the reference carrier frequency f of radar system is set
0, light velocity c, simulating scenes radius D, duration of pulse μ, ρ
a, ρ
rbe respectively azimuth resolution and the range resolution of system expection, obtaining according to formula (12) corner that in the synthetic aperture time, SAR accumulates in plane L is Δ θ, then
(2) set function nextpow2 () as rreturn value be greater than independent variable minimum 2 the function of integer corresponding to integral number power, obtained the quantity N of the sampling pulse launched along flight path direction SAR by formula (13)
a
The time step t of step 4, the emulation of setting runge kutta method
h.Utilize runge kutta method to emulate and obtain body at flight path plane coordinate system E
min sample coordinate matrix Posm.
(1) setting projectile flight speed is V, Texas tower velocity reversal and the locking angle θ played between order oblique distance, the time step t of runge kutta method emulation
h.If synthetic aperture initial time is t
1, finish time is t
2, t
1< t
0< t
2, utilize runge kutta method to generate t
0~ t
2the thick emulation in moment slightly emulates flight path.If (x
f, y
f, 0) and be t
mmoment SAR is at body flight path plane coordinate system E
min coordinate, (x
1f, y
1f, 0) and be t
m+ t
hmoment SAR is at body flight path plane coordinate system E
min coordinate.If x
fy
finitial value be that synthetic aperture central instant SAR platform is at body flight path plane coordinate system E
min coordinate.Intermediate parameters k is calculated according to formula (14)-formula (18)
1xf, k
1yf, k
2xf, k
2yf, k
3xf, k
3yf, k
4xf, k
4yf.
(2) repeat (1) of step 4, utilize runge kutta method to generate t
1~ t
0the thick emulation flight path of moment SAR is at body flight path plane coordinate system E
min coordinate.(1) and (2) of combining step four, to obtain in the synthetic aperture time SAR at body flight path plane coordinate system E
min thick emulation flight path coordinates matrix, be expressed as matrix Posm.
Step 5, integrating step four and the result in step 2, to obtain in the synthetic aperture time SAR at earth rotation coordinate system E after coordinate transformation
gin thick emulation flight path coordinates matrix.
Poss is that this flight path is at earth rotation coordinate system E
gin flight path coordinates matrix.In integrating step two, coordinate system transforms rule, obtains Poss according to formula (19)
Poss=PosmA
mlA
lsA
sg(19)
Step 6, in conjunction with variable element SAR system design interpolation method, interpolation is carried out to the thick emulation flight path obtained in step 5, obtain the flight path of body emulation accurately.
Application claims SAR echo data at wavenumber domain along orientation to being spacedly distributed, as shown in Figure 3.Be located at imaging scene coordinate system E
sin, the angle sequence that SAR departs from Y-axis in plane L is Δ ε, then
Flight path slightly can be emulated at imaging scene coordinate system E by matrix Poss
sin, the angle sequence that SAR departs from Y-axis in plane L is Δ ε
0.With Δ ε
0for raw data position, Poss is raw data value, and Δ ε is Data Position after interpolation, and after obtaining interpolation after carrying out interpolation operation, SAR is at imaging scene coordinate system E
sin coordinate Possl accurately, the SAR flight path namely obtained after the present invention's emulation.
embodiment
Variable element locking mode missile-borne Spotlight SAR Imaging track Simulation performance parameter, systematic parameter, Space geometric parameter, derived parameter are as shown in table 1,
Table 1 variable element locking mode missile-borne Spotlight SAR Imaging track Simulation parameter list
Step one, as Fig. 2 determine emulate in adopt four coordinate systems: earth rotation coordinate system E
g, imaging scene coordinate system E
s, direction of visual lines polar coordinate system E
lwith body flight path plane coordinate system E
m, shown in four coordinate systems.
Step 2, establish SAR at earth rotation coordinate system E
gthe coordinate of middle synthetic aperture central instant is (1253.3 ,-3392.7,5247.6) (km), at imaging scene coordinate system E
smiddle t
0the coordinate in moment is (0,11.491,9.642) (km), and scene center is at E
gcoordinate in coordinate system is (1240.6 ,-3391.5,5239.7) (km),
(1) imaging scene coordinate system E
sto earth rotation coordinate system E
stransition matrix A
sg, earth rotation coordinate system E
gto imaging scene coordinate system E
stransition matrix A
gs.
Base area spherical model can obtain the longitude φ of SAR at synthetic aperture central instant
lowith φ
la:
φ
l0=-69.91°
φ
la=55.42°
(2) imaging scene coordinate system E
sto the polar coordinate system E of direction of visual lines
ltransition matrix A
sl, the polar coordinate system E of direction of visual lines
lto imaging scene coordinate system E
stransition matrix A
ls.
If synthetic aperture central instant SAR departs from imaging scene coordinate system E
sangle [alpha]=0 of Y direction.
(3) the polar coordinate system E of direction of visual lines
lto body flight path plane coordinate system E
mtransition matrix A
lm, body flight path plane coordinate system E
mto the polar coordinate system E of direction of visual lines
ltransition matrix A
ml.
Setting synthetic aperture central instant SAR is β to the incident angle of scene center, and β=50 °, flight path plane relative rotation angle Ω=0 °.
Parameter in step 3, foundation table 1, calculates Δ θ=0.347 °, N according to formula (12) and (13)
a=1024.
The time step t of step 4, the emulation of setting runge kutta method
h=0.001s.Utilize runge kutta method to emulate and obtain body at flight path plane coordinate system E
min sample coordinate matrix Posm.
Step 5, integrating step four and the result in step 2, to obtain in the synthetic aperture time SAR at earth rotation coordinate system E after coordinate transformation
gin thick emulation flight path coordinates matrix.If Poss is that this flight path is at earth rotation coordinate system E
gin flight path coordinates matrix.In integrating step two, coordinate system transforms rule, obtains Poss.
Step 6, in conjunction with simulation parameter list 1, obtaining by formula (20) the angle sequence that SAR departs from Y-axis in plane L is Δ ε.Flight path slightly can be emulated at imaging scene coordinate system E by matrix Poss
sin, the angle sequence that SAR departs from Y-axis in plane L is Δ ε
0.Comprehensive Δ ε
0, Δ ε and Poss, the accurate SAR after interpolation can be obtained at imaging scene coordinate system E
sin coordinate represent Possl.Schematic diagram is as shown in Figure 4 in earth rotation coordinate system for the body flight path of emulation gained.
Process provides a kind of track Simulation method being applicable to the variable element locking mode missile-borne Spotlight SAR Imaging of single flight path plane, simulation result is that subsequent echoes emulation provides the foundation with imaging processing.
Claims (1)
1. a track Simulation method for variable element locking mode missile-borne Spotlight SAR Imaging, comprises following step:
Step one, determine emulate in adopt four coordinate systems: earth rotation coordinate system E
g, imaging scene coordinate system E
s, direction of visual lines polar coordinate system E
lwith body flight path plane coordinate system E
m, be specially:
Set with the earth tangent plane of scene center formation as the flight path plane of L, SAR is for P;
(1) earth rotation coordinate system E
g
True origin: the earth's core is designated as O
g;
X-axis is designated as X
g: under the line in plane, point to zero degree warp direction;
Y-axis is designated as Y
g: under the line in plane, point to east longitude 90 degree of warp directions;
Z axis is designated as Z
g: along earth's axis, point to the positive arctic;
(2) imaging scene coordinate system E
s
True origin: scene center is designated as O
s;
X-axis is designated as X
s: in plane L, point to the east of scene;
Y-axis is designated as Y
s: in plane L, point to the north of scene;
Z axis is designated as Z
s: cross scene central point, deviate from direction, the earth's core perpendicular to the earth tangent plane formed with scene center;
(3) the polar coordinate system E of direction of visual lines
l
True origin is designated as O
l: scene center;
X-axis is designated as X
l: in plane L, with Y
laxle, Z
laxle forms right-handed coordinate system;
Y-axis is designated as Y
l: in plane L, point to SAR sensor along the projection of synthetic aperture central instant oblique distance in this plane;
Z axis is designated as Z
l: cross scene central point, deviate from direction, the earth's core perpendicular to the earth tangent plane formed with scene center;
(4) body flight path plane coordinate system E
m
True origin is designated as O
m: scene center;
X-axis is designated as X
m: in body flight path plane, with Y
maxle, Z
maxle forms right-handed coordinate system;
Y-axis is designated as Y
m: in body flight path plane, point to the direction of target and SAR sensor;
Z axis is designated as Z
m: deviate from direction, the earth's core, directed in orthogonal is in the direction of body flight path plane P;
In step 2, determining step one, the method for transformation of coordinate between four coordinate systems, obtains the transformed matrix between coordinate system;
If SAR is at earth rotation coordinate system E
gin synthetic aperture central instant coordinate be (x
g0, y
g0, z
g0), at imaging scene coordinate system E
sin synthetic aperture central instant coordinate be (x
s0, y
s0, z
s0), scene center is at E
gcoordinate (x in coordinate system
p, y
p, z
p); Certain moment SAR is at earth rotation coordinate system E
gin coordinate be (x
g, y
g, z
g), SAR is at imaging scene coordinate system E in the same time
sin coordinate be (x
s, y
s, z
s), at the polar coordinate system E of direction of visual lines
lin coordinate be (x
l, y
l, z
l), at body flight path plane coordinate system E
min coordinate be (x
m, y
m, z
m);
(1) imaging scene coordinate system E is obtained
sto earth rotation coordinate system E
gtransition matrix A
sg, earth rotation coordinate system E
gto imaging scene coordinate system E
stransition matrix A
gs;
SAR is obtained at synthetic aperture central instant t according to formula (1) and formula (2)
0longitude φ
lowith latitude φ
la:
A
sgmeet the relation that formula (3) represents, can A be obtained according to formula (4)
sg
A
sga is obtained through inversion operation
gs;
(2) imaging scene coordinate system E is obtained
sto the polar coordinate system E of direction of visual lines
ltransition matrix A
sl, the polar coordinate system E of direction of visual lines
lto imaging scene coordinate system E
stransition matrix A
ls;
Obtain synthetic aperture central instant SAR according to formula (5) and depart from imaging scene coordinate system E
sthe angle of Y direction is α, then
A
slmeet the relation that formula (6) represents, can A be obtained according to formula (7)
sl
[x
ly
lz
l]=[x
sy
sz
s]A
sl(6)
A
sla is obtained through inversion operation
ls;
(3) the polar coordinate system E of direction of visual lines is obtained
lto body flight path plane coordinate system E
mtransition matrix A
lm, body flight path plane coordinate system E
mto the polar coordinate system E of direction of visual lines
ltransition matrix A
ml;
Obtaining synthetic aperture central instant SAR to the incident angle of scene center according to formula (8) is β, then
Setting flight path plane P and plane L angle are γ, obtain flight path plane relative rotation angle Ω to be according to formula (9)
A
lmmeet the relation that formula (10) represents, can A be obtained according to formula (11)
lm
[x
my
mz
m]=[x
ly
lz
l]A
lm(10)
A
lma is obtained through inversion operation
ml;
Step 3, setting radar system fundamental simulation parameter, determine the accumulation corner Δ θ of SAR in the synthetic aperture time, determine the quantity N of the sampling pulse launched along flight path direction SAR
a;
(1) the reference carrier frequency f of radar system is set
0, light velocity c, simulating scenes radius D, duration of pulse μ, ρ
afor the azimuth resolution of system expection, obtaining according to formula (12) corner that in the synthetic aperture time, SAR accumulates in plane L is Δ θ, then
(2) set function nextpow2 () as rreturn value be greater than independent variable minimum 2 the function of integer corresponding to integral number power, obtained the quantity N of the sampling pulse launched along flight path direction SAR by formula (13)
a
The time step t of step 4, the emulation of setting runge kutta method
h; Utilize runge kutta method to emulate and obtain body at flight path plane coordinate system E
min sample coordinate matrix Posm;
(1) setting projectile flight speed is V, Texas tower velocity reversal and the locking angle θ played between order oblique distance, the time step t of runge kutta method emulation
h; If synthetic aperture initial time is t
1, finish time is t
2, t
1<t
0<t
2, utilize runge kutta method to generate t
0~ t
2the thick emulation flight path in moment; If (x
f, y
f, 0) and be t
mmoment SAR is at body flight path plane coordinate system E
min coordinate, (x
1f, y
1f, 0) and be t
m+ t
hmoment SAR is at body flight path plane coordinate system E
min coordinate; If x
f, y
finitial value be that synthetic aperture central instant SAR platform is at body flight path plane coordinate system E
min coordinate; Intermediate parameters k is calculated according to formula (14)-formula (18)
1xf, k
1yf, k
2xf, k
2yf, k
3xf, k
3yf, k
4xf, k
4yf;
(2) repeat (1) of step 4, utilize runge kutta method to generate t
1~ t
0the thick emulation flight path of moment SAR is at body flight path plane coordinate system E
min coordinate; (1) and (2) of combining step four, to obtain in the synthetic aperture time SAR at body flight path plane coordinate system E
min thick emulation flight path coordinates matrix, be expressed as matrix Posm;
Step 5, integrating step four and the result in step 2, to obtain in the synthetic aperture time SAR at earth rotation coordinate system E after coordinate transformation
gin thick emulation flight path coordinates matrix;
Poss is that this flight path is at earth rotation coordinate system E
gin flight path coordinates matrix; In integrating step two, coordinate system transforms rule, obtains Poss according to formula (19)
Poss=PosmA
mlA
lsA
sg(19)
Step 6, in conjunction with variable element SAR system design interpolation method, interpolation is carried out to the thick emulation flight path obtained in step 5, obtain the flight path of body emulation accurately;
Application claims SAR echo data at wavenumber domain along orientation to being spacedly distributed;
Be located at imaging scene coordinate system E
sin, the angle sequence that SAR departs from Y-axis in plane L is Δ ε, then
Flight path slightly can be emulated at imaging scene coordinate system E by matrix Poss
sin, the angle sequence that SAR departs from Y-axis in plane L is Δ ε
0; With Δ ε
0for raw data position, Poss is raw data value, and Δ ε is Data Position after interpolation, and after obtaining interpolation after carrying out interpolation operation, SAR is at imaging scene coordinate system E
sin coordinate Poss1 accurately, the SAR flight path obtained after being emulation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310413667.XA CN103487808B (en) | 2013-09-12 | 2013-09-12 | A kind of track Simulation method of variable element locking mode missile-borne Spotlight SAR Imaging |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310413667.XA CN103487808B (en) | 2013-09-12 | 2013-09-12 | A kind of track Simulation method of variable element locking mode missile-borne Spotlight SAR Imaging |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103487808A CN103487808A (en) | 2014-01-01 |
CN103487808B true CN103487808B (en) | 2015-07-29 |
Family
ID=49828177
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310413667.XA Active CN103487808B (en) | 2013-09-12 | 2013-09-12 | A kind of track Simulation method of variable element locking mode missile-borne Spotlight SAR Imaging |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103487808B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103901430B (en) * | 2014-04-21 | 2016-11-09 | 西安电子科技大学 | Double-basis Forward-looking SAR dive section radar transmitter rail design method |
CN105549008B (en) * | 2015-12-02 | 2018-02-06 | 北京航空航天大学 | A kind of parameter optimization method of the big spaceborne Spotlight SAR Imaging system of strabismus of variable element |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102890276A (en) * | 2012-07-20 | 2013-01-23 | 北京航空航天大学 | Location method for missile-borne SAR (synthetic aperture radar) motion error image |
-
2013
- 2013-09-12 CN CN201310413667.XA patent/CN103487808B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102890276A (en) * | 2012-07-20 | 2013-01-23 | 北京航空航天大学 | Location method for missile-borne SAR (synthetic aperture radar) motion error image |
Non-Patent Citations (3)
Title |
---|
"A New Approach of Extended Chirp Scaling Algorithm for High Squint Missile-borne SAR Data Processing";Zhou Qiang et al.;《2008 International Symposium on Computer Science and Computational Technology》;20081231;第133-136页 * |
"弹载聚束SAR平台的PFA算法";彭岁阳 等;《宇航学报》;20110331;第32卷(第3期);第592-598页 * |
"环视SAR成像处理中的几何失真校正算法";李勇 等;《南京航空航天大学学报》;20090430;第41卷(第2期);第232-237页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103487808A (en) | 2014-01-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9007570B1 (en) | Airborne wind profiling algorithm for Doppler Wind LIDAR | |
CN101339244B (en) | On-board SAR image automatic target positioning method | |
CN102749621B (en) | Bistatic synthetic aperture radar (BSAR) frequency domain imaging method | |
CN103675775B (en) | Background ionosphere is to the analytical approach of GEO SAR Imaging | |
CN103197291B (en) | Satellite-borne synthetic aperture radar (SAR) echo signal simulation method based on non-stop walking model | |
EP2759847A1 (en) | Method and apparatus for determining equivalent velocity | |
CN102654406A (en) | Initial alignment method for moving bases based on combination of nonlinear prediction filtering and cubature Kalman filtering | |
CN106646395B (en) | A kind of radar return deduction method of airbound target | |
CN103487803A (en) | Airborne scanning radar imaging method in iteration compression mode | |
CN110823191B (en) | Method and system for determining ocean current measurement performance of mixed baseline dual-antenna squint interference SAR | |
Tang et al. | Focusing hypersonic vehicle-borne SAR data using radius/angle algorithm | |
CN104730506B (en) | A kind of complete zero Doppler attitude guidance method of Synthetic Aperture Radar satellite | |
CN106643804A (en) | Method of pre-determining footprint position of satellite laser altimeter | |
CN102004250A (en) | Frequency domain expansion based spaceborne/airborne hybrid bistatic synthetic aperture radar imaging method | |
CN105241456A (en) | Loitering munition high-precision combination navigation method | |
CN103344958B (en) | Based on the satellite-borne SAR high-order Doppler parameter evaluation method of almanac data | |
CN103823209B (en) | For low cost kinematic error measurement mechanism in small-sized polarization sensitive synthetic aperture radar system | |
CN103885062A (en) | Double-base foresight SAR moving target imaging method and moving target speed estimation method | |
CN102384755A (en) | High-accuracy method for measuring instant navigational speed of airborne phased array whether radar | |
CN103487808B (en) | A kind of track Simulation method of variable element locking mode missile-borne Spotlight SAR Imaging | |
CN103454633A (en) | Interference SAR movable base line processing method based on back-projection algorithm | |
CN104597446A (en) | Space-borne synthetic aperture radar ground range resolution representation and parameter design method | |
CN103245948B (en) | Image match navigation method for double-area image formation synthetic aperture radars | |
Kauffman et al. | Simulation study of UWB-OFDM SAR for navigation with INS integration | |
Łabowski et al. | Inertial navigation system for radar terrain imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |