DESIGN METHOD OF A TRUST-REGION SPACEBORNE DUAL-BASELINE INTERFEROMETRIC SAR SYSTEM

Abstract

A design method of a trust region satellite-borne dual-base interference SAR system is provided, relating to the field of radar measurement. Firstly, error sources of the interference SAR system are classified and analyzed, phase errors caused by image decoherence are calculated according to wave positions, and an analytical expression of baseline-related errors is reserved; a maximum phase error caused by imaging processing and an electronic device are estimated; and a baseline error is estimated. Then, the phase errors and the baseline errors are converted into elevation errors through an elevation fuzzy number, thereby obtaining an analytical expression of the elevation errors and the baseline parameters (length and dip angle). With a relative height measurement accuracy index of the system as an optimization object, a feasible solution interval between the baseline length and the baseline dip angle of the system is solved through a trust region algorithm. A minimum value of the baseline length is corrected according to a flying-around safety distance of a dual-satellite formation. According to the method of the present disclosure, the design of the satellite-borne dual-base interference SAR system can be guided, so that the satellite-borne dual-base interference SAR system is ensured to carry out terrain elevation surveying and mapping tasks with high accuracy and high reliability.

Claims

1. A design method of a trust region satellite-borne dual-base interference SAR system, comprising the following steps: step 1, calculating phase errors caused by image decoherence of the satellite-borne dual-base interference SAR system according to wave positions, and reserving an analytical expression of errors related to baseline parameters, wherein the baseline parameters comprise a baseline length and a baseline dip angle; step 2, estimating a maximum phase error caused by imaging processing of the satellite-borne dual-base interference SAR system; step 3, estimating a maximum phase error caused by an electronic device of the satellite-borne dual-base interference SAR system; step 4, estimating baseline errors of the satellite-borne dual-base interference SAR system; step 5, converting the phase errors caused by image decoherence of the satellite-borne dual-base interference SAR system in step 1, the maximum phase error caused by imaging processing of the satellite-borne dual-base interference SAR system in step 2, the maximum phase error caused by the electronic device of the satellite-borne dual-base interference SAR system in step 3 and the baseline errors of the satellite-borne dual-base interference SAR system in step 4 into elevation errors through an elevation fuzzy number to obtain an analytical expression of a relationship between the elevation errors and the baseline parameters; step 6, with a relative height measurement accuracy index of the system as an optimization object, solving the analytical expression of the relationship between the elevation errors and the baseline parameters obtained in step 5 through a trust region optimization algorithm to obtain a feasible solution interval between the baseline length and the baseline dip angle; and step 7, correcting a minimum value of the baseline length according to a flying-around safety distance of a dual-satellite formation.

2. The design method of the trust region satellite-borne dual-base interference SAR system according to claim 1, wherein the image decoherence of the satellite-borne dual-base interference SAR system in step 1 comprises quantization noise decoherence, Doppler decoherence, sidelobe decoherence, registration error decoherence, ambiguity decoherence, thermal noise decoherence, volume scattering decoherence and baseline decoherence; decoherence coefficients are calculated for typical terrain slopes and typical landforms at each wave position, and then a minimum decoherence coefficient of each wave position is taken for subsequent analysis; the typical terrain slopes are 0 slopes, 20% front slope surfaces, and 20% back slope surfaces; and the typical landforms are rocky soil and forest.

3. The design method of the trust region satellite-borne dual-base interference SAR system according to claim 2, wherein the step 2 comprises: estimating a maximum phase error introduced by an imaging algorithm and a maximum phase error introduced by a multi-channel reconstruction process respectively, and adding the two errors to obtain the maximum phase error caused by imaging processing of the satellite-borne dual-base interference SAR system.

4. The design method of the trust region satellite-borne dual-base interference SAR system according to claim 3, wherein the step 3 comprises: estimating a phase error caused by multiple channels, a phase error caused by an internal calibration residual error and a residual phase error after phase synchronization respectively, and adding the three errors to obtain the maximum phase error caused by the electronic device of the satellite-borne dual-base interference SAR system.

5. The design method of the trust region satellite-borne dual-base interference SAR system according to claim 4, wherein the step 4 comprises: estimating a vertical baseline error and a parallel baseline error respectively according to a platform attitude error, an antenna phase center error and a GPS position measurement error.

6. The design method of the trust region satellite-borne dual-base interference SAR system according to claim 5, wherein the step 6 comprises: calculating an intersection of the feasible solution intervals of the baseline parameters for all wave positions to obtain a final baseline parameter interval.

Description

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS

[0018] FIG. 1 is a flowchart of a design method of a trust region satellite-borne dual-base interference SAR system of the present disclosure.

DETAILED DESCRIPTION

[0019] The following will clearly and completely describe the technical solutions in the embodiments of the present disclosure in conjunction with the accompanying drawings in the embodiments of the present disclosure. Obviously, the described embodiments are only some embodiments of the present disclosure, not all embodiments. Based on the embodiments in the present disclosure, all other embodiments derived by a person of ordinary skill in the art without creative efforts shall fall within the protection scope of the present disclosure.

[0020] In conjunction with one embodiment of the present disclosure, a design method of a trust region satellite-borne dual-base interference SAR system is provided. As shown in FIG. 1, the design method includes the following steps.

[0021] At step 101: phase errors caused by image decoherence of the satellite-borne dual-base interference SAR system are calculated according to wave positions, and an analytical expression of baseline-related errors is reserved. Baseline parameters include a baseline length and a baseline dip angle. The baseline-related errors indicate errors related to the baseline parameters.

[0022] The image decoherence mainly includes quantization noise decoherence, Doppler decoherence, sidelobe decoherence, registration error decoherence, ambiguity decoherence, thermal noise decoherence, volume scattering decoherence, baseline decoherence and other factors. Some decoherence factors are also related to system factors, scene topography and terrain slopes. In the event of performing specific analysis, decoherence coefficients may be calculated for typical terrain slopes (0 slopes, 20% front slope surfaces, and 20% back slope surfaces) and typical landforms (rocky soil and forest) at each wave position (a central incident angle of the wave position is taken). Then, minimum decoherence coefficients of each wave position under six conditions are taken for subsequent analysis. The six conditions refer to conditions of three typical terrain slopes under two landforms, respectively.

[0023] Next, a calculation method for the decoherence of each image is given.

(1) Quantization Noise Decoherence:

[0024] Quantization noise decoherence is derived from the quantization of original recorded signals, which are usually approximated as Gaussian white noise, and a quantization noise decoherence coefficient .sub.SQNR may be expressed as (assuming that quantization noises of two images subjected to interference processing are consistent):

[00001]$\begin{array}{cc}{}_{\mathrm{SQNR}}=\frac{-1}{1+{\mathrm{SQNR}}^{-1}}& \left(1\right)\end{array}$

[0025] wherein, SQNR is a quantitative signal-to-noise ratio.

(2) Doppler Decoherence:

[0026] When interference data is acquired, Doppler frequency shift caused by a difference in beam pointing along an azimuth direction may also cause decoherence. A Doppler decoherence coefficient .sub.fdc may be expressed as:

[00002]$\begin{array}{cc}{}_{\mathrm{fdc}}=1-\frac{.Math.f_{\mathrm{dc}}.Math.}{B_d}& \left(2\right)\end{array}$

[0027] Wherein, f.sub.dc represents a Doppler center frequency difference between main and auxiliary images, and B.sub.d represents a Doppler bandwidth.

(3) Sidelobe Decoherence:

[0028] Since a point target imaging processing result is in the form of a sinc function, there will always be an impact of sidelobes on main lobes of other objects. In comprehensive consideration of the impacts of nearby objects, a sidelobe decoherence coefficient .sub.sidelobe may be expressed as (assuming that integral sidelobe ratios of two images subjected to interference processing are consistent):

[00003]$\begin{array}{cc}{}_{\mathrm{sidelobe}}=\frac{1}{1+\mathrm{ISLR}}& \left(3\right)\end{array}$

[0029] wherein, ISLR is an integral sidelobe ratio.

(4) Registration Error Decoherence:

[0030] A registration error decoherence coefficient .sub.coreg may be expressed as:

[00004]$\begin{array}{cc}{}_{\mathrm{coreg}}=\frac{\sin \left(\right)}{}& \left(4\right)\end{array}$

[0031] wherein, p represents the number of pixels with registration errors in an azimuth or range direction.

(6) Ambiguity Decoherence:

[0032] The ambiguity decoherence is caused by both range ambiguity RASR and azimuth ambiguity AASR. An ambiguity decoherence coefficient .sub.amb may be expressed as:

[00005]$\begin{array}{cc}{}_{\mathrm{amb}}=\frac{1}{1+RASR}.Math.\frac{1}{1+AASR}& \left(5\right)\end{array}$

[0033] The specific expressions of the range ambiguity and the azimuth ambiguity are relatively complicated. Considering that the ambiguity decoherence is not a dominant factor in the overall decoherence and the ambiguity problem can be well suppressed by an ambiguity suppression algorithm, in the present disclosure, typical values of the range ambiguity RASR and the azimuth ambiguity AASR are taken and set as constants. In addition, in view that an ambiguity decoherence coefficient of the rocky soil landform is smaller than that of a forest area under the same conditions, and the influence of decoherence is more serious, in the present disclosure, a typical value of ambiguity under the rocky soil landform is used to calculate the ambiguity decoherence.

(5) Thermal Noise Decoherence:

[0034] Thermal noise decoherence is affected by system parameters, scene topography and terrain slope. A thermal noise decoherence coefficient .sub.SNR may be expressed as:

[00006]$\begin{array}{cc}{}_{\mathrm{SNR}}=\frac{1}{\sqrt{\left(1+{\mathrm{SNR}}_1^{-1}\right)\left(1+{\mathrm{SNR}}_2^{-1}\right)}}& \left(6\right)\end{array}$

[0035] wherein, SNR.sub.1 and SNR.sub.2 represent signal-to-noise ratio of a main interference channel and a signal-to-noise ratio of an auxiliary interference channel, and a calculation formula is

[00007]$\mathrm{SNR}=\frac{{}_0\left({}_i-\right)}{\mathrm{NESZ}\left({}_i-\right)};{}_0$

is a backscattering coefficient; NESZ is a noise equivalent backscattering coefficient; .sub.i is a local incident angle; and is a current slope angle.

[0036] As can be seen from formula (6), 1) regardless of scene topography, the thermal noise decoherence is not obvious at a near end, and the decoherence is more serious at a far end due to the decrease in SNR; 2) the thermal noise decoherence on a front slope is better than that on the flat ground, and the thermal noise decoherence on the flat ground is better than that on a back slope; and 3) the thermal noise decoherence of rocky soil is better than that of forest. As described above, in the event of performing specific analysis, the thermal noise decoherence coefficients may be calculated for typical terrain slopes (0 slopes, 20% front slope surfaces, and 20% back slope surfaces) and typical landforms (rocky soil and forest) at each wave position (a central incident angle of this wave position is taken). Then, minimum thermal noise decoherence coefficients of each wave position under six conditions are taken for subsequent analysis.

(7) Volume Scattering Decoherence:

[0037] Volume scattering decoherence is determined by both system parameters and scene topography. For vegetation areas, when electromagnetic waves propagate in scatterers, decoherence caused by volume scattering is determined by a vegetation scattering elevation and an attenuation function of the electromagnetic waves. An expression of a volume scattering decoherence coefficient .sub.vol is:

[00008]$\begin{array}{cc}{}_{\mathrm{vol}}=\underset{0}{\overset{h_v}{}}{}^0\left(z\right).Math.\exp \left(j2\frac{z}{h_{\mathrm{amb}}}\right).Math.\mathrm{dz}/\underset{0}{\overset{h_v}{}}{}^0\left(z\right).Math.\mathrm{dz}& \left(7\right)\end{array}$

[0038] wherein, z is an integration variable; h.sub.v is a vegetation height; exp( ) represents an exponential function; j represents an imaginary unit; and .sup.0(z) is an attenuation function of vegetation to the electromagnetic waves, expressed as:

[00009]$\begin{array}{cc}{}^0\left(z\right)=\exp \left[-2.Math..Math.\frac{h_{}-z}{\cos \left({}_i\right)}\right]& \left(8\right)\end{array}$

[0039] wherein, is a one-way extinction coefficient of the electromagnetic waves in the vegetation, which is 0.2 dB/m (i.e., 0.023 Np/m) for an L band; .sub.i is a local incident angle (which can be directly calculated based on orbit parameters); and h.sub.amb is an ambiguity height, expressed as:

[00010]$\begin{array}{cc}{h}_{\mathrm{amb}}=\frac{R\sin {}_i}{B_{}}=\frac{R\sin {}_i}{B\cos \left({}_0-\right)}& \left(9\right)\end{array}$

[0040] wherein, is a wavelength; R is a slant range of a main satellite; B.sub.=B cos(.sub.0) is a vertical baseline; B is a length of a baseline on a cross-orbit plane; .sub.0 is a downward viewing angle corresponding to a reference surface (which can be directly calculated based on the orbit parameters); and is a baseline dip angle.

[0041] As can be seen from formula (9), a vegetation height h.sub.v affects volume scattering decoherence. Therefore, a volume scattering decoherence effect of the forest landform is much higher than that of the rocky soil landform, while a volume scattering decoherence effect of the rocky soil landform is almost negligible. In addition, an ambiguity height h.sub.amb also affects the volume scattering decoherence. At one wave position, in the present disclosure, a relationship between the ambiguity height h.sub.amb and the volume scattering decoherence coefficient .sub.vol is fitted through a quadratic polynomial, and the volume scattering decoherence coefficient .sub.vol may be expressed as:

[00011]$\begin{array}{cc}{}_{\mathrm{vol}}=p_2{h}_{\mathrm{amb}}^2+p_1{h}_{\mathrm{amb}}^1+p_0& \left(10\right)\end{array}$

[0042] wherein, p.sub.0-p.sub.2 are fitted polynomial coefficients, and a constant value may be taken at each wave position. h.sub.amb.sup.1 and h.sub.amb.sup.2 represent two ambiguity heights.

(8) Baseline Decoherence:

[0043] The existence of a baseline causes a difference in incident angles of main and auxiliary antennas, resulting in baseline decoherence. The longer a vertical baseline is, the more serious the decoherence is. Baseline decoherence is related to system parameters and terrain slope. A baseline decoherence coefficient .sub.baseline may be expressed as:

[00012]$\begin{array}{cc}{}_{\mathrm{baseline}}=1-\frac{B_{}}{B_{,\mathrm{crit}}}& \left(11\right)\end{array}$

[0044] wherein, B.sub.=B cos(.sub.0) is a vertical baseline; and B.sub.,crit is an extreme baseline, i.e., a maximum vertical baseline with coherence between main and auxiliary antennas, which can be expressed as:

[00013]$\begin{array}{cc}{B}_{,\mathrm{crit}}=\frac{2B_{\mathrm{rg}}R\tan \left({}_i-\right)}{c}& \left(12\right)\end{array}$

[0045] wherein, B.sub.rg is a bandwidth in a range direction; is a wavelength; R is a slant range of a main satellite; .sub.i is a local incident angle; r is a local slope angle; and c is a velocity of light.

[0046] All the above decoherence factors are multiplied to obtain an overall decoherence coefficient :

[00014]$\begin{array}{cc}\begin{array}{c}={}_{\mathrm{SQNR}}.Math.{}_{\mathrm{fdc}}.Math.{}_{\mathrm{sidelobe}}.Math.{}_{\mathrm{coreg}}.Math.{}_{\mathrm{amb}}.Math.{}_{\mathrm{SNR}}.Math.{}_{\mathrm{vol}}.Math.{}_{\mathrm{baseline}}\\ ={}_{\mathrm{other}}.Math.\left(p_2.Math.h_{\mathrm{amb}}^2+p_1.Math.h_{\mathrm{amb}}^1+p_0\right).Math.\left(1-\frac{B_{}}{B_{,\mathrm{crit}}}\right)\\ ={}_{\mathrm{other}}.Math.\left[p_2.Math.{\left(\frac{R\sin {}_i}{B_{}}\right)}^2+p_1.Math.\left(\frac{R\sin {}_i}{B_{}}\right)+p_0\right].Math.\\ \left(1-\frac{{\mathrm{cB}}_{}}{2B_{\mathrm{rg}}R\tan \left({}_i-\right)}\right)\\ ={}_{\mathrm{other}}.Math.\left[p_2.Math.{\left(\frac{R\sin {}_i}{B_{}}\right)}^2+p_1.Math.\left(\frac{R\sin {}_i}{B_{}}\right)+p_0\right].Math.\\ \left(1-p_{-1}{B}_{}\right)\\ ={}_{\mathrm{other}}.Math.\left(p_2^{}\frac{1}{B_{}^2}+p_1^{}\frac{1}{B_{}}+p_0^{}+p_{-1}^{}{B}_{}\right)\\ ={}_{\mathrm{other}}.Math.(p_2^{}\frac{1}{{\left[B\cos \left({}_0-\right)\right]}^2}+p_1^{}\frac{1}{B\cos \left({}_0-\right)}+\\ p_0^{}+p_{-1}^{}B\cos \left({}_0-\right))\end{array}& \left(13\right)\end{array}$

[0047] wherein, intermediate parameters are .sub.other=.sub.SQNR.Math..sub.fdc.Math..sub.sidelobe.Math..sub.coreg.Math..sub.amb.Math..sub.SNR and

[00015]$p_{-1}=\frac{c}{2B_{\mathrm{rg}}R\tan \left({}_i-\right)};$

and all parameters with a superscript are intermediate parameters and have no meaning. p.sub.2=p.sub.2.Math.(R sin .sub.i).sup.2, p.sub.1=p.sub.1.Math.(R sin .sub.i)p.sub.2p.sub.1, p.sub.0=p.sub.0p.sub.1p.sub.1(R sin .sub.i) and p.sub.1=p.sub.0p.sub.1. At each wave position, except for two parameters of baseline length and baseline dip angle to be optimized in formula (13), all other parameters are known.

[0048] In a multi-view case, all the above decoherence factors may be converted into an interference phase error .sub.coherence caused by image decoherence according to the following formula:

[00016]$\begin{array}{cc}{}_{c\mathrm{oherence}}=\sqrt{\frac{1-{}^2}{2L{}^2}}& \left(14\right)\end{array}$

[0049] wherein, L is a multi-view number, which may also change depending on wave positions.

[0050] At step 102: a maximum phase error caused by imaging processing of the satellite-borne dual-base interference SAR system is estimated.

(1) Imaging Algorithm Error:

[0051] A two-dimensional spatial variation of satellite-borne dual-base SAR imaging echoes greatly increases the difficulty of imaging processing. Different dual-base SAR imaging algorithms have different two-dimensional spectrum accuracy, which will introduce different degrees of phase errors and ultimately affect the interference performance of the system.

(2) Multi-Channel Reconstruction Error:

[0052] In the process of dual-satellite interference height measurement, since a single satellite uses an azimuth-direction dual-channel technology to improve an azimuth-direction resolution during data reception, certain phase errors will be introduced in the dual-channel reconstruction process.

[0053] A maximum phase error introduced by imaging processing may be obtained by adding a maximum phase error introduced by an imaging algorithm and a maximum phase error introduced by the multi-channel reconstruction process, which is denoted as .sub.algorithm.

[0054] At step 103: a maximum phase error caused by an electronic device of the satellite-borne dual-base interference SAR system is estimated.

[0055] Phase errors caused by the electronic device may mainly be derived from residual errors after the correction of the electronic device, which may be divided into three categories. [0056] (1) A dual-satellite system adopts a multi-channel mode of one transmission and dual receptions in an azimuth direction, so there will be multi-channel phase errors introduced by a channel phase center position deviation, an attitude measurement error and a phase error between channels themselves. In summary, a total phase error caused by multiple channels may be obtained from the above three error sources. [0057] (2) The system response of a transmitting and receiving loop of the system will change with temperature, and the instability of a transmitting and receiving channel of the system needs to be corrected through a calibration loop, so there will be the phase error introduced by the calibrated residual error. [0058] (3) Since the dual-satellite system has phase errors caused by the instability of frequency sources due to the use of different frequency sources in main and auxiliary satellites, so phase synchronization is required. The phase error will be introduced by the residual phase error after phase synchronization.

[0059] A maximum phase error caused by the electronic device may be obtained by adding maximum phase errors caused by the three categories of residual errors corrected by the electronic device, which is denoted as .sub.equipment.

[0060] At step 104: a baseline error of the satellite-borne dual-base interference SAR system is estimated.

[0061] The baseline measurement accuracy is an important factor affecting interference height measurement, and this error is mainly derived from a platform attitude error, an antenna phase center error and a GPS position measurement error, and is manifested in the form of a large fixed deviation and a small random error. The fixed deviation mainly affects the absolute height measurement accuracy, and can be removed by accurate DEM correction. A residual baseline error mainly affects the relative height measurement accuracy, and is manifested in the form of a random error. B.sub. is denoted as a parallel baseline measurement error of the system, and B.sub. is denoted as a vertical baseline measurement error of the system.

[0062] At step 105: the phase errors and the baseline errors are converted into elevation errors through an elevation fuzzy number. The phase errors include a phase error caused by decoherence, a phase error caused by imaging processing and a phase error caused by the electronic device respectively calculated in step 101 to step 103. The baseline errors are obtained from step 104.

[0063] The phase error caused by the image decoherence obtained in step 101, the maximum phase error caused by imaging processing obtained in step 102 and the maximum phase error caused by the electronic device obtained in step 103 can be converted into elevation errors by multiplying an elevation fuzzy number. Then, the parallel baseline error and the vertical baseline error obtained in step 104 are also converted into elevation errors through an elevation fuzzy number. It is noted that, for the sake of simplicity, the present disclosure does not consider the baseline errors in the event of analyzing the elevation errors caused by the first three. Therefore, a total elevation error Ah caused by all error sources is:

[00017]$\begin{array}{cc}h=\frac{h_{\mathrm{amb}}}{2}\left({}_{\mathrm{coherence}}+{}_{\mathrm{algorithm}}+{}_{\mathrm{equipment}}\right)+\frac{h_{\mathrm{amb}}}{}B+h.Math.\frac{B_{}}{B_{}}& \left(15\right)\end{array}$

[0064] As can be seen from formula (13) and formula (14), the interference phase error .sub.coherence caused by image decoherence is approximated as a vertical baseline increasing function, while the ambiguity height h.sub.amb is a vertical baseline decreasing function. In addition, the elevation error caused by the baseline error is a decreasing function of the vertical baseline. Therefore, Equation (15) mainly compromises the vertical baseline B.sub. from two aspects: a phase error and an elevation error.

[0065] Formula (13) is written as =P(B.sub.) to represent that is a polynomial function of B.sub., and then after the substitution of formula (9), as well as the abbreviated formula (13) and formula (14), this equation may be written as:

[00018]$\begin{array}{cc}h=\frac{R\sin {}_i}{2B_{}}\left(\sqrt{\frac{1-{P\left(B_{}\right)}^2}{2\mathrm{LP}{\left(B_{}\right)}^2}}+{}_{\mathrm{algorithm}}+{}_{\mathrm{equipment}}\right)+\frac{R\sin {}_i}{B_{}}{B}_i+h.Math.\frac{B_{}}{B_{}}& \left(16\right)\end{array}$

[0066] It can be seen that this equation is a height nonlinear equation with respect to the vertical baseline B.sub. and therefore also a height nonlinear equation with respect to the baseline length B and the baseline dip angle . When derivatives of the analytic expression about B and a are calculated, a relationship between the two parameters and B.sub. should be solved as an implicit function to greatly simplify the operation.

[0067] At step 106: with a relative height measurement accuracy index of the system as an optimization target, a feasible solution interval between the baseline parameters of the system is solved through a trust region optimization algorithm.

[0068] That is, the analytical expression of the relationship between the elevation errors and the baseline parameters obtained in step 5 is solved through the trust region optimization algorithm to obtain a feasible solution interval between the baseline length and the baseline dip angle of the system.

[0069] According to the analytic expression (15) or (16) of the elevation errors and the baseline length as well as baseline dip angle obtained in step 105, a parameter interval that meets requirements can be solved with h being less than a relative elevation accuracy index of the system as an optimization object. Because the expression of formula (15) or formula (16) is highly nonlinearized and contains a plurality of parameters, the trust domain optimization algorithm is adopted in the present disclosure for solving. After the feasible parameter intervals of all wave positions are solved one by one, an interval intersection is calculated to obtain a final parameter interval (different parameters may also be used for different wave positions).

[0070] At step 107: a minimum value of the baseline length is corrected according to a flying-around safety distance of a dual-satellite formation.

[0071] Since the height measurement task is performed in a flying-around mode of the dual-satellite formation, it is necessary to ensure that the baseline length is greater than a shortest safety distance (denoted D) between satellites. Therefore, a minimum baseline parameter value (denoted as B.sub.min) solved in step 106 may be corrected as:

[00019]$\begin{array}{cc}{B}_{\min }=\{\begin{array}{c}D,B_{\min }<D\\ B_{\min },B_{\min }D\end{array}& \left(17\right)\end{array}$

[0072] The technical solutions of the present disclosure will be further described below in conjunction with the specific embodiments.

Embodiment 1

[0073] A Lutan-1 (LT-1) L-band differential interference SAR satellite consists of two 3 m-resolution L-band synthetic aperture radar satellites with the same design status. The two satellites are sent into an orbit by two CZ-4C rockets in the mode of one rocket and one satellite, and after entering the orbit, they will fly in turn according to a flying-around mode and a follow-up mode. An engineering task of LT-1 is to operate in a predetermined orbit; establish a ground system for the reception, transmission, processing, distribution and calibration of corresponding satellite remote sensing data; establish application demonstration conditions; promote the construction of multi-user application systems; achieve the functions of single/dual satellite heavy orbit differential interferometric deformation measurement, dual-satellite orbit interference height measurement, single/multi-polarization imaging in accordance with multi-mission planning requirements; complete observation tasks such as surface deformation measurement and digital elevation model acquisition in designated areas; support and serve territorial resources, earthquakes, disaster prevention and mitigation, basic geographic information acquisition, forestry and civil-military integration and other application needs; and form a rapid response capability for major disasters.

[0074] According to the design method of dual-base interference SAR baseline parameters in the specific implementation, the system baseline parameters are designed in combination with a relative elevation accuracy index of less than 5 m in LT-1. Image decoherence is calculated first. SQNR=18 dB is taken to obtain .sub.SQNR=0.98; f.sub.dc=11 Hz and B.sub.d=2426 Hz are obtained according to a 0.01 beam pointing error, and .sub.fdc=0.97 is further obtained; ISLR=15 dB is taken to obtain .sub.sidelobe=0.97; the high-accuracy registration is usually considered as that is 0.1 pixel, and a coherence coefficient obtained by considering two registration directions in a range direction and an azimuth direction is .sub.coreg=0.97; for ambiguity decoherence, an approximate value is .sub.amb=0.995; and the thermal noise decoherence is related to system factors, so a back slope of 20% and .sub.SNR under the forest landform are calculated under different wave positions, respectively. Volume scattering decoherence, as well as baseline decoherence, are reserved in the form of analytic expressions for parameter optimization.

[0075] A maximum phase error caused by imaging processing is calculated. According to the simulation analysis, a phase error value caused by Generic SAR Processor Algorithm is approximately 0.5, and the magnitude of a multi-channel reconstruction error is approximately 0.015.

[0076] Electronic device errors are further calculated according to error categories. Through the simulation analysis, when the channel phase center position deviation is 1 mm, a phase error after imaging is 0.8; and when an attitude measurement error is 0.003, the phase error after imaging is 0.13, and a phase error between channels themselves is 0.35, so a phase error caused by multiple channels is totally 1.28. For a single satellite, a residual error after internal calibration correction is approximately 0.7, and an error after imaging processing is approximately 0.7. A residual error after synchronous phase correction is approximately 2, and an error after imaging processing is 0.8.

[0077] For the baseline error, the simulation analysis verifies that a random error of GPS position measurement is 1 mm, a SAR phase center offset is 1 mm, and an attitude measurement error is 0.003; and meanwhile, in conjunction with spatial variation characteristics of a parallel baseline and a vertical baseline, measurement errors of the parallel baseline and the vertical baseline are B.sub.|=2 mm and B.sub.=2.24 mm, respectively.

[0078] Based on the above error sources, according to a requirement that the relative elevation accuracy of LT-1 is less than 5 m, a finally calculated vertical baseline range is about 1 km to 5.2 km.

[0079] Although the specific illustrative embodiments of the present disclosure are described above, the present disclosure may be understood by a person skilled in the art. In addition, it should be clear that the present disclosure is not limited to the scope of specific embodiments. For a person of ordinary skill in the art, as long as various changes are within the spirit and scope of the present disclosure limited and determined by the attached claims, these changes are obvious, and all disclosures and creations that make use of the present disclosure are protected.

[0080] The various embodiments described above can be combined to provide further embodiments. Aspects of the embodiments can be modified, if necessary to employ concepts of the various embodiments to provide yet further embodiments.

[0081] These and other changes can be made to the embodiments in light of the above-detailed description. In general, in the following claims, the terms used should not be construed to limit the claims to the specific embodiments disclosed in the specification and the claims, but should be construed to include all possible embodiments along with the full scope of equivalents to which such claims are entitled. Accordingly, the claims are not limited by the disclosure.