CN102680973A - Method for removing repeat-pass InSAR (Interferometric Synthetic Aperture Radar) nonsynchronous signals based on burst mode - Google Patents

Abstract

The invention discloses a method for removing repeat-pass InSAR (Interferometric Synthetic Aperture Radar) nonsynchronous signals based on a burst mode. The method is used for interfering removal of nonsynchronous signals in combination with azimuth time domain convolution imaging, and can be applied to interference processing of mode data such as ScanSAR and the like. The method comprises the following steps of: calculating the burst synchronism by using SAR original signal data, additional pass data and the like through a similar InSAR space base line calculating method, and setting a time node for distinguishing synchronous data from nonsynchronous data; and after azimuth time domain convolution imaging, extracting a phase only containing a synchronous signal convolution result for serving as a phase of an output picture element, and extracting an amplitude only containing a convolution result of synchronous signals and nonsynchronous signals for serving as an amplitude of the output picture element, so that amplitude information is ensured while nonsynchronous signals are removed. The method can further be used for processing a group of interference data; and under the condition, a group of synchronisms is required to be calculated and a group of time nodes is required to be set, so that a position of an imaging result corresponds to a group of output picture elements at last.

Description

Repeat track InSAR nonsynchronous signal removal method based on the burst pattern

Technical field

The present invention relates to a kind of interfering synthetic aperture radar and measure (InSAR) data processing method, specifically is a kind of nonsynchronous signal removal method of the repeat track InSAR data based on the burst pattern.

Background technology

Satellite-borne synthetic aperture radar in recent years (SAR) system almost all possesses the mode of operation based on burst, like the ScanSAR pattern, has expanded the spatial observation scope of SAR greatly.Band pattern with conventional is the same, and the data of burst pattern also can be used for interfering to be handled to obtain DEM or face of land deformation data.

Burst is the set of pulses that SAR launches and receives in flight course, be made up of mutual discontinuous burst one by one based on the data of burst pattern, thereby the orientation of terrain object is discontinuous to Doppler frequency spectrum.If the burst data are interfered processing; It is certain overlapping then to require the Doppler frequency spectrum of major-minor image burst to have; Promptly require major-minor image when data are obtained, certain scan-synchronized property will be arranged, this control to satellite orbit and attitude thereof has proposed very high requirement.In fact, more or less all there is certain non-synchronous data in the burst data of obtaining for repeat track SAR, and it also can bring noise to interferogram on the contrary to the not contribution of resolution of interferogram, so should it be removed to improve the quality of interferogram.Present burst nonsynchronous signal removal method mainly contains two kinds:

(1) in raw data, directly non-synchronous data is removed before the imaging.This method has also been removed its corresponding amplitude information when removing asynchronous phase information, be unfavorable for the image registration of back.

(2) remove nonsynchronous signal through filtering before interfering processing the imaging back, but the burst fringe region can produce noise because of the part convolution after the filtering.

Summary of the invention

The present invention proposes a kind of new method, is the burst nonsynchronous signal to be removed in imaging process in the orientation, can effectively overcome the problem that above-mentioned two kinds of methods exist.

Technical scheme of the present invention is following:

A kind of repeat track InSAR nonsynchronous signal removal method based on the burst pattern comprises the steps:

1) original SAR signal data of input and subsidiary orbital data thereof etc.;

2) utilize major-minor image corresponding SAR locus and speed, calculate the synchronism of major-minor image burst through the temporal information of interative computation and burst;

3) with step 2) in the synchronism that calculates and pre-set threshold relatively, if synchronism is higher than threshold value, so burst is carried out distance to compression with apart from proofreading and correct to migration; Then data are carried out the orientation to compression, simultaneously nonsynchronous signal is removed.Otherwise, abandon interfering and handle.

Described repeat track InSAR nonsynchronous signal removal method based on the burst pattern is characterized in that, in the step 1), can import a pair of SAR data, also can import one group of SAR data; Simultaneously, also need import the subsidiary orbital data of SAR data etc., comprising locus and the speed of SAR.

Described repeat track InSAR nonsynchronous signal removal method based on the burst pattern is characterized in that step 2) in, if handle a pair of SAR data, so only estimate this synchronism to data; If handle one group of SAR data; To estimate

individual synchronism so; X presentation video scape number wherein, x>2.

Described repeat track InSAR nonsynchronous signal removal method based on the burst pattern is characterized in that, in the step 3), process (as shown in Figure 2, corresponding to the operation in the figure l frame of broken lines) is as follows removed to compression and nonsynchronous signal in concrete orientation:

3.1) according to synchronism result of calculation timing node is set on burst, to distinguish synchronous and asynchronous data;

3.2) burst and orientation are multiplied each other to matched filter; And multiplied result added up; Get the phase place that only contains the synchronizing signal accumulation result phase place as the output pixel, get contain synchronizing signal and nonsynchronous signal accumulation result amplitude as the amplitude of exporting pixel.

Described repeat track InSAR nonsynchronous signal removal method based on the burst pattern is characterized in that step 3.1) in, if processing is a pair of SAR data, this timing node to data so only is set; If what handle is one group of SAR data, the right timing node of corresponding SAR data is set respectively then.

Described repeat track InSAR nonsynchronous signal removal method based on the burst pattern; It is characterized in that; Step 2) method of the synchronism of the major-minor image burst of said calculating is: for the data line among burst of master image, calculate its corresponding SAR locus (R _M) and speed, it is capable to seek burst corresponding in the auxilliary image through interative computation then; The capable pairing SAR position (P of burst when auxilliary image _S) on the given velocity reversal of master image with P _MThe pairing deviation of master image corresponding row stops iteration during less than pre-set threshold; Find out the pairing burst of major-minor image capable after, according to the zero-time of burst separately, calculate corresponding row time of burst zero-time relatively separately, can confirm the burst synchronism of major-minor image.

Advantage of the present invention mainly contains:

With respect to the method for directly in raw data, removing non-synchronous data, this method has kept whole amplitude informations, and this will help the registration between image; With respect to the method for imaging back through filtering removal nonsynchronous signal, the noise of bringing for the burst fringe region because of the part convolution when this method has been avoided filtering.

Aspect calculated amount, with respect to the method that remove nonsynchronous signal through filtering before interfering processing the back of forming images, when the processing single interference was right, because the time domain convolution algorithm plays a leading role, calculated amount of the present invention was the former several times; But the process mass interference data to the time, because of filtering plays a leading role, calculated amount of the present invention is on the contrary less than the former.

Description of drawings

Fig. 1 is for adopting the interference processing flow chart based on the burst pattern of the inventive method.

Fig. 2 is the concrete operations synoptic diagram in Fig. 1 frame of broken lines.

Embodiment

Below in conjunction with accompanying drawing the said technical scheme of summary of the invention is elaborated.Fig. 1 is for adopting the interference processing flow chart based on the burst pattern of the inventive method.Fig. 2 is the concrete operations synoptic diagram in Fig. 1 frame of broken lines.

When the original SAR signal data of input, handle equally with conventional interference, need to extract and interfere the parameters of using in the processing, also to calculate information such as doppler centroid simultaneously.

The synchronism of burst can be calculated through the subsidiary orbital data of SAR data, and its process is similar to the calculating of interfering Space Baseline: for the data line among burst of master image, calculate its corresponding SAR locus (P _M) and speed, it is capable to seek burst corresponding in the auxilliary image through interative computation then.The capable pairing SAR position (P of burst when auxilliary image _S) on the given velocity reversal of master image with P _MDeviation during less than pre-set threshold, stop iteration; Find out the pairing burst of major-minor image capable after, according to the zero-time of burst separately, calculate corresponding row time of burst zero-time relatively separately, can confirm the burst synchronism of major-minor image.

Compare with conventional band pattern signal, the burst signal with its distance to characteristics identical, the difference of the two be embodied in the orientation to.Be higher than the interference data that preestablishes threshold value for synchronism, burst carried out distance compression and distance after migration is proofreaied and correct, the present invention combines the orientation to carry out the nonsynchronous signal removal to imaging processing.

Ignore the orientation to antenna radiation pattern, the orientation of a burst can be described as to signal

$s_a(\mathrm{\η},{\mathrm{\η}}_{\mathrm{dn}},{\mathrm{\η}}_{\mathrm{cn}})=\underset{n=0}{\overset{+\∞}{\mathrm{\Σ}}}{s}_{\mathrm{an}}(\mathrm{\η},{\mathrm{\η}}_{\mathrm{dn}},{\mathrm{\η}}_{\mathrm{cn}})$ $\left(1\right)$

$=\underset{n=0}{\overset{+\∞}{\mathrm{\Σ}}}\mathrm{rect}\left(\frac{\mathrm{\η}-{\mathrm{\η}}_{\mathrm{cn}}}{T_b}\right)\·\exp \{-j\frac{4\mathrm{\π}\sqrt{R_0^2+V_r^2(\mathrm{\η}-{\mathrm{\η}}_{\mathrm{dn}})}}{\mathrm{\λ}}\}$

Wherein, η be the orientation to the time, η _DnBe zero Doppler's moment of point target, η _CnBe the burst center of point target, T _bBe the burst duration, R ₀Be low coverage, λ is a radar wavelength, V _rBe effective radar speed, n is a point target number, s _An() be a single point target direction to signal, rect () is a rectangular function.

The complete signal copy of a point target does

$g_{\mathrm{an}}(\mathrm{\η},{\mathrm{\η}}_{\mathrm{dn}},{\mathrm{\η}}_{\mathrm{dcn}})=h_{\mathrm{an}}^{*}(-\mathrm{\η},{\mathrm{\η}}_{\mathrm{dn}},{\mathrm{\η}}_{\mathrm{dcn}})$

$=\mathrm{rect}\left(\frac{\mathrm{\η}-{\mathrm{\η}}_{\mathrm{dcn}}}{T_a}\right)\·\exp \{j\frac{4\mathrm{\π}[R_c-\sqrt{R_0^2+V_r^2(\mathrm{\η}-{\mathrm{\η}}_{\mathrm{dn}})}]}{\mathrm{\λ}}\}---\left(2\right)$

Wherein, η _DcnFor the pairing orientation of beam center to the time, R _cBe the corresponding oblique distance of beam center, T _aBe the whole synthetic aperture time, h _An() is corresponding matched filter, and * representes conjugation.In the following formula with R _cThe phase place of bringing is regarded as constant term, and this phase term is the phase place of imaging back target, it is not included in the point target signal here, and all the other are approximately linear frequency modulation phase place.Through the time domain convolution algorithm burst is carried out the orientation to imaging, so can be through following relevant or convolution completion

${\∫}_{-\∞}^{+\∞}{s}_{\mathrm{an}}(u,{\mathrm{\η}}_{\mathrm{an}},{\mathrm{\η}}_{\mathrm{cn}})g_{\mathrm{an}}^{*}(u-\mathrm{\η},{\mathrm{\η}}_{\mathrm{dn}},{\mathrm{\η}}_{\mathrm{dcn}})\mathrm{du}$ $\left(3\right)$

$={\∫}_{-\∞}^{+\∞}{s}_{\mathrm{an}}(u,{\mathrm{\η}}_{\mathrm{dn}},{\mathrm{\η}}_{\mathrm{cn}})h_{\mathrm{an}}(\mathrm{\η}-u,{\mathrm{\η}}_{\mathrm{dn}},{\mathrm{\η}}_{\mathrm{dcn}})\mathrm{du}.$

The convolutional calculation result of a point adding up after to be exactly that matched filter and burst are point-to-point multiply each other, adding up is that pointwise is carried out.Suppose the synchronism of known major-minor image; So timing node can be set before adding up, to distinguish synchronously and asynchronous data, cumulative process begins from an end of synchrodata then; When cumulative process has reached this timing node; Accumulation result leaves among the variable p, adds up then and proceeds, and last accumulation result leaves among the variable m.Only comprised synchronous signal among the p, but its amplitude resolution has reduced; Both comprised synchronizing signal among the m, and also comprised nonsynchronous signal, its amplitude resolution is corresponding to whole burst bandwidth.What final imaging results needed is the phase place of p and the amplitude of m, so the amplitude of p is replaced to the amplitude of m, its advantage is both to have removed asynchronous signal in the interferometric phase, has kept the amplitude of full resolution (corresponding to whole burst bandwidth) again.Can find out that the extra computation of removing corresponding to nonsynchronous signal in the whole process only is the amplitude replacement operation.

At present; Often want the interferogram (Interferogram Stack) of process mass in the interventional applications, one group of SAR data is promptly arranged, each the scape data in the data all might be right with these other data composition interference of organizing in the data; To interfere processing, the present invention also can be used for this situation.At this moment, for other pairing interference data, the nonsynchronous signal of scape data both maybe be at the head end of burst, also maybe be at the tail end of burst.For scape SAR data, suppose to have k and l to interfering to carrying out the removal of nonsynchronous signal from burst head end and tail end respectively.Before adding up, a k and l timing node is set respectively, it is right to interfering to correspond respectively to above-mentioned k and l.Suppose to begin to add up from the burst head end of correspondence, add up after the completion, the part accumulation result can exist in the following vector

P=(p′ ₁，p′ ₂，…，p′ _k，p ₁，p ₂，…，p _l，) （4）

Wherein, (p ' ₁, p ' ₂..., p ' _k) and (p ₁, p ₂..., p _l) correspond respectively to k and l is right to interfering.And complete convolution results leaves among the variable m.Being used for next step so interferes the imaging results of handling to be expressed as

O=(P+A)B （5.1）

Wherein

For a burst, its calculated amount comprises that time domain convolution algorithm and the inferior amplitude of z (k+l) substitute to be calculated, wherein z be expressed as picture back burst the orientation to the pixel number.

Claims (6)

1. the repeat track InSAR nonsynchronous signal removal method based on the burst pattern comprises the steps:

1) original SAR signal data of input and subsidiary orbital data thereof;

2) utilize major-minor image corresponding SAR locus and speed, calculate the synchronism of major-minor image burst through the temporal information of interative computation and burst;

3) with step 2) in the synchronism that calculates and pre-set threshold relatively, if synchronism is higher than threshold value, so burst is carried out distance to compression with apart from proofreading and correct to migration; Then data are carried out the orientation to compression, simultaneously nonsynchronous signal is removed; Otherwise, abandon interfering and handle.

2. the repeat track InSAR nonsynchronous signal removal method based on the burst pattern as claimed in claim 1 is characterized in that, in the step 1), can import a pair of SAR data, also can import one group of SAR data; Simultaneously, also need import the subsidiary orbital data of SAR data, comprising locus and the speed of SAR.

3. like the described repeat track InSAR nonsynchronous signal removal method of claim l, it is characterized in that step 2 based on the burst pattern) in, if handle a pair of SAR data, so only estimate this synchronism to data; If handle one group of SAR data; To estimate

individual synchronism so; X presentation video scape number wherein, x>2.

4. the repeat track InSAR nonsynchronous signal removal method based on the burst pattern as claimed in claim 1 is characterized in that, in the step 3), concrete orientation is following to compression and nonsynchronous signal removal process:

3.1) according to synchronism result of calculation timing node is set on burst, to distinguish synchronous and asynchronous data;

3.2) burst and orientation are multiplied each other to matched filter; And multiplied result added up; Get the phase place that only contains the synchronizing signal accumulation result phase place as the output pixel, get contain synchronizing signal and nonsynchronous signal accumulation result amplitude as the amplitude of exporting pixel.

5. the repeat track InSAR nonsynchronous signal removal method based on the burst pattern as claimed in claim 4 is characterized in that step 3.1) in, if processing is a pair of SAR data, this timing node to data so only is set; If what handle is one group of SAR data, the right timing node of corresponding SAR data is set respectively then.

6. like the described repeat track InSAR nonsynchronous signal removal method of claim l based on the burst pattern; It is characterized in that; Step 2) method of the synchronism of the major-minor image burst of said calculating is: for the data line among burst of master image, calculate its corresponding SAR locus P _MAnd speed, it is capable to seek burst corresponding in the auxilliary image through interative computation then; The capable pairing SAR position P of burst when auxilliary image _SOn the given velocity reversal of master image with P _MThe pairing deviation of master image corresponding row stops iteration during less than pre-set threshold; Find out the pairing burst of major-minor image capable after, according to the zero-time of burst separately, calculate corresponding row time of burst zero-time relatively separately, can confirm the burst synchronism of major-minor image.

Images