CN103954964B - The method that multi-angle data of synthetic aperture radar obtains - Google Patents

Abstract

The present invention provides a kind of method that multi-angle data of synthetic aperture radar obtains. the method utilize band synthetic-aperture radar mode to enroll advantage that atural object scene multi-angle scatter echo data can break through the limited constraint of scene, in Texas tower moving process, by the multiple transmit/receive antennas adopting different azimuth to point to, form multiple antenna irradiation wave beam, launch and receive the pulse signal with orthogonal encoding characteristics, multiple visual angles scatter echo data of Same Scene or target are enrolled with striped patterns, and obtained by image processing method comprise multi-angle scattered information on a large scale, continuous print target stripe scene two dimension radar image, be applicable to scattering properties with azimuthal variation target scene as buildings carry out to observe.

Description

The method that multi-angle data of synthetic aperture radar obtains

Technical field

The present invention relates to Synthetic Aperture Radar Technique field, particularly relate to a kind of method that multi-angle data of synthetic aperture radar obtains.

Background technology

Stripmap synthetic aperture radar data acquiring mode adopts single antenna wave beam to irradiate target scene with positive side-looking or stravismus mode usually, enrolls the echo data of scene scattering in point with striped patterns. Owing to major part target especially man-made target such as buildings exists more complicated limit, angle structure, thus there is the scattering properties with orientation visual angle change, and traditional band synthetic-aperture radar wave beam orientation is pointed to fixing, the echo data comprising target fixing visual angle, orientation scattering properties can only be obtained.

In order to obtain the scattering properties of target scene multi-angle, just circle mark synthetic-aperture radar and wide angle synthetic-aperture radar obtain target multi-angle scattering properties data and expand some research (Lin both at home and abroad, Yun, etal. " AirbornecircularSARimaging:ResultsatP-band. " 2012IEEEInternationalGeoscienceandRemoteSensingSymposium (IGARSS), 2012; Qi, Liu, etal. " StudyonSARimageformationforaspect-dependentscatterers. " 2ndAsian-PacificConferenceonSyntheticApertureRadar, 2009.APSAR2009.IEEE, 2009.), but circle mark synthetic-aperture radar and wide angle synthetic-aperture radar are all obtain echo data for a certain range limited atural object scene with circle/curved path and pack mode, have certain limitation.

Summary of the invention

(1) technical problem solved

In view of above-mentioned technical problem, the present invention provides a kind of method that multi-angle data of synthetic aperture radar obtains, with realize to scattering properties with azimuthal variation atural object scene as buildings carry out to observe, obtain on a large scale, continuous print target stripe scene two dimension radar image.

(2) technical scheme

The method that multi-angle data of synthetic aperture radar of the present invention obtains comprises: steps A: along orientation to X on Texas tower, arranges the transmit/receive antenna T that row are total to N number of transmit-receive sharing mode_{X, i}/R_{X, i}, make Texas tower with speed V along orientation to carrying out linear uniform motion, wherein, T_{X, i}Represent that i-th transmit-receive sharing Antenna Operation is in transmitting antenna state, R_{X, i}Represent that i-th transmit-receive sharing Antenna Operation is in receiving antenna state, i=1,2 ..., N; Step B: Texas tower often moves V/PRF interval along orientation to X, N number of transmit/receive antenna T_{X, i}/R_{X, i}Just carrying out primary emission receiving course, the scatter echo data that N number of receiving antenna receives respectively form raw radar data collection R after treatment:

R={r₁(��, t; ��₁, �� ��) ..., r_i(��, t; ��_i, �� ��) ... r_N(��, t; ��_N, �� ��) }

Wherein, r_i(��, t; ��_i, �� ��) and it is i-th receiving antenna R_{X, i}The scatter echo data received are via the data after down coversion and AD sampling, and PRF is the pulse-repetition frequency of radar work; Step C: to each receiving antenna R corresponding in raw radar data collection R_{X, i}Two-dimentional echo data part r_i(��, t; ��_i, �� ��), isolate only by transmitting antenna T according to the mode demodulation of orthogonal coding_{X, i}Launch, received antenna R after target scattering_{X, i}The part receivedForm new echo data collectionWherein $\hat{R}=\{{\hat{r}}_1(\mathrm{\η},t;{\mathrm{\θ}}_1,\mathrm{\Δ\θ}),...{\hat{r}}_i(\mathrm{\η},t;{\mathrm{\θ}}_i,\mathrm{\Δ\θ}),...{\hat{r}}_N(\mathrm{\η},t;{\mathrm{\θ}}_N,\mathrm{\Δ\θ})\};$ Step D: to new echo data collectionIn 2-D data after each separationUse two-dimensional imaging algorithm to carry out imaging processing, obtain the band two dimension radar image g of N number of same target area_i(x, y; ��_i, �� ��), form image collection G:

G={g₁(x, y; ��_i, �� ��) ... g_i(x, y; ��_i, �� ��) ... g_N(x, y; ��_i, �� ��) }

Step e: to the band two dimension radar image g of the N number of same target area in image collection G_i(x, y; ��_i, �� ��) and carry out image registration, obtain images after registration collection

$\hat{G}=\{{\hat{g}}_1(x,y;{\mathrm{\θ}}_i,\mathrm{\Δ\θ}),...{\hat{g}}_i(x,y;{\mathrm{\θ}}_i,\mathrm{\Δ\θ}),...{\hat{g}}_N(x,y;{\mathrm{\θ}}_i,\mathrm{\Δ\θ})\};$ And step F, it may also be useful to multi-angle information fusion and extracting method image collection after registrationMiddle extraction target stripe scene range of scatter angles [��_N-�� ��/2, ��₁+ �� ��/2] multi-angle scattered information, obtain target stripe scene two dimension radar image M (x, y) by this multi-angle scattered information.

(3) useful effect

From technique scheme it may be seen that the method that multi-angle data of synthetic aperture radar of the present invention obtains has following useful effect:

(1) by adopting multiple transmit/receive antennas of different azimuth sensing, form multiple antenna irradiation wave beam, overcome tradition band synthetic-aperture radar wave beam orientation and point to the single limitation of scattering data angle information fixing, that obtain, can obtaining the scattered information of target scene scatter point multi-angle, scattered information is abundanter.

(2) it is undertaken by striped patterns to the acquisition of multi-angle scattering data, the situation range limited compared to other data of multiple angles obtain manner atural object scenes such as circle mark synthetic-aperture radar and wide angle synthetic-aperture radar, it is possible to obtain comprise multi-angle scattered information on a large scale, continuous print target stripe scene two dimension radar image.

Accompanying drawing explanation

Fig. 1 is the schema of multi-angle data of synthetic aperture radar acquisition methods according to embodiments of the present invention;

Fig. 2 is the schematic diagram that multi-angle data of synthetic aperture radar obtains;

Fig. 3 is the schematic diagram of radar scanner and wave beam;

Fig. 4 is that building target is irradiated schematic diagram by different visual angles wave beam.

Embodiment

For making the object, technical solutions and advantages of the present invention clearly understand, below in conjunction with specific embodiment, and with reference to accompanying drawing, the present invention is described in more detail. It should be noted that, in accompanying drawing or specification sheets describe, similar or identical part all uses identical figure number. The implementation not illustrated in accompanying drawing or describe is form known to a person of ordinary skill in the art in art. In addition, although herein can providing package containing the demonstration of the parameter of particular value, it should be understood that parameter is without the need to definitely equaling corresponding value, but can be similar to corresponding value in acceptable error tolerance limit or design constraint. The direction term mentioned in embodiment, such as " on ", D score, "front", "rear", "left", "right" etc., be only the direction with reference to accompanying drawing. Therefore, it may also be useful to direction term be used to illustrate not be used for limiting the scope of the invention.

The present invention's utilize band synthetic-aperture radar mode to enroll advantage that atural object scene multi-angle scatter echo data can break through the limited constraint of scene, in Texas tower moving process, by the multiple transmit/receive antennas adopting different azimuth to point to, form multiple antenna irradiation wave beam, launch and receive the pulse signal with orthogonal encoding characteristics, multiple visual angles scatter echo data of Same Scene or target are enrolled with striped patterns, and obtained by image processing method comprise multi-angle scattered information on a large scale, continuous print target stripe scene two dimension radar image, be applicable to scattering properties with azimuthal variation target scene as buildings carry out to observe.

In one exemplary embodiment of the present invention, it provides a kind of method that multi-angle data of synthetic aperture radar obtains. Fig. 1 is the schema of multi-angle data of synthetic aperture radar acquisition methods according to embodiments of the present invention, wherein X, Y and Z represent imageable target region three-dimensional space rectangular coordinates axle, and wherein X is along the flight path direction of Texas tower, for orientation to, Y be distance to, Z be elevation to; X, y and z represent the coordinate figure along coordinate axis respectively. As shown in Figure 1, the present embodiment multi-angle data of synthetic aperture radar acquisition methods comprises:

Steps A: along orientation to X on Texas tower, arranges the transmit/receive antenna T that row are total to N number of transmit-receive sharing mode_{X, i}/R_{X, i}, wherein, i=1,2 ..., N, T_{X, i}Represent that i-th transmit-receive sharing Antenna Operation is in transmitting antenna state, R_{X, i}Represent that i-th transmit-receive sharing Antenna Operation is in receiving antenna state, make Texas tower with speed V along orientation to carrying out linear uniform motion;

In the present embodiment, the form of antenna can be single loudspeaker antenna, it is also possible to be that microstrip antenna forms antenna submatrix;

Step B: Texas tower often moves V/PRF interval along orientation to X, N number of transmit/receive antenna T_{X, i}/R_{X, i}Just carrying out primary emission receiving course, form raw radar data collection R, wherein, this original two dimensional echo data R comprises same target by multiple wave beam with the scattered information after different visual angles irradiation, and PRF is the pulse-repetition frequency of radar work;

Fig. 2 is the schematic diagram that multi-angle data of synthetic aperture radar obtains, for five different sensing wave beams in figure. Fig. 3 is the schematic diagram of radar scanner and wave beam. Please refer to Fig. 2 and Fig. 3, in this step B, for launch and accept process once, comprising:

Sub-step B1, in slow time-sampling moment ��, N number of transmitting antenna forms different azimuth and points to ��₁... ��_i... ��_N, N number of beam of an antenna of bandwidth �� ��;

Wherein ��_iBeing the sensing angle, wave beam center of i-th antenna irradiation wave beam, �� �� is its bandwidth, meets ��_i=��_i-1-�� ��, i=1,2 ..., N, ��_i> 0 represents forward sight, ��_i< 0 looks after representing, so the total visual angle range Theta �� of formation [-| ��_N|-�� ��/2, ��₁+ �� ��/2], ��₁And ��_NThe orientation being respectively the 1st transmit/receive antenna and N number of transmit/receive antenna is pointed to;

Sub-step B2, the pulse signal s with orthogonal encoding characteristics launched respectively by N number of transmitting antenna₁(t) ... s_i(t) ... s_N(t);

Wherein, s_iT () represents the pulse signal of i-th transmission antennas transmit, t represents the fast time coordinate of echo data, pulse signal s₁(t) ... s_i(t) ... s_NT () is launched by corresponding transmitting antenna simultaneously, and have orthogonal encoding characteristics, that is:

$\underset{0}{\overset{T_P}{\&Integral;}}{s}_i\left(t\right)s_j^{*}(t+\mathrm{\&tau;})\mathrm{dt}={\mathrm{\&delta;}}_{\mathrm{ij}},$ When ��=0 (1) wherein, T_pFor transponder pulse time width, �� is time delay variable, ��_ijFor Kronecker punching swashs function, it is defined as:

Sub-step B3, N number of receiving antenna R_{X, 1}... R_{X, i}... R_{X, N}Receive the scatter echo data of the described pulse signal of Same Scene or target simultaneously;

Sub-step B4, the scatter echo data received by each receiving antenna carry out down coversion and AD sampling, and result are saved in raw radar data collection R and go:

R={r₁(��, t; ��₁, �� ��) ..., r_i(��, t; ��_i, �� ��) ... r_N(��, t; ��_N, �� ��) } (3)

Wherein, r_i(��, t; ��_i, �� ��) and it is the data after the echo data that i-th receiving antenna receives is sampled via down coversion and AD, its expression formula is:

$r_i(\mathrm{\&eta;},t;{\mathrm{\&theta;}}_i,\mathrm{\&Delta;\&theta;})=\underset{{\mathrm{\&Omega;}}_i}{\&Integral;}{\mathrm{\&sigma;}}_P({\mathrm{\&theta;}}_i,\mathrm{\&Delta;\&theta;})w_a(\mathrm{\&eta;}-{\mathrm{\&eta;}}_0+{\mathrm{\&eta;}}_c)s_i(t-\frac{2R_P^i}{c})\exp (-j\frac{4\mathrm{\&pi;}{R}_P^i}{\mathrm{\&lambda;}})\mathrm{dxdydz}---\left(4\right)$

Wherein, ��_iRepresent i-th wave beam wave cover district, in it, the scattering properties of arbitrary scattering target P is with orientation visual angle change, but in the scope of �� �� approximate constant, ��_P(��_i, �� ��) and it is scattering coefficient when target P is irradiated by i-th wave beam, c is electromagnetic wave propagation speed, and �� is radar work mid-frequency, w_a(��) represent that orientation is to antenna pattern, ��₀Represent that target P is in slow moment time of zero doppler position, ��_cRepresent zero doppler position of target P and the slow time migration of wave beam central position,Represent i-th transmit/receive antenna and and target P between distance.

So far, launch and accept process once is introduced complete.

In addition, the needs of choosing of PRF ensure can not to produce in data gathering process that orientation is fuzzy and range ambiguity. Radar data obtains to terminate ordering and sends by user, or automatically terminate after the observation of certain target stripe scene is terminated, radar data obtains the judgement condition that process terminates: if described multi-angle synthetic-aperture radar needs to obtain the target stripe scene multi-angle information that a certain length is L, then data acquisition time T should be greater than and is oriented to �� from first₁Wave beam start this band scene inswept, N number of be oriented to �� to_NWave beam leave time of this band scene completely, namely

$T>\frac{L+2R_c\tan \left(\frac{{\mathrm{\&theta;}}_1-{\mathrm{\&theta;}}_N+\mathrm{\&Delta;\&theta;}}{2}\right)}{V}---\left(5\right)$

Wherein, R_cFor Texas tower is to the shortest oblique distance of band scene;

After radar data obtains process termination, i-th receiving antenna R_{X, i}The scatter echo data received are r via the data after down coversion and AD sampling_i(��, t; ��_i, �� ��), N number of such data form raw radar data collection R={r₁(��, t; ��₁, �� ��) ... r_i(��, t; ��_i, �� ��) ... r_N(��, t; ��_N, �� ��) }, raw radar data collection comprise same target by multiple wave beam with different visual angles irradiate after scattered information, be illustrated in figure 4 building target by different visual angles wave beam irradiate schematic diagram;

Step C: to each receiving antenna R corresponding in raw radar data collection R_{X, i}Two-dimentional echo data part r_i(��, t; ��_i, �� ��), isolate only by transmitting antenna T according to the mode demodulation of orthogonal coding_{X, i}Launch, received antenna R after target scattering_{X, i}The part receivedForm new echo data collectionWherein $\hat{R}=\{{\hat{r}}_1(\mathrm{\&eta;},t;{\mathrm{\&theta;}}_1,\mathrm{\&Delta;\&theta;}),...{\hat{r}}_i(\mathrm{\&eta;},t;{\mathrm{\&theta;}}_i,\mathrm{\&Delta;\&theta;}),...{\hat{r}}_N(\mathrm{\&eta;},t;{\mathrm{\&theta;}}_N,\mathrm{\&Delta;\&theta;})\};$

To each receiving antenna R corresponding in raw radar data collection R_{X, i}Two-dimentional echo data part r_i(��, t; ��_i, �� ��), owing to it not only comprises i-th back scattered s of beam coverage internal object_iT () signal echo, also comprises adjacent i.e. the i-th-1 and the i-th+1 wave beam and s that the territory scattering in of i-th beam overlap area is returned_i-1(t) and s_i+1(t) signal echo, and the echo composition that other beam side lobe causes. Therefore, in this step, r is isolated in mode demodulation according to orthogonal coding_i(��, t; ��_i, �� ��) in only by transmitting antenna T_{X, i}Launch, by R after target scattering_{X, i}The part receivedThe new echo data collection of final formation

Step D: to new echo data collectionIn 2-D data after each separationUse two-dimensional imaging algorithm to carry out imaging processing, obtain the band two dimension radar image g of N number of same target area_i(x, y; ��_i, �� ��), form image collection G, wherein G={g₁(x, y; ��_i, �� ��) ... g_i(x, y; ��_i, �� ��) ... g_N(x, y; ��_i, �� ��) };

Wherein, two-dimensional imaging algorithm can be distance range and Doppler, Chirp-Scaling algorithm or range migration algorithm.

Step e: to the band two dimension radar image g of the N number of same target area in image collection G_i(x, y; ��_i, �� ��) and carry out image registration, obtain images after registration collection $\hat{G}=\{{\hat{g}}_1(x,y;{\mathrm{\&theta;}}_i,\mathrm{\&Delta;\&theta;}),...{\hat{g}}_i(x,y;{\mathrm{\&theta;}}_i,\mathrm{\&Delta;\&theta;}),...{\hat{g}}_N(x,y;{\mathrm{\&theta;}}_i,\mathrm{\&Delta;\&theta;})\};$

Wherein, the method for image registration can adopt correlation coefficient process of the prior art, coefficient of coherence method etc.

Step F, it may also be useful to multi-angle information fusion and extracting method image collection after registrationMiddle extraction target stripe scene range of scatter angles [��_N-�� ��/2, ��₁+ �� ��/2] multi-angle scattered information, obtain target stripe scene two dimension radar image M (x, y) by this multi-angle scattered information.

Wherein, multi-angle information fusion and extracting method extract target stripe scene range of scatter angles [�� by means such as image co-registration, statistics, optimizations_N-�� ��/2, ��₁+ �� ��/2] multi-angle scattered information. In addition, what obtain that target stripe scene two dimension radar image M (x, y) adopts by multi-angle scattered information is the method for prior art, and the method is well known to those skilled in the art, illustrates no longer in detail herein.

So far, by reference to the accompanying drawings the present embodiment has been described in detail. Describing according to above, the method that multi-angle data of synthetic aperture radar of the present invention obtains should have been had and clearly recognized by those skilled in the art.

In addition, above-mentioned various concrete structures, shape or the mode that the definition of each element and method is not limited in mention in embodiment, it can be carried out simply changing or replacing by those of ordinary skill in the art.

In sum, multiple transmit/receive antennas that the method that multi-angle data of synthetic aperture radar of the present invention obtains is pointed to by adopting different azimuth, form multiple antenna irradiation wave beam, overcome tradition band synthetic-aperture radar wave beam orientation and point to the single limitation of scattering data angle information fixing, that obtain, can obtaining the scattered information of target scene scatter point multi-angle, scattered information is abundanter. In addition, be undertaken by striped patterns to the acquisition of multi-angle scattering data, it is possible to obtain comprise multi-angle scattered information on a large scale, continuous print target stripe scene two dimension radar image, range of application is more extensive.

Above-described specific embodiment; the object of the present invention, technical scheme and useful effect have been further described; it is it should be understood that; the foregoing is only specific embodiments of the invention; it is not limited to the present invention; within the spirit and principles in the present invention all, any amendment of making, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (4)

1. the method that a multi-angle data of synthetic aperture radar obtains, it is characterised in that, comprising:

Steps A: along orientation to X on Texas tower, arranges the transmit/receive antenna T that row are total to N number of transmit-receive sharing mode_X,i/R_X,i, make Texas tower with speed V along orientation to carrying out linear uniform motion, wherein, T_X,iRepresent that i-th transmit-receive sharing Antenna Operation is in transmitting antenna state, R_X,iRepresent that i-th transmit-receive sharing Antenna Operation is in receiving antenna state, i=1,2 ..., N;

Step B: Texas tower often moves V/PRF interval along orientation to X, N number of transmit/receive antenna T_X,i/R_X,iJust carrying out primary emission receiving course, the scatter echo data that N number of receiving antenna receives respectively form raw radar data collection R after treatment:

R={r₁(��, t; ��₁,����),��,r_i(��, t; ��_i,����),��r_N(��, t; ��_N,����)}

Wherein, r_i(��, t; ��_i, �� ��) and it is i-th receiving antenna R_X,iThe scatter echo data received are via the data after down coversion and AD sampling, and PRF is the pulse-repetition frequency of radar work;

Step C: to each receiving antenna R corresponding in raw radar data collection R_X,iTwo-dimentional echo data part r_i(��, t; ��_i, �� ��), isolate only by transmitting antenna T according to the mode demodulation of orthogonal coding_X,iLaunch, received antenna R after target scattering_X,iThe part receivedForm new echo data collectionWherein $\hat{R}=\{{\hat{r}}_1(\mathrm{\&eta;},t;{\mathrm{\&theta;}}_1,\mathrm{\&Delta;}\mathrm{\&theta;}),...{\hat{r}}_i(\mathrm{\&eta;},t;{\mathrm{\&theta;}}_i,\mathrm{\&Delta;}\mathrm{\&theta;}),...{\hat{r}}_N(\mathrm{\&eta;},t;{\mathrm{\&theta;}}_N,\mathrm{\&Delta;}\mathrm{\&theta;})\};$

Step D: to new echo data collectionIn 2-D data after each separationUse two-dimensional imaging algorithm to carry out imaging processing, obtain the band two dimension radar image g of N number of same target area_i(x, y; ��_i, �� ��), form image collection G:

G={g₁(x, y; ��_i,����),��g_i(x, y; ��_i,����),��g_N(x, y; ��_i,����)}

Step e: to the band two dimension radar image g of the N number of same target area in image collection G_i(x, y; ��_i, �� ��) and carry out image registration, obtain images after registration collection

$\hat{G}=\{{\hat{g}}_1(x,y;{\mathrm{\&theta;}}_i,\mathrm{\&Delta;}\mathrm{\&theta;}),...{\hat{g}}_i(x,y;{\mathrm{\&theta;}}_i,\mathrm{\&Delta;}\mathrm{\&theta;}),...{\hat{g}}_N(x,y;{\mathrm{\&theta;}}_i,\mathrm{\&Delta;}\mathrm{\&theta;})\};$ And

Step F, it may also be useful to multi-angle information fusion and extracting method image collection after registrationMiddle extraction target stripe scene range of scatter angles [��_N-����/2,��₁+ �� ��/2] multi-angle scattered information, obtain target stripe scene two dimension radar image M (x, y) by this multi-angle scattered information

Wherein, in described step B, launch and accept process comprises each time:

Sub-step B1, in slow time-sampling moment ��, N number of transmitting antenna forms different azimuth and points to ��₁,����_i,����_N, N number of beam of an antenna of bandwidth �� ��;

Wherein ��_iBeing the sensing angle, wave beam center of i-th antenna irradiation wave beam, �� �� is its bandwidth, meets ��_i=��_i-1-�� ��, i=1,2 ..., N, ��_i> 0 expression forward sight, ��_i< 0 represent after look, so the total visual angle range Theta �� of formation [-| ��_N|-����/2,��₁+ �� ��/2], ��₁And ��_NThe orientation being respectively the 1st transmit/receive antenna and N number of transmit/receive antenna is pointed to;

Sub-step B2, the pulse signal s with orthogonal encoding characteristics launched respectively by N number of transmitting antenna₁(t),��s_i(t),��s_N(t), this pulse signal s₁(t),��s_i(t),��s_NT () meets:

$\underset{0}{\overset{T_P}{\&Integral;}}{s}_i\left(t\right)s_j^{*}(t+\mathrm{\&tau;})dt={\mathrm{\&delta;}}_{ij},$ When ��=0

Wherein, T_pFor transponder pulse time width, �� is time delay variable, ��_ijFor Kronecker punching swashs function, it is defined as:

Sub-step B3, N number of receiving antenna R_X,1,��R_X,i,��R_X,NReceive the scatter echo data of the described pulse signal of Same Scene or target simultaneously;

Sub-step B4, the scatter echo data received by each receiving antenna carry out down coversion and AD sampling, and result are saved in raw radar data collection R:

R={r₁(��, t; ��₁,����),��,r_i(��, t; ��_i,����),��r_N(��, t; ��_N,����)}(3)

Wherein, r_i(��, t; ��_i, �� ��) and it is the data after the echo data that i-th receiving antenna receives is sampled via down coversion and AD, its expression formula is:

$r_i(\mathrm{\&eta;},t;{\mathrm{\&theta;}}_i,\mathrm{\&Delta;}\mathrm{\&theta;})=\underset{{\mathrm{\&Omega;}}_i}{\&Integral;}{\mathrm{\&sigma;}}_P({\mathrm{\&theta;}}_i,\mathrm{\&Delta;}\mathrm{\&theta;})w_a(\mathrm{\&eta;}-{\mathrm{\&eta;}}_0+{\mathrm{\&eta;}}_c)s_i(t-\frac{2R_P^i}{c})\exp (-j\frac{4{\mathrm{\&pi;R}}_P^i}{\mathrm{\&lambda;}})dxdydz$

Wherein, ��_iRepresent the wave cover district of i-th wave beam, ��_P(��_i, �� ��) and it is scattering coefficient when target P is irradiated by i-th wave beam, c is electromagnetic wave propagation speed, and �� is radar work mid-frequency, w_aRepresent that orientation is to antenna pattern, ��₀Represent that target P is in slow moment time of zero doppler position, ��_cRepresent zero doppler position of target P and the slow time migration of wave beam central position,Represent i-th transmit/receive antenna and and target P between distance;

Wherein, radar data obtains to terminate ordering and sends by user, or automatically terminate after the observation of certain target stripe scene is terminated, radar data obtains the judgement condition that process terminates: if described multi-angle synthetic-aperture radar needs to obtain the target stripe scene multi-angle information that a certain length is L, then data acquisition time T should be greater than and is oriented to �� from first₁Wave beam start this band scene inswept, N number of be oriented to �� to_NWave beam leave time of this band scene completely, namely

$T>\frac{L+2R_{c}tan\left(\frac{{\mathrm{\&theta;}}_1-{\mathrm{\&theta;}}_N+\mathrm{\&Delta;}\mathrm{\&theta;}}{2}\right)}{V}$

Wherein, R_cFor Texas tower is to the shortest oblique distance of band scene.

2. method according to claim 1, it is characterised in that, the two-dimensional imaging algorithm in described step D is distance range and Doppler, Chirp-Scaling algorithm or range migration algorithm.

3. method according to claim 1, it is characterised in that, in described step e, the method for image registration is correlation coefficient process or coefficient of coherence method.

4. method according to claim 1, it is characterised in that, described transmit/receive antenna is that single loudspeaker antenna or microstrip antenna form antenna submatrix.