CN103278819A - Onboard high-resolution strabismus bunching synthetic aperture radar (SAR) imaging method based on sliding receiving window - Google Patents
Onboard high-resolution strabismus bunching synthetic aperture radar (SAR) imaging method based on sliding receiving window Download PDFInfo
- Publication number
- CN103278819A CN103278819A CN2013101672694A CN201310167269A CN103278819A CN 103278819 A CN103278819 A CN 103278819A CN 2013101672694 A CN2013101672694 A CN 2013101672694A CN 201310167269 A CN201310167269 A CN 201310167269A CN 103278819 A CN103278819 A CN 103278819A
- Authority
- CN
- China
- Prior art keywords
- complex data
- row
- orientation
- dimentional
- data
- 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
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses an onboard high-resolution strabismus bunching synthetic aperture radar (SAR) imaging method based on a sliding receiving window. The onboard high-resolution strabismus bunching SAR imaging method comprises the following steps of 1, reading original echo data and relevant imaging parameters; 2, performing azimuth linear frequency demodulation processing; 3, performing azimuth Fourier transformation processing; 4, performing azimuth linear frequency demodulation residual phase error compensation processing; 5, performing range Fourier transformation processing; 6, performing uniform compression processing; 7, performing stolt interpolation processing; 8, performing azimuth inverse Fourier transformation processing; 9, performing geometric correction processing; and 10, performing range Fourier transformation processing. The invention provides the onboard high-resolution strabismus bunching SAR imaging method based on the sliding receiving window in order to solve the problem that a method for imaging the onboard high-resolution strabismus bunching SAR original echo data based on a sliding receiving window technology does not exist at present.
Description
Technical field
The invention belongs to the signal process field, particularly a kind of Airborne High-resolution stravismus spot beam SAR synthetic-aperture radar (Synthetic Aperture Radar, SAR) formation method based on the slip receiver window.
Background technology
SAR is a kind of active remote sensing device that is operated in microwave frequency band, overcome the defective that optical imagery is limited by weather and illumination condition, can round-the-clock, round-the-clock, carry out remote sensing of the earth observation at a distance, and can penetrate natural vegetation, artificial camouflage etc., improved the information capture ability of radar greatly.Therefore, SAR has become the popular research field of Radar Technology, by more and more countries is paid attention to.Than traditional side-looking SAR system, the Airborne Squint SAR imaging has very high dirigibility and maneuverability in actual applications, and by the adjustment controlling antenna wave beam to point, but the observation area is selected on SAR system freedom and flexibility ground, and can heavily visit the sensitizing range fast, improve the observing capacity of SAR greatly.In addition, along with the raising of SAR resolution, make SAR target reconnaissance and recognition capability be significantly improved.Therefore, in the last few years, the imaging of high resolving power stravismus had become an important developing direction.
But the big Squint SAR imaging of high resolving power has also brought new technological difficulties.On the one hand, along with the increase of stravismus angle, SAR original echo orientation causes the high precision imaging more difficult to more serious to coupling phenomenon with distance; On the other hand, the range migration amount forms geometric growth with the increase of angle of squint and the increase of resolution, causes the remarkable increase of SAR original echo data volume, has strengthened the difficulty of data storage with real time imagery.Slip receiver window technology refers to open constantly by changing the echo window, eliminates once item, i.e. range walk in the range migration.Slip receiver window technology can reduce SAR original echo range migration amount, and then reduces SAR original echo data volume.But, the technological adjustment of slip receiver window the admission zero-time of echo data, in the observation area Doppler's course of each target also with the orientation to variation has taken place, make traditional formation method no longer suitable.
Summary of the invention
To the objective of the invention is in order addressing the above problem, based on slip receiver window technical characterstic, in conjunction with conventional wave number field formation method, to have proposed a kind of Airborne High-resolution stravismus Spotlight SAR Imaging formation method based on the slip receiver window.
A kind of Airborne High-resolution stravismus Spotlight SAR Imaging formation method based on the slip receiver window comprises following step:
Step 1: read in original echo data and dependent imaging parameter;
Read in the Airborne High-resolution stravismus Spotlight SAR Imaging two dimension original echo emulation complex data S based on the slip receiver window
StartAnd corresponding imaging parameters, specifically comprise: the orientation is to sampling number N
a, distance is to sampling number N
r, signal sampling rate f
s, signal bandwidth Bw, pulse width τ, chirp rate b, pulse repetition rate PRF is with reference to oblique distance R
Ref, doppler centroid fd
0, doppler frequency rate f
R0, satellite velocities P
v, equivalent angle of squint
Signal wavelength lambda, signal carrier frequency f
0, signal velocity c;
Step 2: the orientation is handled to separating linear frequency modulation;
With two-dimentional echo simulation complex data S
StartCarrying out the orientation handles to separating linear frequency modulation: complex data S at first
StartWith factor H
1Multiply each other, obtain complex data S
1_1Secondly to complex data S
1_1Do the orientation to Fourier transform, namely carry out fast Fourier along each distance to (by row) and change (FFT), obtain complex data S
1_2At last with complex data S
1_2With factor H
2Multiply each other, obtain final orientation to the complex data S that separates behind the linear frequency modulation
1
Step 3: the orientation is handled to Fourier transform;
The complex data S that step 2 is obtained
1Carry out Fast Fourier Transform (FFT) (FFT) along each distance to (by row), obtain the orientation to frequency spectrum complex data S
2
Step 4: the orientation is handled to separating linear frequency modulation residual phase error compensation;
The complex data S that step 3 is obtained
2With orientation constantly, corresponding orientation to separating linear frequency modulation residual phase error compensation factor Ω
1Multiply each other the complex data S after being compensated
3
Step 5: distance is handled to Fourier transform;
The complex data S that step 4 is obtained
3(by row) carries out Fast Fourier Transform (FFT) (FFT) along each orientation constantly, obtains 2-d spectrum complex data S
4
Step 6: consistently compress processing;
The complex data S that step 5 is obtained
4With corresponding consistent compressibility factor Ω
2Multiply each other, slightly focused on complex data S
5
Step 7: Stolt (stolt) interpolation processing;
The complex data S that obtains for step 6
5, utilize Singh (sinc) method of interpolation to carry out the stolt interpolation processing, obtain being mapped to by two-dimensional frequency the complex data S of two-dimentional wavenumber domain
6
Step 8: the orientation is handled to inverse Fourier transform;
The complex data S that step 7 is obtained
6Carry out inverse fast Fourier transform (IFFT) along each distance to (by row), obtain the orientation time domain apart from wavenumber domain complex data S
7
Step 9: geometry correction is handled;
The complex data S that step 8 is obtained
7With geometry correction factor Ω
3Multiply each other, obtain the complex data S after geometry correction
8
Step 10: distance is handled to Fourier transform;
The complex data S that step 9 is obtained
9(by row) carries out inverse fast Fourier transform (IFFT) along each orientation constantly, obtains final imaging results S
End
The invention has the advantages that:
(1) the present invention proposes a kind of Airborne High-resolution stravismus Spotlight SAR Imaging formation method based on the slip receiver window, solved the present situation that does not have formation method at present based on the Airborne High-resolution stravismus Spotlight SAR Imaging original echo data of slip receiver window technology.
(2) the present invention proposes a kind of Airborne High-resolution stravismus Spotlight SAR Imaging formation method based on the slip receiver window, have the characteristics of high precision focal imaging.Because the formation method that the present invention proposes is a kind of improved wavenumber domain formation method, and the advantage of wavenumber domain formation method is, as long as satisfy this condition of constant airspeed (carried SAR satisfies this condition just), just can realize that high precision focuses on.Therefore, utilize the present invention can realize scene high precision focal imaging.
(3) the present invention proposes a kind of Airborne High-resolution stravismus Spotlight SAR Imaging formation method based on the slip receiver window, have characteristic of strong applicability.On the one hand, because the formation method that the present invention proposes is a kind of improved wavenumber domain formation method, and the wavenumber domain formation method is not looked side ways the restriction of angle, therefore, under the very big condition of stravismus angle, the present invention can realize the vernier focusing of scene equally, obtains high-quality SAR image.On the other hand, the present invention can realize the imaging of 0.1m ultrahigh resolution, and obtains vernier focusing, and therefore, the present invention is applicable to the imaging requirements of present various resolution.
Description of drawings
Fig. 1 is a kind of Airborne High-resolution stravismus Spotlight SAR Imaging formation method process flow diagram based on the slip receiver window that the present invention proposes;
Fig. 2 is embodiment simulating scenes synoptic diagram;
Fig. 3 is the embodiment imaging results;
Fig. 4 is the upper left point target sectional view of embodiment;
Fig. 5 is embodiment intermediate point object profile figure;
Fig. 6 is embodiment lower-right most point object profile figure.
Embodiment
The present invention is described in further detail below in conjunction with drawings and Examples.
The present invention proposes a kind of airborne resolution stravismus Spotlight SAR Imaging formation method based on the slip receiver window, and the object of processing is based on the Airborne High-resolution stravismus Spotlight SAR Imaging original echo data of slip receiver window, and the result who obtains is a panel height resolution oblique view picture.
The slip receiver window refers to SAR when receiving echo, and the echo receiver window is opened constantly and constantly changed with the orientation, thereby reduces the range migration of target, eliminates the range walk (linear segment in the range migration) of target in other words.The SAR system of employing slip receiver window is each orientation echo window unlatching constantly moment T'(i at work) be:
Wherein, T
0Be traditional SAR system's fixed echo window unlatching moment under the same terms, fd
0Refer to doppler centroid, λ refers to signal wavelength, and c refers to signal velocity, and t (i) refers to the orientation constantly, and
, i=0,1,2 ..., N
a-1.Slip receiver window technology has reduced the range migration amount of echo, and then has reduced the echo data amount, has alleviated the pressure of data storage and real time imagery.
The present invention is a kind of Airborne High-resolution stravismus Spotlight SAR Imaging formation method based on the slip receiver window, and idiographic flow may further comprise the steps as shown in Figure 1:
Step 1: read in the Airborne High-resolution stravismus Spotlight SAR Imaging two dimension original echo complex data S based on the slip receiver window
StartAnd corresponding imaging parameters.Wherein, S
StartBe a two-dimentional plural groups, size is N
a* N
r, and the imaging parameter specifically comprises: the orientation is to sampling number N
a, the distance to sampling number N
r, signal sampling rate f
s, signal bandwidth Bw, chirp rate b, pulse repetition rate PRF, with reference to oblique distance R
Ref, doppler centroid fd
0, doppler frequency rate f
R0, satellite velocities P
v, equivalent angle of squint
, signal wavelength lambda, signal carrier frequency f
0, signal velocity c;
Step 2: with two-dimentional original echo complex data S
StartCarry out the orientation and handle to separating linear frequency modulation, specifically can be divided into following step:
(a) two one-dimensional sequence i of structure, j, wherein i represents the orientation to sequence (OK), and j represents distance to sequence (row);
i=[1,2,…,N
a] (2)
j=[1,2,…,N
r]
(b) obtain two-dimentional original echo complex data S
StartThe orientation moment t (i) that each row is corresponding;
(c) with two-dimentional original echo complex data S
StartWith removing twiddle factor H
1Multiply each other, obtain data S
1_1(i, j), factor H wherein
1(i) be that size is N
a* 1 one dimension plural groups, formula is:
H
1(i)=exp{jπ(fr
0·t
2(i)+2·fd
0·t(i))} (4)
Two-dimentional plural groups S then
1_1Can be drawn by following formula:
S
1_1(i,j)=S
start(i,j)·H
1(i) (5)
(d) to complex data S
1_1Carry out fast Fourier along each distance to (by row) and change (FFT), obtain complex data S
1_2
S
1_2(:,j)=FFT(S
1_1(:,j)) (6)
Wherein, S
1_2(:, j) expression S
1_2N row, S
1_1(:, j) expression S
1_1N row, FFT () expression is carried out Fast Fourier Transform (FFT) to one-dimension array.
(e) obtain postrotational equivalent pulse repetition frequency PRF';
(g) convolution (7) obtains the corresponding orientation of the every row of two-dimentional complex data t constantly
1(i);
(h) with complex data S
1_2With factor H
2Multiply each other, obtain final orientation to the data S that separates behind the linear frequency modulation
1Factor H wherein
2(i) be that size is N
a* 1 one dimension plural groups, its formula is:
Two-dimentional plural groups S then
1Can be drawn by following formula:
S
1(i,j)=S
1_2(i,j)·H
2(i) (10)
Step 3: the complex data S that step 2 is obtained
1(i j) carries out Fast Fourier Transform (FFT) (FFT) along each distance to (by row), obtains the orientation to frequency spectrum complex data S
2(i, j);
S
2(:,j)=FFT(S
1(:,j)) (11)
Wherein, S
2(:, j) expression S
2J row, S
1(:, j) expression S
1J row, FFT () expression is carried out Fast Fourier Transform (FFT) to one-dimension array.
Step 4: the complex data S that step 3 is obtained
2(i, j) with orientation constantly, corresponding orientation to separating linear frequency modulation residual phase error compensation factor Ω
1(i) multiply each other the complex data S after being compensated
3(i, j);
(a) convolution (7) obtain two-dimentional orientation to frequency domain distance to time domain complex data S
2(i, j) the corresponding orientation frequency f of every row
a(i);
(b) convolution (12) obtains the big or small N that is
a* 1 one dimension compensating factor Ω
1(i);
(c) obtain two-dimentional complex data S after the compensation
3(i, j);
S
3(i,j)=S
2(i,j)·Ω
1(i) (14)
Step 5: the complex data S that step 4 is obtained
3(i, j) (by row) carries out Fast Fourier Transform (FFT) (FFT) along each orientation constantly, obtains 2-d spectrum complex data S
4(i, j);
S
4(i,:)=FFT(S
3(i,:)) (15)
Wherein, S
3(i :) expression S
3I capable, S
4(i :) expression S
4I capable.
Step 6: the complex data S that step 5 is obtained
4(i is j) with corresponding consistent compressibility factor Ω
2(i j) multiplies each other, and is slightly focused on complex data S
5(i, j).
(a) according to reference oblique distance R
RefObtain the shortest oblique distance R
Min
(b) obtain two-dimensional frequency complex data S
4(i, j) every row correspondence apart from frequency domain f
τ(j);
(c) convolution (12) obtains two-dimensional frequency complex data S with formula (17)
4(i, j) orientation of every row correspondence is to wave number k
x(i) corresponding with every row distance is to wave number k
Rc(j);
(d) convolution (16)~formula (18) is obtained the big or small N that is
a* N
rTwo-dimentional consistent compressibility factor Ω
2(i, j);
(e) convolution (19) obtains the two-dimentional complex data S after consistent the compression
5(i, j);
S
5(i,j)=S
4(i,j)·Ω
2(i,j) (20)
Step 7: the complex data S that obtains for step 6
5(i j), utilizes the sinc method of interpolation to carry out the stolt interpolation processing, obtains being mapped to by two-dimensional frequency the complex data S of two-dimentional wavenumber domain
6(i, j);
(a) according to the distance frequency domain to apart from the mapping relations of wavenumber domain, convolution (18) obtains apart from wavenumber domain wave number k'
Rc(i, j);
(b) traversal is obtained apart from wavenumber domain wave number k'
Rc(i, maximal value k' j)
Rc, maxWith minimum value k'
Rc, max, and obtain apart from branch interval delta k' such as wavenumber domain wave numbers
Rc
(c) obtain two-dimentional wavenumber domain data uniformly apart from the wavenumber domain wave number
(d) obtain two-dimentional wavenumber domain complex data each uniformly apart from the wavenumber domain wave number
At the corresponding inhomogeneous k' of every row
RcPosition p in (i :) (i, j).
Method is for carrying out following operation by row: be example with the position p (1,1) that obtains first row, first row, at first obtain absolute difference
N=[1,2 ..., N
r], obtain minimum absolute difference Δ k
MinWith the n of correspondence position, if
P (1,1)=n-1, if
p(1,1)=n,
By that analogy, obtain each position p (i, j).
(e) the position p that obtained in conjunction with the last step (i, j), obtain the required sampled point of sinc interpolation at interval q (i, j, n);
(f) convolution (24) utilizes the sinc method of interpolation, obtains out the 2-D data S after the stolt interpolation
6, because 2-D data is complex data, need respectively to S
6(i, real part S j)
6_re(i is j) with imaginary part S
6_im(i j) carries out the sinc method of interpolation respectively and obtains and draw.
Wherein, N is interpolation kernel length, and sinc () refers to interpolating function
S
5_re(i j) refers to 2-D data S
5_reThe real part of the capable j row of i, S
5_im(i j) refers to 2-D data S
5_imThe imaginary part of the capable j row of i.
Step 8: the complex data S that step 7 is obtained
6(i j) carries out inverse fast Fourier transform (IFFT) along each distance to (by row), obtains the orientation time domain apart from wavenumber domain complex data S
7(i, j);
S
7(:,j)=IFFT(S
6(:,j)) (27)
Wherein, S
6(:, j) expression S
6J row, S
7(:, j) expression S
7J row, IFFT () expression is carried out inverse fast Fourier transform to one-dimension array.
Step 9: the complex data S that step 8 is obtained
7(i is j) with geometry correction factor Ω
3(i j) multiplies each other, and obtains the complex data S after geometry correction
8(i, j);
(a) convolution (8) and formula (12) obtain geometry correction factor Ω
4(i, j);
(b) utilize formula (30) to obtain complex data S after geometry correction
8(i, j);
S
8(i,j)=S
7(i,j)·Ω
3(i,j) (29)
Step 10: the complex data S that step 9 is obtained
8(i, j) (by row) carries out inverse fast Fourier transform (IFFT) along each orientation constantly, obtains final imaging results S
End(i, j);
S
end(i,:)=IFFT(S
8(i,:)) (30)
Wherein, S
8(i :) expression S
8I capable, S
End(i :) expression S
EndI capable.
Embodiment:
Present embodiment proposes a kind of Airborne High-resolution stravismus Spotlight SAR Imaging formation method based on the slip receiver window, simulating scenes is 3 * 3 dot matrix, its point and the 100m that is spaced apart that puts, finally respectively upper left in the simulating scenes, middle, three points in bottom right are assessed, concrete simulating scenes as shown in Figure 2, the imaging parameters that relates in its imaging process is as shown in table 1.
Table 1 embodiment parameter
Present embodiment specifically may further comprise the steps:
Step 1: read in the Airborne High-resolution stravismus Spotlight SAR Imaging two dimension original echo complex data S based on the slip receiver window
StartAnd corresponding imaging parameters.Wherein, S
StartBe two-dimentional plural groups, size is 65536 * 16384, and concrete imaging parameters is as shown in table 1;
Step 2: with two-dimentional original echo complex data S
StartCarry out the orientation and handle to separating linear frequency modulation, the concrete operations step is:
(a) structuring one-dimensional sequence, as the formula (2), i=[1,2 ..., 65536], j=[1,2 ..., 16384];
(b) obtain two-dimentional original echo complex data S
StartThe orientation moment t (i) that every row is corresponding, detailed process is undertaken by formula (3);
(c) with data S
StartWith removing twiddle factor H
1(i) multiply each other, obtain data S
1_1(i, j).Factor H wherein
1Be that size is 65536 * 1 one dimension plural groups, concrete acquisition process is undertaken by formula (4), and two-dimentional plural groups S
1_1(i, j) acquisition process is undertaken by formula (5);
(d) to complex data S
1_1(i j) changes (FFT) along each distance to carrying out fast Fourier by formula (6), obtains complex data S
1_2(i, j);
(e) obtain postrotational equivalent pulse repetition frequency PRF', acquisition process is undertaken by formula (7);
(g) convolution (7) obtains the corresponding orientation of the every row of two-dimentional complex data t constantly
1(i), acquisition process is undertaken by formula (8);
(h) with data S
1_2(i is j) with factor H
2(i) multiply each other, obtain final orientation to the data S that separates behind the linear frequency modulation
1(i, j).Factor H wherein
2(i) be that size is 65536 * 1 one dimension plural groups, its acquisition process is undertaken by formula (9), and two-dimentional plural groups S
1Acquisition process is undertaken by formula (10);
Step 3: the complex data S that step 2 is obtained
1(i j) carries out Fast Fourier Transform (FFT) (FFT) along each distance to (by row), obtains the orientation to frequency spectrum complex data S
2(i, j), specific operation process is undertaken by formula (11);
Step 4: the complex data S that step 3 is obtained
2(i, j) with orientation constantly, corresponding orientation to separating linear frequency modulation residual phase error compensation factor Ω
1(i) multiply each other the complex data S after being compensated
3(i, j);
(a) obtain two-dimentional orientation to frequency domain distance to time domain complex data S
2(i, j) the corresponding orientation frequency f of every row
a, concrete acquisition process is undertaken by formula (12);
(b) to obtain size be 65536 * 1 one dimension compensating factor Ω to convolution (7)
1(i), concrete acquisition process is undertaken by formula (13);
(c) obtain two-dimentional complex data S after the compensation
3(i, j), concrete acquisition process is undertaken by formula (14);
Step 5: the complex data S that step 4 is obtained
3(i, j) (by row) carries out Fast Fourier Transform (FFT) (FFT) along each orientation constantly, obtains 2-d spectrum complex data S
4(i, j), specific operation process is undertaken by formula (15);
Step 6: the complex data S that step 5 is obtained
4(i is j) with corresponding consistent compressibility factor Ω
2(i j) multiplies each other, and is slightly focused on complex data S
5(i, j).
(a) according to reference oblique distance R
Ref=19.31km obtains the shortest oblique distance R
Min, concrete acquisition process is undertaken by formula (16);
(b) obtain two-dimensional frequency complex data S
4(i, j) every row correspondence apart from frequency domain f
τ(j), concrete acquisition process is undertaken by formula (17);
(c) convolution (7) obtains two-dimensional frequency complex data S with formula (11)
4(i, j) orientation of every row correspondence is to wave number k
x(i) corresponding with every row distance is to wave number k
Rc(j), concrete acquisition process is undertaken by formula (18);
(d) to obtain size be 65536 * 16384 two-dimentional consistent compressibility factor Ω to convolution (10)~formula (12)
2(i, j), concrete acquisition process is undertaken by formula (19);
(e) convolution (13) obtains the two-dimentional complex data S after consistent the compression
5(i, j), concrete acquisition process is undertaken by formula (20);
Step 7: the complex data S that obtains for step 6
5(i j), utilizes the sinc method of interpolation to carry out the stolt interpolation processing, obtains being mapped to by two-dimensional frequency the complex data S of two-dimentional wavenumber domain
6(i, j);
(a) according to the distance frequency domain to apart from the mapping relations of wavenumber domain, convolution (18) obtains apart from wavenumber domain wave number k'
Rc(i, j), concrete acquisition process is undertaken by formula (21);
(b) traversal is obtained apart from wavenumber domain wave number k'
Rc(i, maximal value k' j)
Rc, max=69.12rad/s and minimum value k'
Rc, max=-105.06rad/s, and obtain apart from branch interval delta k' such as wavenumber domain wave numbers
Rc, concrete operations are undertaken by formula (22);
(c) obtain two-dimentional wavenumber domain complex data uniformly apart from the wavenumber domain wave number
Concrete operations are undertaken by formula (23);
(d) obtain two-dimentional wavenumber domain complex data each uniformly apart from the wavenumber domain wave number
At the corresponding inhomogeneous k' of every row
RcPosition p in (i :) (i, j).Method is for carrying out following operation by row: obtain the position p (1,1) of first row, first row earlier, obtain absolute difference
N=[1,2 ..., 16384], obtain minimum absolute difference Δ k
MinWith the n of correspondence position, if
P (1,1)=n-1, if
p(1,1)=n,
By that analogy, obtain each position p (i, j).
(e) the position p that obtained in conjunction with the last step (i, j), obtain the required sampled point of sinc interpolation at interval q (n), concrete operations are undertaken by formula (24) for i, j;
(f) convolution (18) utilizes the sinc method of interpolation, and selecting sinc interpolation kernel length is N=8, obtains out the two-dimentional complex data S after the stolt interpolation
6(i j), because 2-D data is complex data, needs respectively to S
6(i, real part S j)
6_re(i is j) with imaginary part S
6_im(i j) carries out the sinc method of interpolation respectively and obtains and draw, and concrete operations are undertaken by formula (25) (26).
Step 8: the complex data S that step 7 is obtained
6(i j) carries out inverse fast Fourier transform (IFFT) along each distance to (by row), obtains the orientation time domain apart from wavenumber domain complex data S
7(i, j), concrete operations are undertaken by formula (27);
Step 9: the complex data S that step 8 is obtained
7(i is j) with geometry correction factor Ω
3(i j) multiplies each other, and obtains the complex data S after geometry correction
8(i, j);
(a) convolution (8) and formula (12) obtain geometry correction factor Ω
3(i, j), concrete operations are undertaken by formula (28);
(b) utilize formula (28) to obtain complex data S after geometry correction
8(i, j), concrete operations are undertaken by formula (29);
Step 10: the complex data S that step 9 is obtained
8(i, j) (by row) carries out inverse fast Fourier transform (IFFT) along each orientation constantly, obtains final imaging results S
End(i, j), concrete operations are undertaken by formula (30);
Imaging processing through above-mentioned steps obtains final scene imaging results as shown in Figure 3.Table 2 has provided that scene is upper left, middle, the Imaging Evaluation result of three point targets in bottom right, and Fig. 4, Fig. 5, Fig. 6 have provided respectively that scene is upper left, middle, the sectional view of three point targets in bottom right.
Table two Imaging Evaluation result
According to table 2 assessment result and Fig. 4~sectional view shown in Figure 6, can draw: on the one hand, this formation method is 70 still can vernier focusing when spending in the stravismus angle, illustrate that the method for the present invention's proposition do not looked side ways the restriction of angle; On the other hand, this formation method still can vernier focusing for the 0.1m super-resolution, illustrates that the method that the present invention proposes can realize vernier focusings to present various resolution.Therefore, method proposed by the invention can realize the Airborne High-resolution angle of squint Spotlight SAR Imaging accurately image based on the slip receiver window, has obtained high-precision imaging results.
Claims (1)
1. the Airborne High-resolution based on the slip receiver window is looked side ways the Spotlight SAR Imaging formation method, may further comprise the steps:
Step 1: read in the Airborne High-resolution stravismus Spotlight SAR Imaging two dimension original echo complex data S based on the slip receiver window
StartAnd corresponding imaging parameters;
S
StartBe a two-dimentional plural groups, size is N
a* N
r, imaging parameters comprises: the orientation is to sampling number N
a, the distance to sampling number N
r, signal sampling rate f
s, signal bandwidth Bw, chirp rate b, pulse repetition rate PRF, with reference to oblique distance R
Ref, doppler centroid fd
0, doppler frequency rate f
R0, satellite velocities P
v, equivalent angle of squint
Signal wavelength lambda, signal carrier frequency f
0, signal velocity c;
Step 2: with two-dimentional original echo complex data S
StartCarry out the orientation and handle to separating linear frequency modulation, specifically comprise following step:
(a) two one-dimensional sequence i of structure, j, wherein i represents the orientation to sequence, and j represents distance to sequence;
i=[1,2,…,N
a] (2)
j=[1,2,…,N
r]
(b) obtain two-dimentional original echo complex data S
StartThe orientation moment t (i) that each row is corresponding;
(c) with two-dimentional original echo complex data S
StartWith remove twiddle factor H
1Multiply each other, obtain data S
1_1(i, j), factor H wherein
1(i) be that size is N
a* 1 one dimension plural groups, formula is:
H
1(i)=exp{jπ(fr
0·t
2(i)+2·fd
0·t(i))} (4)
Two-dimentional plural groups S then
1_1Drawn by following formula:
S
1_1(i,j)=S
start(i,j)·H
1(i) (5)
(d) to complex data S
1_1Change to carrying out fast Fourier along each distance, obtain complex data S
1_2
S
1_2(:,j)=FFT(S
1_1(:,j)) (6)
Wherein, S
1_2(:, j) expression S
1_2N row, S
1_1(:, j) expression S
1_1N row, FFT () expression is carried out Fast Fourier Transform (FFT) to one-dimension array;
(e) obtain postrotational equivalent pulse repetition frequency PRF';
(g) convolution (7) obtains the corresponding orientation of the every row of two-dimentional complex data t constantly
1(i);
(h) with complex data S
1_2With factor H
2Multiply each other, obtain final orientation to the data S that separates behind the linear frequency modulation
1Factor H wherein
2(i) be that size is N
a* 1 one dimension plural groups, its formula is:
Two-dimentional plural groups S then
1Drawn by following formula:
S
1(i,j)=S
1_2(i,j)·H
2(i) (10)
Step 3: the complex data S that step 2 is obtained
1(i, j) along each the distance to carrying out Fast Fourier Transform (FFT), obtain the orientation to frequency spectrum complex data S
2(i, j);
S
2(:,j)=FFT(S
1(:,j)) (11)
Wherein, S
2(:, j) expression S
2J row, S
1(:, j) expression S
1J row, FFT () expression is carried out Fast Fourier Transform (FFT) to one-dimension array;
Step 4: the complex data S that step 3 is obtained
2(i, j) with orientation constantly, corresponding orientation to separating linear frequency modulation residual phase error compensation factor Ω
1(i) multiply each other the complex data S after being compensated
3(i, j);
(a) convolution (7) obtain two-dimentional orientation to frequency domain distance to time domain complex data S
2(i, j) the corresponding orientation frequency f of every row
a(i);
(b) convolution (12) obtains the big or small N that is
a* 1 one dimension compensating factor Ω
1(i);
(c) obtain two-dimentional complex data S after the compensation
3(i, j);
S
3(i,j)=S
2(i,j)·Ω
1(i) (14)
Step 5: the complex data S that step 4 is obtained
3(i j) carries out Fast Fourier Transform (FFT) constantly along each orientation, obtain 2-d spectrum complex data S
4(i, j);
S
4(i,:)=FFT(S
3(i,:)) (15)
Wherein, S
3(i :) expression S
3I capable, S
4(i :) expression S
4I capable;
Step 6: the complex data S that step 5 is obtained
4(i is j) with corresponding consistent compressibility factor Ω
2(i j) multiplies each other, and is slightly focused on complex data S
5(i, j);
(a) according to reference oblique distance R
RefObtain the shortest oblique distance R
Min
(b) obtain two-dimensional frequency complex data S
4(i, j) every row correspondence apart from frequency domain f
τ(j);
(c) convolution (12) obtains two-dimensional frequency complex data S with formula (17)
4(i, j) orientation of every row correspondence is to wave number k
x(i) corresponding with every row distance is to wave number k
Rc(j);
(d) convolution (16)~formula (18) is obtained the big or small N that is
a* N
rTwo-dimentional consistent compressibility factor Ω
2(i, j);
(e) convolution (19) obtains the two-dimentional complex data S after consistent the compression
5(i, j);
S
5(i,j)=S
4(i,j)·Ω
2(i,j) (20)
Step 7: the complex data S that obtains for step 6
5(i j), utilizes the sinc method of interpolation to carry out the stolt interpolation processing, obtains being mapped to by two-dimensional frequency the complex data S of two-dimentional wavenumber domain
6(i, j);
(a) according to the distance frequency domain to apart from the mapping relations of wavenumber domain, convolution (18) obtains apart from wavenumber domain wave number k'
Rc(i, j);
(b) traversal is obtained apart from wavenumber domain wave number k'
Rc(i, maximal value k' j)
Rc, maxWith minimum value k'
Rc, max, and obtain apart from branch interval delta k' such as wavenumber domain wave numbers
Rc
(c) obtain two-dimentional wavenumber domain data uniformly apart from the wavenumber domain wave number
(d) obtain two-dimentional wavenumber domain complex data each uniformly apart from the wavenumber domain wave number
At the corresponding inhomogeneous k' of every row
RcPosition p in (i :) (i, j);
Be specially: be example with the position p (1,1) that obtains first row first row, at first obtain absolute difference
N=[1,2 ..., N
r], obtain minimum absolute difference Δ k
MinWith the n of correspondence position, if
P (1,1)=n-1, if
P (1,1)=n, in like manner, by that analogy, obtain each position p (i, j);
(e) the position p that obtained in conjunction with the last step (i, j), obtain the required sampled point of sinc interpolation at interval q (i, j, n);
(f) convolution (24) utilizes the sinc method of interpolation, obtains out the 2-D data S after the stolt interpolation
6, because 2-D data is complex data, need respectively to S
6(i, real part S j)
6_re(i is j) with imaginary part S
6_im(i j) carries out the sinc method of interpolation respectively and obtains and draw;
Wherein, N is interpolation kernel length, and sinc () refers to interpolating function
S
5_re(i j) refers to 2-D data S
5_reThe real part of the capable j row of i, S
5_im(i j) refers to 2-D data S
5_imThe imaginary part of the capable j row of i;
Step 8: the complex data S that step 7 is obtained
6(i, j) along each the distance to carrying out inverse fast Fourier transform, obtain the orientation time domain apart from wavenumber domain complex data S
7(i, j);
S
7(:,j)=IFFT(S
6(:,j)) (27)
Wherein, S
6(:, j) expression S
6J row, S
7(:, j) expression S
7J row, IFFT () expression is carried out inverse fast Fourier transform to one-dimension array;
Step 9: the complex data S that step 8 is obtained
7(i is j) with geometry correction factor Ω
3(i j) multiplies each other, and obtains the complex data S after geometry correction
8(i, j);
(a) convolution (8) and formula (12) obtain geometry correction factor Ω
4(i, j);
(b) utilize formula (30) to obtain complex data S after geometry correction
8(i, j);
S
8(i,j)=S
7(i,j)·Ω
3(i,j) (29)
Step 10: the complex data S that step 9 is obtained
8(i j) carries out inverse fast Fourier transform constantly along each orientation, obtain final imaging results S
End(i, j);
S
end(i,:)=IFFT(S
8(i,:)) (30)
Wherein, S
8(i :) expression S
8I capable, S
End(i :) expression S
EndI capable.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310167269.4A CN103278819B (en) | 2013-05-08 | 2013-05-08 | Onboard high-resolution strabismus bunching synthetic aperture radar (SAR) imaging method based on sliding receiving window |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310167269.4A CN103278819B (en) | 2013-05-08 | 2013-05-08 | Onboard high-resolution strabismus bunching synthetic aperture radar (SAR) imaging method based on sliding receiving window |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103278819A true CN103278819A (en) | 2013-09-04 |
CN103278819B CN103278819B (en) | 2015-04-15 |
Family
ID=49061392
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310167269.4A Active CN103278819B (en) | 2013-05-08 | 2013-05-08 | Onboard high-resolution strabismus bunching synthetic aperture radar (SAR) imaging method based on sliding receiving window |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103278819B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103728619A (en) * | 2014-01-16 | 2014-04-16 | 北京航空航天大学 | Variable repetition frequency based airborne high-squint stripe Synthetic Aperture Radar (SAR) imaging method |
CN104898119A (en) * | 2015-05-19 | 2015-09-09 | 北京航空航天大学 | Correlation function-based moving-target parameter estimation method |
CN105607057A (en) * | 2014-11-21 | 2016-05-25 | 中国航空工业集团公司雷华电子技术研究所 | Airborne SAR real echo data transformation method |
WO2022155921A1 (en) * | 2021-01-22 | 2022-07-28 | 京东方科技集团股份有限公司 | Display panel and driving method therefor, compensation data compression method, and compensation data decompression method |
CN116148856A (en) * | 2023-04-17 | 2023-05-23 | 南京邮电大学 | SAR moving target two-dimensional self-focusing imaging processing method |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102176016A (en) * | 2011-01-25 | 2011-09-07 | 北京航空航天大学 | Large squint sliding spotlight SAR (synthetic aperture radar) imaging processing method |
-
2013
- 2013-05-08 CN CN201310167269.4A patent/CN103278819B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102176016A (en) * | 2011-01-25 | 2011-09-07 | 北京航空航天大学 | Large squint sliding spotlight SAR (synthetic aperture radar) imaging processing method |
Non-Patent Citations (2)
Title |
---|
杨威 等: "星载SAR滑动聚束模式三步扩展算法", 《北京航空航天大学学报》 * |
王鹏波 等: "基于改进Deramp处理的滑动聚束SAR成像处理算法", 《仪器仪表学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103728619A (en) * | 2014-01-16 | 2014-04-16 | 北京航空航天大学 | Variable repetition frequency based airborne high-squint stripe Synthetic Aperture Radar (SAR) imaging method |
CN103728619B (en) * | 2014-01-16 | 2016-02-17 | 北京航空航天大学 | Based on the airborne large stravismus stripmap SAR formation method becoming repetition technology |
CN105607057A (en) * | 2014-11-21 | 2016-05-25 | 中国航空工业集团公司雷华电子技术研究所 | Airborne SAR real echo data transformation method |
CN104898119A (en) * | 2015-05-19 | 2015-09-09 | 北京航空航天大学 | Correlation function-based moving-target parameter estimation method |
CN104898119B (en) * | 2015-05-19 | 2018-08-24 | 北京航空航天大学 | A kind of moving target parameter estimation method based on correlation function |
WO2022155921A1 (en) * | 2021-01-22 | 2022-07-28 | 京东方科技集团股份有限公司 | Display panel and driving method therefor, compensation data compression method, and compensation data decompression method |
CN116148856A (en) * | 2023-04-17 | 2023-05-23 | 南京邮电大学 | SAR moving target two-dimensional self-focusing imaging processing method |
Also Published As
Publication number | Publication date |
---|---|
CN103278819B (en) | 2015-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | A frequency-domain imaging algorithm for highly squinted SAR mounted on maneuvering platforms with nonlinear trajectory | |
EP3144702B1 (en) | Method and device for synthethic aperture radar imaging based on non-linear frequency modulation signal | |
EP2650695B1 (en) | Imaging method for synthetic aperture radar in high squint mode | |
An et al. | Performance evaluation of frequency-domain algorithms for chirped low frequency UWB SAR data processing | |
CN110632594B (en) | Long-wavelength spaceborne SAR imaging method | |
CN104597447B (en) | A kind of big stravismus of sub-aperture SAR improves Omega K imaging method | |
WO1993002368A1 (en) | Synthetic aperture radar processor to handle large squint with high phase and geometric accuracy | |
Wu et al. | An azimuth frequency non-linear chirp scaling (FNCS) algorithm for TOPS SAR imaging with high squint angle | |
CN103278819B (en) | Onboard high-resolution strabismus bunching synthetic aperture radar (SAR) imaging method based on sliding receiving window | |
CN104931966A (en) | DCS algorithm-based satellite-borne video SAR (synthetic aperture radar) imaging processing method | |
CN103728619B (en) | Based on the airborne large stravismus stripmap SAR formation method becoming repetition technology | |
CN105093224A (en) | High squint synthetic aperture radar imaging processing method | |
CN102288961A (en) | Imaging method for synthetic aperture radar nonlinear frequency modulation label change | |
CN102226841A (en) | Synchronous orbit SAR imaging method based on high-order polynomial range equation | |
CN104898119A (en) | Correlation function-based moving-target parameter estimation method | |
CN102819020A (en) | Synthetic aperture radar imaging method for azimuth-direction nonlinear chirp scaling of diving model | |
Yang et al. | A novel three-step image formation scheme for unified focusing on spaceborne SAR data | |
Qiu et al. | Bistatic SAR data processing algorithms | |
CN104062657A (en) | Generalized polar coordinate imaging method for synthetic aperture radar (SAR) | |
CN102788978A (en) | Squint spaceborne/airborne hybrid bistatic synthetic aperture radar imaging method | |
CN114325704B (en) | Rapid time domain imaging method of synthetic aperture radar based on wave number spectrum stitching | |
Chen et al. | Iterative minimum entropy algorithm for refocusing of moving targets in SAR images | |
CN105676223B (en) | A kind of non-at the uniform velocity Terahertz radar platform SAR imaging methods | |
CN109143236B (en) | Bistatic bunching SAR large-scene imaging method suitable for complex flight trajectory | |
CN105759264A (en) | Micro-motion target defect echo high-resolution imaging method based on time-frequency dictionary |
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 |