METHOD AND DEVICE FOR SUPPRESSING RANGE AMBIGUITY IN SYNTHETIC APERTURE RADAR

Abstract

In a method and device for suppressing range ambiguity in Synthetic Aperture Radar (SAR), azimuth transmission-phase modulation is performed respectively, using a preset rule for phase modulation, on pulse signal data transmitted on a channel of Horizontal (H) polarization and pulse signal data transmitted on a channel of Vertical (V) polarization. Full polarimetric echo data corresponding to the transmission-phase modulated pulse signal data transmitted on the channel of H polarization and the transmission-phase modulated pulse signal data transmitted on the channel of V polarization are acquired. Azimuth phase demodulation is performed on the full polarimetric echo data using a preset rule for phase demodulation. Echo data are acquired by filtering out, using a preset azimuth filter, range-ambiguity energy from the phase demodulated full polarimetric echo data.

Claims

1. A method for suppressing range ambiguity in Synthetic Aperture Radar (SAR), comprising: performing, using a preset rule for phase modulation, azimuth transmission-phase modulation respectively on pulse signal data transmitted on a channel of Horizontal (H) polarization and pulse signal data transmitted on a channel of Vertical (V) polarization; acquiring full polarimetric echo data corresponding to the transmission-phase modulated pulse signal data transmitted on the channel of H polarization and the transmission-phase modulated pulse signal data transmitted on the channel of V polarization, and performing azimuth phase demodulation on the full polarimetric echo data using a preset rule for phase demodulation; and acquiring echo data by filtering out, using a preset azimuth filter, range-ambiguity energy from the phase demodulated full polarimetric echo data.

2. The method according to claim 1, wherein the performing, using a preset rule for phase modulation, azimuth transmission-phase modulation respectively on pulse signal data transmitted on a channel of Horizontal (H) polarization and pulse signal data transmitted on a channel of Vertical (V) polarization comprises: performing no phase modulation on the pulse signal data transmitted on the channel of H polarization; and performing linear phase modulation on the pulse signal data transmitted on the channel of V polarization by setting a product of and a sequence number of the pulse signal data being transmitted on the channel of V polarization as a modulation phase for performing linear phase modulation on the pulse signal data being transmitted on the channel of V polarization.

3. The method according to claim 1, wherein the performing azimuth phase demodulation on the full polarimetric echo data using a preset rule for phase demodulation comprises: performing phase demodulation on the full polarimetric echo data using a transmission phase consistent with a main signal.

4. The method according to claim 1, wherein the filtering out, using a preset azimuth filter, range-ambiguity energy from the phase demodulated full polarimetric echo data comprises: converting the phase demodulated full polarimetric echo data into phase demodulated full polarimetric echo data in a Doppler frequency domain; and filtering the phase demodulated full polarimetric echo data in the Doppler frequency domain using an azimuth frequency-domain Wiener filter.

5. The method according to claim 4, wherein the azimuth frequency-domain Wiener filter is associated with a formula of: $H\left(f\right)=\frac{S_0\left(f\right)}{S_0\left(f\right)+S_{-1}\left(f\right)+S_{+1}\left(f\right)+{}_t},$ wherein the S.sub.0(f) is an azimuth spectrum of a payload signal, each of the S.sub.1(f) and the S.sub.+1 (f) is an azimuth spectrum of a first-order ambiguous signal, and the .sub.t is a noise power spectrum.

6. The method according to claim 2, wherein the filtering out, using a preset azimuth filter, range-ambiguity energy from the phase demodulated full polarimetric echo data comprises: converting the phase demodulated full polarimetric echo data into phase demodulated full polarimetric echo data in a Doppler frequency domain; and filtering the phase demodulated full polarimetric echo data in the Doppler frequency domain using an azimuth frequency-domain Wiener filter.

7. The method according to claim 6, wherein the azimuth frequency-domain Wiener filter is associated with a formula of: $H\left(f\right)=\frac{S_0\left(f\right)}{S_0\left(f\right)+S_{-1}\left(f\right)+S_{+1}\left(f\right)+{}_t},$ wherein the S.sub.0(f) is an azimuth spectrum of a payload signal, each of the S.sub.1(f) and the S.sub.+1(f) is an azimuth spectrum of a first-order ambiguous signal, and the .sub.t is a noise power spectrum.

8. The method according to claim 3, wherein the filtering out, using a preset azimuth filter, range-ambiguity energy from the phase demodulated full polarimetric echo data comprises: converting the phase demodulated full polarimetric echo data into phase demodulated full polarimetric echo data in a Doppler frequency domain; and filtering the phase demodulated full polarimetric echo data in the Doppler frequency domain using an azimuth frequency-domain Wiener filter.

9. The method according to claim 8, wherein the azimuth frequency-domain Wiener filter is associated with a formula of: $H\left(f\right)=\frac{S_0\left(f\right)}{S_0\left(f\right)+S_{-1}\left(f\right)+S_{+1}\left(f\right)+{}_t},$ wherein the S.sub.0(f) is an azimuth spectrum of a payload signal, each of the S.sub.1(f) and the S.sub.+1(f) is an azimuth spectrum of a first-order ambiguous signal, and the .sub.t is a noise power spectrum.

10. A device for suppressing range ambiguity in Synthetic Aperture Radar (SAR), comprising: a processor; and memory storing instructions executable by the processor, wherein when executed by the processor, the instructions cause the processor to perform a method for suppressing range ambiguity in Synthetic Aperture Radar (SAR), the method comprising: performing, using a preset rule for phase modulation, azimuth transmission-phase modulation respectively on pulse signal data transmitted on a channel of Horizontal (H) polarization and pulse signal data transmitted on a channel of Vertical (V) polarization; acquiring full polarimetric echo data corresponding to the transmission-phase modulated pulse signal data transmitted on the channel of H polarization and the transmission-phase modulated pulse signal data transmitted on the channel of V polarization, and performing azimuth phase demodulation on the full polarimetric echo data using a preset rule for phase demodulation; and acquiring echo data by filtering out, using a preset azimuth filter, range-ambiguity energy from the phase demodulated full polarimetric echo data.

11. The device according to claim 10, wherein the performing, using a preset rule for phase modulation, azimuth transmission-phase modulation respectively on pulse signal data transmitted on a channel of Horizontal (H) polarization and pulse signal data transmitted on a channel of Vertical (V) polarization comprises: performing no phase modulation on the pulse signal data transmitted on the channel of H polarization; and performing linear phase modulation on the pulse signal data transmitted on the channel of V polarization by setting a product of and a sequence number of the pulse signal data being transmitted on the channel of V polarization as a modulation phase for performing linear phase modulation on the pulse signal data being transmitted on the channel of V polarization.

12. The device according to claim 10, wherein the performing azimuth phase demodulation on the full polarimetric echo data using a preset rule for phase demodulation comprises: performing phase demodulation on the full polarimetric echo data using a transmission phase consistent with a main signal.

13. The device according to claim 10, wherein the filtering out, using a preset azimuth filter, range-ambiguity energy from the phase demodulated full polarimetric echo data comprises: converting the phase demodulated full polarimetric echo data into phase demodulated full polarimetric echo data in a Doppler frequency domain; and filtering the phase demodulated full polarimetric echo data in the Doppler frequency domain using an azimuth frequency-domain Wiener filter.

14. The device according to claim 13, wherein the azimuth frequency-domain Wiener filter is associated with a formula of: $H\left(f\right)=\frac{S_0\left(f\right)}{S_0\left(f\right)+S_{-1}\left(f\right)+S_{+1}\left(f\right)+{}_t},$ wherein the S.sub.0(f) is an azimuth spectrum of a payload signal, each of the S.sub.1(f) and the S.sub.+1(f) is an azimuth spectrum of a first-order ambiguous signal, and the .sub.t is a noise power spectrum.

15. The device according to claim 11, wherein the filtering out, using a preset azimuth filter, range-ambiguity energy from the phase demodulated full polarimetric echo data comprises: converting the phase demodulated full polarimetric echo data into phase demodulated full polarimetric echo data in a Doppler frequency domain; and filtering the phase demodulated full polarimetric echo data in the Doppler frequency domain using an azimuth frequency-domain Wiener filter.

16. The device according to claim 15, wherein the azimuth frequency-domain Wiener filter is associated with a formula of: $H\left(f\right)=\frac{S_0\left(f\right)}{S_0\left(f\right)+S_{-1}\left(f\right)+S_{+1}\left(f\right)+{}_t},$ wherein the S.sub.0(f) is an azimuth spectrum of a payload signal, each of the S.sub.1(f) and the S.sub.+1(f) is an azimuth spectrum of a first-order ambiguous signal, and the .sub.t is a noise power spectrum.

17. The device according to claim 12, wherein the filtering out, using a preset azimuth filter, range-ambiguity energy from the phase demodulated full polarimetric echo data comprises: converting the phase demodulated full polarimetric echo data into phase demodulated full polarimetric echo data in a Doppler frequency domain; and filtering the phase demodulated full polarimetric echo data in the Doppler frequency domain using an azimuth frequency-domain Wiener filter.

18. The device according to claim 17, wherein the azimuth frequency-domain Wiener filter is associated with a formula of: $H\left(f\right)=\frac{S_0\left(f\right)}{S_0\left(f\right)+S_{-1}\left(f\right)+S_{+1}\left(f\right)+{}_t},$ wherein the S.sub.0(f) is an azimuth spectrum of a payload signal, each of the S.sub.1(f) and the S.sub.+1(f) is an azimuth spectrum of a first-order ambiguous signal, and the .sub.t is a noise power spectrum.

19. A computer-readable storage medium having stored therein instructions that, when executed by a processor, cause the processor to perform a method for suppressing range ambiguity in Synthetic Aperture Radar (SAR), the method comprising: performing, using a preset rule for phase modulation, azimuth transmission-phase modulation respectively on pulse signal data transmitted on a channel of Horizontal (H) polarization and pulse signal data transmitted on a channel of Vertical (V) polarization; acquiring full polarimetric echo data corresponding to the transmission-phase modulated pulse signal data transmitted on the channel of H polarization and the transmission-phase modulated pulse signal data transmitted on the channel of V polarization, and performing azimuth phase demodulation on the full polarimetric echo data using a preset rule for phase demodulation; and acquiring echo data by filtering out, using a preset azimuth filter, range-ambiguity energy from the phase demodulated full polarimetric echo data.

20. The storage medium according to claim 19, wherein the performing, using a preset rule for phase modulation, azimuth transmission-phase modulation respectively on pulse signal data transmitted on a channel of Horizontal (H) polarization and pulse signal data transmitted on a channel of Vertical (V) polarization comprises: performing no phase modulation on the pulse signal data transmitted on the channel of H polarization; and performing linear phase modulation on the pulse signal data transmitted on the channel of V polarization by setting a product of and a sequence number of the pulse signal data being transmitted on the channel of V polarization as a modulation phase for performing linear phase modulation on the pulse signal data being transmitted on the channel of V polarization.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0041] FIG. 1 is a flowchart of a method for suppressing range ambiguity in SAR according to an embodiment herein.

[0042] FIG. 2 is a diagram of a modulated transmission phase of full polarimetric SAR according to an embodiment herein.

[0043] FIG. 3 is a diagram of an azimuth spectrum of an azimuth phase demodulated echo of full polarimetric SAR according to an embodiment herein.

[0044] FIG. 4 is a diagram of Range-Ambiguity-to-Signal Ratios (RASR) of channels of co-polarization and of cross polarization according to an embodiment herein.

[0045] FIG. 5 is a diagram of results measured on channels of co-polarization and of cross polarization according to an embodiment herein.

[0046] FIG. 6 is a diagram of a structure of a device for suppressing range ambiguity in SAR according to an embodiment herein.

DETAILED DESCRIPTION

[0047] According to embodiments herein, azimuth transmission-phase modulation is performed respectively, using a preset rule for phase modulation, on pulse signal data transmitted on a channel of Horizontal (H) polarization and pulse signal data transmitted on a channel of Vertical (V) polarization. Full polarimetric echo data corresponding to the transmission-phase modulated pulse signal data transmitted on the channel of H polarization and the transmission-phase modulated pulse signal data transmitted on the channel of V polarization are acquired. Azimuth phase demodulation is performed on the full polarimetric echo data using a preset rule for phase demodulation. Echo data are acquired by filtering out, using a preset azimuth filter, range-ambiguity energy from the phase demodulated full polarimetric echo data.

[0048] The present disclosure will be elaborated below with reference to an embodiment herein.

[0049] As shown in FIG. 1, a method for suppressing range ambiguity in SAR according to an embodiment herein includes step(s) as follows.

[0050] In step 101, azimuth transmission-phase modulation is performed respectively, using a preset rule for phase modulation, on pulse signal data transmitted on a channel of H polarization and pulse signal data transmitted on a channel of V polarization.

[0051] In full polarimetric SAR, a pulse may be transmitted alternatively on a channel of H polarization (H channel) and a channel of V polarization (V channel). Azimuth transmission-phase modulation may be performed on the pulse signal data corresponding to the pulse transmitted alternatively on the channel of H polarization and the channel of V polarization. The full polarimetric SAR may be spaceborne.

[0052] As shown in FIG. 2, the preset rule for phase modulation may be as follows. No phase modulation may be performed on the pulse signal data transmitted on the channel of H polarization. Linear phase modulation may be performed on the pulse signal data transmitted on the channel of V polarization, such as by setting a product of and a sequence number of the pulse signal data being transmitted on the channel of V polarization as a modulation phase for performing linear phase modulation on the pulse signal data being transmitted on the channel of V polarization, as shown by a formula (1).

.sub.mod.sub._.sub.H=0

.sub.mod.sub._.sub.V=n(1)

[0053] The .sub.mod.sub._.sub.H may be a modulation phase of the pulse signal transmitted on the channel of H polarization. The .sub.mod.sub._.sub.V may be the modulation phase of the pulse signal transmitted on the channel of V polarization. The n may be the sequence number of the pulse signal data that specifies an nth pulse transmitted.

[0054] The modulated pulse signal data may then be transmitted. A polarized signal H+V.Math.exp(j) and a polarized signal HV.Math.exp(j) may be transmitted, alternatively.

[0055] In step 102, full polarimetric echo data corresponding to the transmission-phase modulated pulse signal data transmitted on the channel of H polarization and the transmission-phase modulated pulse signal data transmitted on the channel of V polarization are acquired. Azimuth phase demodulation is performed on the full polarimetric echo data using a preset rule for phase demodulation.

[0056] In full polarimetric SAR, a pulse may be transmitted respectively on the channel of H polarization and the channel of V polarization. The full polarimetric echo data may then be formed by polarimetric echo data of four channels HH, HV, VH, and VV received respectively on the channel of H polarization and the channel of V polarization.

[0057] Phase demodulation may be performed, using the preset rule for phase demodulation, on the overall full polarimetric echo signal of the four-channel polarimetric echo data received on the channel of H polarization and the channel of V polarization of full polarimetric SAR. According to the preset rule for phase demodulation, phase demodulation may be performed on the overall full polarimetric echo signal using a transmission phase consistent with a main signal, as shown by a formula (2). The main signal may refer to a payload (i.e., useful) signal of full polarimetric SAR.

.sub.demo(n)=.sub.mod(nm)(2)

[0058] The .sub.demo(n) may be a demodulation phase. The n may specify nth pulse received. The m may account for a round-trip delay between a pulse transmitted and a pulse received.

[0059] The modulation and demodulation may have no impact on the main signal. The demodulation phase used is the transmission phase used for the main signal. Therefore, a residual phase of the main signal after the demodulation is 0. Accordingly, there will be no shift in an azimuth spectrum of the main signal.

[0060] A residual phase of a kth-order range ambiguous signal in the four-channel polarimetric echo data of full polarimetric SAR having gone through the modulation and the demodulation may be expressed by a formula (3).

[00003]$\begin{array}{cc}{}_{{\mathrm{res}}_{}.Math.p}\left(n,k\right)=\{\begin{array}{cc}0.Math.& p=\mathrm{HH}.Math.\text{/}.Math.\mathrm{VH},k=\mathrm{even}\\ \frac{k}{2}.Math..Math.& p=\mathrm{HV}.Math.\text{/}.Math.\mathrm{VV},k=\mathrm{even}\\ n.Math..Math.+\frac{k-1}{2}.Math.& k=\mathrm{odd}.Math.\end{array}& \left(3\right)\end{array}$

[0061] The p may be one of the four channels of polarization HH, HV, VH, and VV. An even k may indicate an even ambiguous region. An odd k may indicate an odd ambiguous region.

[0062] According to the formula (3), after the modulation and the demodulation of full polarimetric SAR, there is no residual phase for the channels HH and VH in an even ambiguous region, and a residual phase of

[00004]$\frac{k}{2}.Math.$

for the channels HV and VV in an even ambiguous region. Therefore, a fixed residual phase due to the modulation and the demodulation may have no impact on a spectrum of a signal in an even ambiguous region. On the other hand, an azimuth linear residual phase

[00005]$n.Math..Math.+\frac{k-1}{2}.Math.$

may be produced by odd-order signals of all components of the H channel and the V channel. A formula (4) may be acquired according to an azimuth time t.sub.m=n/PRF.

[00006]$\begin{array}{cc}\begin{array}{c}f\left({}_{{\mathrm{res}}_{}.Math.\mathrm{odd}}\right)=\exp .Math.\left\{\mathrm{jn}.Math..Math.+j.Math.\frac{k-1}{2}.Math.\right\}\\ =\exp .Math.\left\{j.Math.\mathrm{PRF}.Math.t_m.Math.+j.Math.\frac{k-1}{2}.Math.\right\}\\ =\exp .Math.\left\{j.Math..Math.2.Math.\left(\frac{\mathrm{PRF}}{2}.Math.t_m+\frac{k-1}{4}\right)\right\}\end{array}& \left(4\right)\end{array}$

[0063] According to the formula (4), a frequency shift of PRF/2 in a Doppler frequency domain may be produced for a signal in an odd ambiguous region. FIG. 3 is a diagram of an azimuth spectrum of the phase demodulated echo data.

[0064] In step 103, echo data are acquired by filtering out, using a preset azimuth filter, range-ambiguity energy from the phase demodulated full polarimetric echo data.

[0065] An azimuth filter may be preset according to the rule for phase modulation and the rule for phase demodulation. The phase demodulated full polarimetric echo data may then be filtered using the preset azimuth filter. The preset azimuth filter may be an azimuth frequency-domain Wiener filter. The azimuth frequency-domain Wiener filter may be constructed based on a Minimum-Mean-Square-Error (MMSE) theory, i.e., by minimizing a mean square error of an unambiguous signal and a filtered signal. The phase demodulated full polarimetric echo data may be converted into phase demodulated full polarimetric echo data in the Doppler frequency domain. The phase demodulated full polarimetric echo data in the Doppler frequency domain may then be filtered using the azimuth frequency-domain Wiener filter. Out-of-band ambiguity energy may thus be removed. Part of in-band ambiguity energy may be removed as well.

[0066] A key reason the azimuth frequency-domain Wiener filter may be used is that the spectrum of an ambiguous signal is inconsistent with that of the main signal. The modulation and the demodulation may convert the range-ambiguity energy into azimuth-ambiguity energy. A range ambiguous spectrum in the Doppler frequency domain may be shifted with respect to the spectrum of the main signal. A Wiener filter may show high performance in solving an MMSE-related problem. The Wiener filter may be expressed by a formula (5).

[00007]$\begin{array}{cc}H\left(f\right)=\frac{S_0\left(f\right)}{S_0\left(f\right)+S_{-1}\left(f\right)+S_{+1}\left(f\right)+{}_t}& \left(5\right)\end{array}$

[0067] The S.sub.0(f) may be the azimuth spectrum of the main signal, namely, the payload signal. The S.sub.1(f) may be an azimuth spectrum of an ambiguous signal of an order 1. The S.sub.+1(f) may be an azimuth spectrum of an ambiguous signal of an order 1. The .sub.t may be a noise power spectrum. Each of the S.sub.0(f), the S.sub.1(f), the S.sub.+1(f), and the .sub.t may be acquired using an existing computational model according to a parameter of full polarimetric SAR.

[0068] With the step, the range ambiguity in full polarimetric SAR may be suppressed, and the echo data with the range-ambiguity energy removed may be acquired.

[0069] The present disclosure is elaborated below with reference to one or more examples.

[0070] In an example, a Range-Ambiguity-To-Signal Ratio (RASR) of full polarimetric SAR may be simulated and compared. Table 1 shows a parameter of radar in a band L.

TABLE-US-00001 TABLE 1 carrier frequency of 1.26 GHz bandwidth of 140 MHz transmitted signal transmitted signal width of transmitted signal 70 us pitch angle 35 height of platform 607 km speed of flight of 7,500 m/s platform system Doppler bandwidth 1,235.9 Hz length of antenna 2.9 m

[0071] After the azimuth filtering according to an embodiment herein, the RASR of full polarimetric SAR may be expressed by a formula (6).

[00008]$\begin{array}{cc}\mathrm{RASR}=\frac{\underset{k0}{.Math.}.Math.C_{\mathrm{ak}}.Math.S_{\mathrm{ak}}}{C_0.Math.S_0}& \left(6\right)\end{array}$

[0072] The C.sub.ak may be a coefficient for filtering an azimuth spectrum of a kth-order ambiguous signal. The C.sub.0 may be an azimuth filtering ratio of the payload signal. The S.sub.ak may be energy of the kth-order ambiguous signal. The S.sub.0 may be energy of the payload signal.

[0073] Using the Wiener filter, the coefficient C.sub.ak and the ratio C.sub.0 may be computed respectively using a formula (7) and a formula (8).

[00009]$\begin{array}{cc}{C}_{\mathrm{ak}}=\frac{{}_{-B_a.Math.\text{/}.Math.2}^{B_a.Math.\text{/}.Math.2}.Math.H^2\left(f\right).Math.G_{\mathrm{ak}}^2\left(f+\frac{k.Math.\mathrm{PRF}}{2}\right).Math.\mathrm{df}}{{}_{-\mathrm{PRF}.Math.\text{/}.Math.2}^{\mathrm{PRF}.Math.\text{/}.Math.2}.Math.G_{\mathrm{ak}}^2\left(f+\frac{k.Math.\mathrm{PRF}}{2}\right).Math.\mathrm{df}}& \left(7\right)\\ C_0=\frac{{}_{-B_a.Math.\text{/}.Math.2}^{B_a.Math.\text{/}.Math.2}.Math.H^2\left(f\right).Math.G_0^2\left(f\right).Math.\mathrm{df}}{{}_{-\mathrm{PRF}.Math.\text{/}.Math.2}^{\mathrm{PRF}.Math.\text{/}.Math.2}.Math.G_{\mathrm{ak}}^2\left(f\right).Math.\mathrm{df}}& \left(8\right)\end{array}$

[0074] The H(f) may be the frequency-domain Wiener filter. The G.sub.ak may be an azimuth radiation pattern of the system antenna. The G.sub.0 may be a pattern of the payload signal.

[0075] In FIG. 4, panel (a) panel (b) show the RASR of a channel of co-polarization and the RASR of a channel of cross polarization, respectively, when applying the present disclosure to conventional full polarimetric SAR in the band L. As simulated, data on the channel of co-polarization of full polarimetric SAR in an odd ambiguous region may be of an ambiguous signal of cross polarization, so that ambiguity performance of the channel of co-polarization is not degraded and may meet a system requirement. However, data on the channel of cross polarization in an odd ambiguous region may be of an ambiguous signal of co-polarization. The ambiguous signal thus may have a major impact on the payload signal. The RASR of the channel of cross polarization may drop dramatically. The channel of co-polarization may include a low-energy signal of HV polarization in an odd ambiguous region. Therefore, with the method according to an embodiment herein, the RASR of the channel of co-polarization may be improved but marginally, by about 2 dB. The channel of cross polarization may include a high-energy signal of co-polarization in an odd ambiguous region. Therefore, with the method according to an embodiment herein, more energy may be filtered out compared with the case of the channel of co-polarization. Accordingly, improvement to the range ambiguity performance of the channel of cross polarization is significant, compared to improvement to that of the channel of co-polarization, with a RASR of the channel of cross polarization improved by about 7 dB.

[0076] FIG. 5 shows images measured by SAR on the channel of VH polarization and the channel of HH polarization in the band L when applying and not applying the present disclosure. FIG. 5 shows significant suppression of range ambiguity of an area object by the method according to an embodiment herein. In FIG. 5, panel (a) is an image of echo of VH cross polarization by conventional full polarimetric SAR. Significant contamination by range ambiguity of co-polarization may sharply lower quality of the image. Compared with FIG. 5, panel (a) of VH cross polarization, there is barely any impact of range ambiguity on an image of echo of HH co-polarization by conventional full polarimetric SAR in FIG. 5, panel (b). In FIG. 5, panel (c) and panel (d) are an image of echo of VH cross polarization and an image of echo of HH co-polarization, respectively, with range ambiguity suppressed using the method according to an embodiment herein. Most ambiguity energy on the channels HH and HV of conventional full polarimetric SAR may be removed by applying the range ambiguity suppression according to an embodiment herein, greatly improving performance of images of the channels VH and HH.

[0077] As shown in FIG. 6, a device for suppressing range ambiguity in SAR according to an embodiment herein includes a modulating module 61, a demodulating module 62, and a filtering module 63.

[0078] The modulating module 61 is arranged for: performing, using a preset rule for phase modulation, azimuth transmission-phase modulation respectively on pulse signal data transmitted on a channel of Horizontal (H) polarization and pulse signal data transmitted on a channel of Vertical (V) polarization.

[0079] In full polarimetric SAR, a pulse may be transmitted alternatively on a channel of H polarization (H channel) and a channel of V polarization (V channel). Azimuth transmission-phase modulation may be performed on the pulse signal data corresponding to the pulse transmitted alternatively on the channel of H polarization and the channel of V polarization. The full polarimetric SAR may be spaceborne.

[0080] As shown in FIG. 2, the preset rule for phase modulation may be as follows. No phase modulation may be performed on the pulse signal data transmitted on the channel of H polarization. Linear phase modulation may be performed on the pulse signal data transmitted on the channel of V polarization, such as by setting a product of and a sequence number of the pulse signal data being transmitted on the channel of V polarization as a modulation phase for performing linear phase modulation on the pulse signal data being transmitted on the channel of V polarization, as shown by the formula (1).

[0081] The .sub.mod.sub._.sub.H may be a modulation phase of the pulse signal transmitted on the channel of H polarization. The .sub.mod.sub._.sub.V may be the modulation phase of the pulse signal transmitted on the channel of V polarization. The n may be the sequence number of the pulse signal data that specifies an nth pulse transmitted.

[0082] The modulated pulse signal data may then be transmitted. A polarized signal H+V.Math.exp(j) and a polarized signal HV.Math.exp(j) may be transmitted, alternatively.

[0083] The demodulating module 62 is arranged for: acquiring full polarimetric echo data corresponding to the transmission-phase modulated pulse signal data transmitted on the channel of H polarization and the transmission-phase modulated pulse signal data transmitted on the channel of V polarization, and performing azimuth phase demodulation on the full polarimetric echo data using a preset rule for phase demodulation.

[0084] In full polarimetric SAR, a pulse may be transmitted respectively on the channel of H polarization and the channel of V polarization. The full polarimetric echo data may then be formed by polarimetric echo data of four channels HH, HV, VH, and VV received respectively on the channel of H polarization and the channel of V polarization.

[0085] Phase demodulation may be performed, using the preset rule for phase demodulation, on the overall full polarimetric echo signal of the four-channel polarimetric echo data received on the channel of H polarization and the channel of V polarization of full polarimetric SAR. According to the preset rule for phase demodulation, phase demodulation may be performed on the overall full polarimetric echo signal using a transmission phase consistent with a main signal, as shown by the formula (2). The main signal may refer to a payload (i.e., useful) signal of full polarimetric SAR.

[0086] The .sub.demo(n), an may be a demodulation phase. The n may specify nth pulse received. The m may account for a round-trip delay between a pulse transmitted and a pulse received.

[0087] The modulation and demodulation may have no impact on the main signal. The demodulation phase used is the transmission phase used for the main signal. Therefore, a residual phase of the main signal after the demodulation is 0. Accordingly, there will be no shift in an azimuth spectrum of the main signal.

[0088] A residual phase of a kth-order range ambiguous signal in the four-channel polarimetric echo data of full polarimetric SAR having gone through the modulation and the demodulation may be expressed by the formula (3).

[0089] The p may be one of the four channels of polarization HH, HV, VH, and VV. An even k may indicate an even ambiguous region. An odd k may indicate an odd ambiguous region.

[0090] According to the formula (3), after the modulation and the demodulation of full polarimetric SAR, there is no residual phase for the channels HH and VH in an even ambiguous region, and a residual phase of

[00010]$\frac{k}{2}.Math.$

for the channels HV and VV in an even ambiguous region. Therefore, a fixed residual phase due to the modulation and the demodulation may have no impact on a spectrum of a signal in an even ambiguous region. On the other hand, an azimuth linear residual phase

[00011]$n.Math..Math.+\frac{k-1}{2}.Math.$

may be produced by odd-order signals of all components of the H channel and the V channel. The formula (4) may be acquired according to an azimuth time t.sub.m=n/PRF.

[0091] According to the formula (4), a frequency shift of PRF/2 in a Doppler frequency domain may be produced for a signal in an odd ambiguous region. FIG. 3 is a diagram of an azimuth spectrum of the phase demodulated echo data.

[0092] The filtering module 63 is arranged for: acquiring echo data by filtering out, using a preset azimuth filter, range-ambiguity energy from the phase demodulated full polarimetric echo data.

[0093] An azimuth filter may be preset according to the rule for phase modulation and the rule for phase demodulation. The phase demodulated full polarimetric echo data may then be filtered using the preset azimuth filter. The preset azimuth filter may be an azimuth frequency-domain Wiener filter. The azimuth frequency-domain Wiener filter may be constructed based on a Minimum-Mean-Square-Error (MMSE) theory, i.e., by minimizing a mean square error of an unambiguous signal and a filtered signal. The phase demodulated full polarimetric echo data may be converted into phase demodulated full polarimetric echo data in the Doppler frequency domain. The phase demodulated full polarimetric echo data in the Doppler frequency domain may then be filtered using the azimuth frequency-domain Wiener filter. Out-of-band ambiguity energy may thus be removed. Part of in-band ambiguity energy may be removed as well.

[0094] A key reason the azimuth frequency-domain Wiener filter may be used is that the spectrum of an ambiguous signal is inconsistent with that of the main signal. The modulation and the demodulation may convert the range-ambiguity energy into azimuth-ambiguity energy. A range ambiguous spectrum in the Doppler frequency domain may be shifted with respect to the spectrum of the main signal. A Wiener filter may show high performance in solving an MMSE-related problem. The Wiener filter may be expressed by the formula (5).

[0095] The S.sub.0(f) may be the azimuth spectrum of the main signal, namely, the payload signal. The S.sub.1(f) may be an azimuth spectrum of an ambiguous signal of an order 1. The S.sub.+1(f) may be an azimuth spectrum of an ambiguous signal of an order 1. The .sub.t may be a noise power spectrum. Each of the S.sup.0 (f), the S.sub.1(f), the S.sub.+1(f), and the .sub.t may be acquired using an existing computational model according to a parameter of full polarimetric SAR.

[0096] With the step, the range ambiguity in full polarimetric SAR may be suppressed, and the echo data with the range-ambiguity energy removed may be acquired.

[0097] Each of the modulating module 61, the demodulating module 62, and the filtering module 63 may be implemented by at least one of a Central Processing Unit (CPU), a Micro Processing Unit (MPU), a Digital Signal Processor (DSP), a Field-Programmable Gate Array (FPGA), etc., in SAR.

[0098] Those skilled in the art may understand that all or part of the steps of the method may be completed by instructing a related hardware through a program. The program may be stored in a transitory or non-transitory computer-readable storage medium. The computer-readable storage medium may be Read-Only Memory (ROM), a magnetic disk, a CD, and/or the like. All or part of the steps in an embodiment herein may be implemented using one or more Integrated Circuits. A unit or module in an embodiment herein may be implemented in form of hardware or a software functional module. The present disclosure is not limited to any specific combination of hardware and software.

[0099] According to an embodiment herein, a computer-readable storage medium has stored therein instructions that, when executed by a processor, cause the processor to perform a method for suppressing range ambiguity in SAR. As shown in FIG. 1, the method includes step(s) as follows.

[0100] In step 101, azimuth transmission-phase modulation is performed respectively, using a preset rule for phase modulation, on pulse signal data transmitted on a channel of H polarization and pulse signal data transmitted on a channel of V polarization.

[0101] In full polarimetric SAR, a pulse may be transmitted alternatively on a channel of H polarization (H channel) and a channel of V polarization (V channel). Azimuth transmission-phase modulation may be performed on the pulse signal data corresponding to the pulse transmitted alternatively on the channel of H polarization and the channel of V polarization. The full polarimetric SAR may be spaceborne.

[0102] As shown in FIG. 2, the preset rule for phase modulation may be as follows. No phase modulation may be performed on the pulse signal data transmitted on the channel of H polarization. Linear phase modulation may be performed on the pulse signal data transmitted on the channel of V polarization, such as by setting a product of 7E and a sequence number of the pulse signal data being transmitted on the channel of V polarization as a modulation phase for performing linear phase modulation on the pulse signal data being transmitted on the channel of V polarization, as shown by the formula (1).

[0103] The .sub.mod.sub._.sub.H may be a modulation phase of the pulse signal transmitted on the channel of H polarization. The .sub.mod.sub._.sub.V may be the modulation phase of the pulse signal transmitted on the channel of V polarization. The n may be the sequence number of the pulse signal data that specifies an nth pulse transmitted.

[0104] The modulated pulse signal data may then be transmitted. A polarized signal H+V.Math.exp(j) and a polarized signal HV.Math.exp(j) may be transmitted, alternatively.

[0105] In step 102, full polarimetric echo data corresponding to the transmission-phase modulated pulse signal data transmitted on the channel of H polarization and the transmission-phase modulated pulse signal data transmitted on the channel of V polarization are acquired. Azimuth phase demodulation is performed on the full polarimetric echo data using a preset rule for phase demodulation.

[0106] In full polarimetric SAR, a pulse may be transmitted respectively on the channel of H polarization and the channel of V polarization. The full polarimetric echo data may then be formed by polarimetric echo data of four channels HH, HV, VH, and VV received respectively on the channel of H polarization and the channel of V polarization.

[0107] Phase demodulation may be performed, using the preset rule for phase demodulation, on the overall full polarimetric echo signal of the four-channel polarimetric echo data received on the channel of H polarization and the channel of V polarization of full polarimetric SAR. According to the preset rule for phase demodulation, phase demodulation may be performed on the overall full polarimetric echo signal using a transmission phase consistent with a main signal, as shown by the formula (2). The main signal may refer to a payload (i.e., useful) signal of full polarimetric SAR.

[0108] The .sub.demo(n) may be a demodulation phase. The n may specify an nth pulse received. The m may account for a round-trip delay between a pulse transmitted and a pulse received.

[0109] The modulation and demodulation may have no impact on the main signal. The demodulation phase used is the transmission phase used for the main signal. Therefore, a residual phase of the main signal after the demodulation is 0. Accordingly, there will be no shift in an azimuth spectrum of the main signal.

[0110] A residual phase of a kth-order range ambiguous signal in the four-channel polarimetric echo data of full polarimetric SAR having gone through the modulation and the demodulation may be expressed by the formula (3).

[0111] The p may be one of the four channels of polarization HH, HV, VH, and VV. An even k may indicate an even ambiguous region. An odd k may indicate an odd ambiguous region.

[0112] According to the formula (3), after the modulation and the demodulation of full polarimetric SAR, there is no residual phase for the channels HH and VH in an even ambiguous region, and a residual phase of

[00012]$\frac{k}{2}.Math.$

for the channels HV and VV in an even ambiguous region. Therefore, a fixed residual phase due to the modulation and the demodulation may have no impact on a spectrum of a signal in an even ambiguous region. On the other hand, an azimuth linear residual phase

[00013]$n.Math..Math.+\frac{k-1}{2}.Math.$

may be produced by odd-order signals of all components of the H channel and the V channel. The formula (4) may be acquired according to an azimuth time t.sub.m=n/PRF.

[0113] According to the formula (4), a frequency shift of PRF/2 in a Doppler frequency domain may be produced for a signal in an odd ambiguous region. FIG. 3 is a diagram of an azimuth spectrum of the phase demodulated echo data.

[0114] In step 103, echo data are acquired by filtering out, using a preset azimuth filter, range-ambiguity energy from the phase demodulated full polarimetric echo data.

[0115] An azimuth filter may be preset according to the rule for phase modulation and the rule for phase demodulation. The phase demodulated full polarimetric echo data may then be filtered using the preset azimuth filter. The preset azimuth filter may be an azimuth frequency-domain Wiener filter. The azimuth frequency-domain Wiener filter may be constructed based on a Minimum-Mean-Square-Error (MMSE) theory, i.e., by minimizing a mean square error of an unambiguous signal and a filtered signal. The phase demodulated full polarimetric echo data may be converted into phase demodulated full polarimetric echo data in the Doppler frequency domain. The phase demodulated full polarimetric echo data in the Doppler frequency domain may then be filtered using the azimuth frequency-domain Wiener filter. Out-of-band ambiguity energy may thus be removed. Part of in-band ambiguity energy may be removed as well.

[0116] A key reason the azimuth frequency-domain Wiener filter may be used is that the spectrum of an ambiguous signal is inconsistent with that of the main signal. The modulation and the demodulation may convert the range-ambiguity energy into azimuth-ambiguity energy. A range ambiguous spectrum in the Doppler frequency domain may be shifted with respect to the spectrum of the main signal. A Wiener filter may show high performance in solving an MMSE-related problem. The Wiener filter may be expressed by the formula (5).

[0117] The S.sub.0(f) may be the azimuth spectrum of the main signal, namely, the payload signal. The S.sub.1(f) may be an azimuth spectrum of an ambiguous signal of an order 1. The S.sub.+1(f) may be an azimuth spectrum of an ambiguous signal of an order 1. The .sub.t may be a noise power spectrum. Each of the S.sub.0(f), the S.sub.1(f), the S.sub.+1(f), and the .sub.t may be acquired using an existing computational model according to a parameter of full polarimetric SAR.

[0118] With the step, the range ambiguity in full polarimetric SAR may be suppressed, and the echo data with the range-ambiguity energy removed may be acquired.

[0119] According to an embodiment herein, a device for suppressing range ambiguity in SAR includes a processor and memory storing instructions executable by the processor. When executed by the processor, the instructions cause the processor to perform a method for suppressing range ambiguity in SAR. As shown in FIG. 1, the method includes step(s) as follows.

[0120] In step 101, azimuth transmission-phase modulation is performed respectively, using a preset rule for phase modulation, on pulse signal data transmitted on a channel of H polarization and pulse signal data transmitted on a channel of V polarization.

[0121] In full polarimetric SAR, a pulse may be transmitted alternatively on a channel of H polarization (H channel) and a channel of V polarization (V channel). Azimuth transmission-phase modulation may be performed on the pulse signal data corresponding to the pulse transmitted alternatively on the channel of H polarization and the channel of V polarization. The full polarimetric SAR may be spaceborne.

[0122] As shown in FIG. 2, the preset rule for phase modulation may be as follows. No phase modulation may be performed on the pulse signal data transmitted on the channel of H polarization. Linear phase modulation may be performed on the pulse signal data transmitted on the channel of V polarization, such as by setting a product of and a sequence number of the pulse signal data being transmitted on the channel of V polarization as a modulation phase for performing linear phase modulation on the pulse signal data being transmitted on the channel of V polarization, as shown by the formula (1).

[0123] The .sub.mod.sub._.sub.H may be a modulation phase of the pulse signal transmitted on the channel of H polarization. The .sub.mod.sub._.sub.V may be the modulation phase of the pulse signal transmitted on the channel of V polarization. The n may be the sequence number of the pulse signal data that specifies an nth pulse transmitted.

[0124] The modulated pulse signal data may then be transmitted. A polarized signal H+V.Math.exp(j) and a polarized signal HV.Math.exp(j) may be transmitted, alternatively.

[0125] In step 102, full polarimetric echo data corresponding to the transmission-phase modulated pulse signal data transmitted on the channel of H polarization and the transmission-phase modulated pulse signal data transmitted on the channel of V polarization are acquired. Azimuth phase demodulation is performed on the full polarimetric echo data using a preset rule for phase demodulation.

[0126] In full polarimetric SAR, a pulse may be transmitted respectively on the channel of H polarization and the channel of V polarization. The full polarimetric echo data may then be formed by polarimetric echo data of four channels HH, HV, VH, and VV received respectively on the channel of H polarization and the channel of V polarization.

[0127] Phase demodulation may be performed, using the preset rule for phase demodulation, on the overall full polarimetric echo signal of the four-channel polarimetric echo data received on the channel of H polarization and the channel of V polarization of full polarimetric SAR. According to the preset rule for phase demodulation, phase demodulation may be performed on the overall full polarimetric echo signal using a transmission phase consistent with a main signal, as shown by the formula (2). The main signal may refer to a payload (i.e., useful) signal of full polarimetric SAR.

[0128] The .sub.demo(n) may be a demodulation phase. The n may specify an nth pulse received. The m may account for a round-trip delay between a pulse transmitted and a pulse received.

[0129] The modulation and demodulation may have no impact on the main signal. The demodulation phase used is the transmission phase used for the main signal. Therefore, a residual phase of the main signal after the demodulation is 0. Accordingly, there will be no shift in an azimuth spectrum of the main signal.

[0130] A residual phase of a kth-order range ambiguous signal in the four-channel polarimetric echo data of full polarimetric SAR having gone through the modulation and the demodulation may be expressed by the formula (3).

[0131] The p may be one of the four channels of polarization HH, HV, VH, and VV. An even k may indicate an even ambiguous region. An odd k may indicate an odd ambiguous region.

[0132] According to the formula (3), after the modulation and the demodulation of full polarimetric SAR, there is no residual phase for the channels HH and VH in an even ambiguous region, and a residual phase of

[00014]$\frac{k}{2}.Math.$

for the channels HV and VV in an even ambiguous region. Therefore, a fixed residual phase due to the modulation and the demodulation may have no impact on a spectrum of a signal in an even ambiguous region. On the other hand, an azimuth linear residual phase

[00015]$n.Math..Math.+\frac{k-1}{2}.Math.$

may be produced by odd-order signals of all components of the H channel and the V channel. The formula (4) may be acquired according to an azimuth time t.sub.m=n/PRF.

[0133] According to the formula (4), a frequency shift of PRF/2 in a Doppler frequency domain may be produced for a signal in an odd ambiguous region. FIG. 3 is a diagram of an azimuth spectrum of the phase demodulated echo data.

[0134] In step 103, echo data are acquired by filtering out, using a preset azimuth filter, range-ambiguity energy from the phase demodulated full polarimetric echo data.

[0135] An azimuth filter may be preset according to the rule for phase modulation and the rule for phase demodulation. The phase demodulated full polarimetric echo data may then be filtered using the preset azimuth filter. The preset azimuth filter may be an azimuth frequency-domain Wiener filter. The azimuth frequency-domain Wiener filter may be constructed based on a Minimum-Mean-Square-Error (MMSE) theory, i.e., by minimizing a mean square error of an unambiguous signal and a filtered signal. The phase demodulated full polarimetric echo data may be converted into phase demodulated full polarimetric echo data in the Doppler frequency domain. The phase demodulated full polarimetric echo data in the Doppler frequency domain may then be filtered using the azimuth frequency-domain Wiener filter. Out-of-band ambiguity energy may thus be removed. Part of in-band ambiguity energy may be removed as well.

[0136] A key reason the azimuth frequency-domain Wiener filter may be used is that the spectrum of an ambiguous signal is inconsistent with that of the main signal. The modulation and the demodulation may convert the range-ambiguity energy into azimuth-ambiguity energy. A range ambiguous spectrum in the Doppler frequency domain may be shifted with respect to the spectrum of the main signal. A Wiener filter may show high performance in solving an MMSE-related problem. The Wiener filter may be expressed by the formula (5).

[0137] The S.sub.0(f) may be the azimuth spectrum of the main signal, namely, the payload signal. The S.sub.1(f) may be an azimuth spectrum of an ambiguous signal of an order 1. The S.sub.+1(f) may be an azimuth spectrum of an ambiguous signal of an order 1. The .sub.t may be a noise power spectrum. Each of the S.sub.0(f), the S.sub.1(f), the S.sub.+1(f), and the .sub.t may be acquired using an existing computational model according to a parameter of full polarimetric SAR.

[0138] With the step, the range ambiguity in full polarimetric SAR may be suppressed, and the echo data with the range-ambiguity energy removed may be acquired.

[0139] What described are merely embodiments of the present disclosure, and are not intended to limit the scope of the present disclosure. Any modification, equivalent replacement, improvement, and/or the like made within the spirit and principle of the present disclosure should be included in the scope of the present disclosure.