Method and apparatus for computer-assisted processing of SAR raw data

Abstract

The present disclosure relates to a method for computer-implemented processing of SAR raw data, which comprises radar echoes from radar pulse. An interference radar echo and an interference pulse are associated with a respective radar pulse, wherein the interference pulse and the respective radar pulse have orthogonal waveforms. SAR raw data are focused by a first focusing on the interference pulses including a range compression and an azimuth compression, to obtain first focused data, where a filtering is used for the range compression of a respective radar pulse, which filtering is matched to the waveform of the associated interference pulse. Thereafter, the first focused data undergo a signal suppression, which at least partially suppresses the interference radar echo, as a result of which second focused data are obtained. These second focused data finally undergo a defocusing including range decompression and azimuth decompression to obtain modified SAR raw data.

Claims

1. A method for computer-implemented processing of Synthetic Aperture Radar (SAR) raw data, which comprise radar echoes from radar pulses, which have been emitted by a radar device on at least one flying object which is moving in an azimuth direction over the earth's surface, wherein the radar echoes of the radar pulses have been received and recorded by the radar device and wherein the SAR raw data contain data samples for a plurality of range positions for each radar pulse, wherein the range positions represent different distances between the earth's surface and the azimuth direction, wherein an interference radar echo and an interference pulse are associated with a respective radar pulse, wherein the interference radar echo is a radar echo of the interference pulse for a first range position and the interference radar echo is received by the radar device simultaneously with a radar echo of the respective radar pulse for a second range position, wherein the interference pulse is a radar pulse other than the respective radar pulse and wherein the interference pulse and the respective radar pulse have orthogonal waveforms, wherein a) the SAR raw data are focused by a first focusing on the interference pulses associated with the respective radar pulses, as a result of which first focused data are obtained, wherein a range compression and an azimuth compression are carried out as part of the first focusing where a filtering is used for the range compression of a respective radar pulse, which filtering is matched to the waveform of the associated interference pulse, b) the first focused data undergo a signal suppression, which at least partially suppresses the associated interference radar echo for each radar pulse, as a result of which second focused data are obtained; c) the second focused data undergo defocusing, as a result of which modified SAR raw data are obtained, wherein as part of the defocusing a range decompression and an azimuth decompression are carried out, which are inverse to the range compression and the azimuth compression of the first focusing.

2. The method according to claim 1, wherein the modified raw data are focused on the respective radar pulses by a second focusing, as a result of which third focused data are obtained, wherein a range compression and an azimuth compression are carried out as part of the second focusing, where a filtering is used for the range compression of a respective radar pulse, which filtering is matched to the waveform of the respective radar pulse.

3. The method according to claim 1, wherein the interference radar echo of the interference pulse is a nadir echo, the first range position thereof corresponds to the shortest distance between the earth's surface and the azimuth direction.

4. The method according to claim 1, wherein the radar pulses are frequency-modulated and preferably linearly frequency-modulated waveforms.

5. The method according to claim 1, wherein all radar pulses have the same pulse duration and the waveform of an emitted radar pulse corresponds to the waveform of the previously emitted radar pulse shifted cyclically within the pulse duration.

6. The method according to claim 1, wherein the radar pulses are emitted in repeating sequences of successive radar pulses wherein waveforms s.sub.i(t) of the radar pulses within a sequence read as follows: $s_{i} (t) = {\begin{matrix} e^{j \frac{B}{T} {(t - t_{i} - T (.Math. \frac{t + \frac{T}{2} - t_{i}}{T} .Math.))}^{2},} & - \frac{T}{2} t \frac{T}{2} \\ 0, & otherwise \end{matrix}$ wherein i indicates the successive radar pulses of the sequence and i=0, 2BT1 holds, wherein B represents the bandwidth of the radar pulse and T represents the pulse duration; $t_{i} = \frac{i (i + 1)}{2 B} - T .Math. \frac{i (i + 1) + BT}{2 BT} .Math.$ wherein holds.

7. The method according to claim 1, wherein the signal suppression is such that the second range position is determined for a respective radar pulse and the signal is at least partially suppressed in the first focused data essentially only at the second range position.

8. The method according to claim 1, wherein the signal suppression is such that a threshold value is applied to the signals in the entire first focused data, so that all signals above the threshold value are cut off at the threshold value or so that the signals are set to zero only at the positions which exceed the threshold value.

9. An apparatus for the computer-implemented processing of Synthetic Aperture Radar (SAR) raw data, which comprise radar echoes of radar pulses which have been emitted by a radar device on at least one flying object, which is moving in an azimuth direction over the earth's surface, wherein the radar echoes of the radar pulses have been received and recorded by the radar device and wherein the SAR raw data contain data samples for a plurality of range positions for each radar pulse, wherein the range positions represent different distances between the earth's surface and the azimuth direction, wherein an interference radar echo and an interference pulse are associated with a respective radar pulse, wherein the interference radar echo is a radar echo of the interference pulse for a first range position and the interference radar echo is received by the radar device simultaneously with a radar echo of the respective radar pulse for a second range position, wherein the interference pulse is a radar pulse other than the respective radar pulse and wherein the interference pulse and the respective radar pulse have orthogonal waveforms, wherein the apparatus is configured to perform a method, wherein a) the SAR raw data are focused by a first focusing on the interference pulses associated with the respective radar pulses, as a result of which first focused data are obtained, wherein a range compression and an azimuth compression are carried out as part of the first focusing, where a filtering is used for the range compression (RC1) of a respective radar pulse, which filtering is matched to the waveform of the associated interference pulse, b) the first focused data undergo a signal suppression which at least partially suppresses the associated interference radar echo for each radar pulse, as a result of which second focused data are obtained; c) the second focused data undergo defocusing, as a result of which modified SAR raw data are obtained, wherein as part of the defocusing a range decompression and an azimuth decompression are carried out, which are inverse to the range compression and the azimuth compression of the first focusing.

10. The apparatus according to claim 9, wherein the modified raw data are focused on the respective radar pulses by a second focusing, as a result of which third focused data are obtained, wherein a range compression and an azimuth compression are carried out as part of the second focusing, where a filtering is used for the range compression of a respective radar pulse, which filtering is matched to the waveform of the respective radar pulse.

11. The apparatus according to claim 9, wherein the interference radar echo of the interference pulse is a nadir echo, the first range position thereof corresponds to the shortest distance between the earth's surface and the azimuth direction.

12. The apparatus according to claim 9, wherein at least one of: the radar pulses are frequency-modulated or linearly frequency-modulated waveforms; or all radar pulses have the same pulse duration and the waveform of an emitted radar pulse corresponds to the waveform of the previously emitted radar pulse shifted cyclically within the pulse duration.

13. The apparatus according to claim 9, wherein the radar pulses are emitted in repeating sequences of successive radar pulses, wherein waveforms s.sub.i(t) of the radar pulses within a sequence read as follows: $s_{i} (t) = {\begin{matrix} e^{j \frac{B}{T} {(t - t_{i} - T (.Math. \frac{t + \frac{T}{2} - t_{i}}{T} .Math.))}^{2},} & - \frac{T}{2} t \frac{T}{2} \\ 0, & otherwise \end{matrix}$ wherein i indicates the successive radar pulses of the sequence and i=0, 2BT1 holds, wherein B represents the bandwidth of the radar pulse and T represents the pulse duration; $t_{i} = \frac{i (i + 1)}{2 B} - T .Math. \frac{i (i + 1) + BT}{2 BT} .Math.$ wherein holds.

14. The apparatus according to claim 9, wherein the signal suppression is such that at least one of: the second range position is determined for a respective radar pulse and the signal is at least partially suppressed in the first focused data essentially only at the second range position; or a threshold value is applied to the signals in the entire first focused data, so that all signals above the threshold value are cut off at the threshold value or so that the signals are set to zero only at the positions which exceed the threshold value.

15. A non-transitory computer-readable medium storing computer code executable by a processor, comprising one or more codes executable to process Synthetic Aperture Radar (SAR) raw data, which comprise radar echoes from radar pulses, which have been emitted by a radar device on at least one flying object which is moving in an azimuth direction over the earth's surface, wherein the radar echoes of the radar pulses have been received and recorded by the radar device and wherein the SAR raw data contain data samples for a plurality of range positions for each radar pulse, wherein the range positions represent different distances between the earth's surface and the azimuth direction, wherein an interference radar echo and an interference pulse are associated with a respective radar pulse, wherein the interference radar echo is a radar echo of the interference pulse for a first range position and the interference radar echo is received by the radar device simultaneously with a radar echo of the respective radar pulse for a second range position, wherein the interference pulse is a radar pulse other than the respective radar pulse and wherein the interference pulse and the respective radar pulse have orthogonal waveforms, wherein a) the SAR raw data are focused by a first focusing on the interference pulses associated with the respective radar pulses, as a result of which first focused data are obtained, wherein a range compression and an azimuth compression are carried out as part of the first focusing, where a filtering is used for the range compression of a respective radar pulse, which filtering is matched to the waveform of the associated interference pulse, b) the first focused data undergo a signal suppression, which at least partially suppresses the associated interference radar echo for each radar pulse, as a result of which second focused data are obtained; c) the second focused data undergo defocusing, as a result of which modified SAR raw data are obtained, wherein as part of the defocusing a range decompression and an azimuth decompression are carried out, which are inverse to the range compression and the azimuth compression of the first focusing.

16. The non-transitory computer-readable medium according to claim 15, wherein the modified raw data are focused on the respective radar pulses by a second focusing, as a result of which third focused data are obtained, wherein a range compression and an azimuth compression are carried out as part of the second focusing where a filtering is used for the range compression of a respective radar pulse, which filtering is matched to the waveform of the respective radar pulse.

17. The non-transitory computer-readable medium according to claim 15, wherein the interference radar echo of the interference pulse is a nadir echo, the first range position thereof corresponds to the shortest distance between the earth's surface and the azimuth direction.

18. The non-transitory computer-readable medium according to claim 15, wherein at least one of: the radar pulses are frequency-modulated or linearly frequency-modulated waveforms; or all radar pulses have the same pulse duration and the waveform of an emitted radar pulse corresponds to the waveform of the previously emitted radar pulse shifted cyclically within the pulse duration.

19. The non-transitory computer-readable medium according to claim 15, wherein the radar pulses are emitted in repeating sequences of successive radar pulses (RP), wherein waveforms s.sub.i(t) of the radar pulses within a sequence read as follows: $s_{i} (t) = {\begin{matrix} e^{j \frac{B}{T} {(t - t_{i} - T (.Math. \frac{t + \frac{T}{2} - t_{i}}{T} .Math.))}^{2},} & - \frac{T}{2} t \frac{T}{2} \\ 0, & otherwise \end{matrix}$ wherein i indicates the successive radar pulses of the sequence and i=0, 2BT1 holds, wherein B represents the bandwidth of the radar pulse and T represents the pulse duration; $t_{i} = \frac{i (i + 1)}{2 B} - T .Math. \frac{i (i + 1) + BT}{2 BT} .Math.$ wherein holds.

20. The non-transitory computer-readable medium according to claim 15, wherein the signal suppression is such that at least one of: the second range position is determined for a respective radar pulse and the signal is at least partially suppressed in the first focused data essentially only at the second range position; or a threshold value is applied to the signals in the entire first focused data, so that all signals above the threshold value are cut off at the threshold value or so that the signals are set to zero only at the positions which exceed the threshold value.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The novel features believed to be characteristic of the disclosure are set forth in the appended claims. In the descriptions that follow, like parts are marked throughout the specification and drawings with the same numerals, respectively. The drawing figures are not necessarily drawn to scale and certain figures may be shown in exaggerated or generalized form in the interest of clarity and conciseness. The disclosure itself, however, as well as a preferred mode of use, further objects and advances thereof, will be best understood by reference to the following detailed description of illustrative aspects of the disclosure when read in conjunction with the accompanying drawings, wherein:

(2) FIG. 1 shows a diagrammatic representation, based on which the SAR principle used within the scope of the present disclosure is explained;

(3) FIG. 2 shows a flow chart, which illustrates the steps of a variant of the method according to the present disclosure; and

(4) FIG. 3 shows a diagram which shows the waveforms of radar pulses which are used in an embodiment of the present disclosure.

DETAILED DESCRIPTION

(5) FIG. 1 shows in a diagrammatic representation a SAR radar device 1, which is indicated as a rectangular aperture and, depending on the embodiment, comprises one or more antennas. In the represented scenario, the radar device is located on a satellite (not shown), which is moving at a height h along a radar path RT above the earth's surface GR. The direction of the radar path corresponds to the azimuth direction known per se, which is denoted by reference symbol x in FIG. 1. Instead of moving the radar device by means of a satellite, another flying object may be used for this purpose, e.g. an aircraft.

(6) The radar device 1 emits radar pulses RP in successive pulse repetition intervals at a predetermined pulse repetition frequency in an oblique direction to the earth's surface GR. The major part of the energy of a respective radar pulse is directed onto the elliptical area FP on the earth's surface. This area is usually referred to as a footprint of the radar device or of the associated radar antenna.

(7) According to FIG. 1, a respective radar pulse RP has a predetermined pulse duration T, so that the radar pulse has the spatial extension c.sub.0T, wherein c.sub.0 corresponds to the speed of light. As part of an SAR measurement, radar echoes of radar pulses RP scattered back from the Earth's surface are received and recorded by radar device 1 during its movement along path RT. Information about the earth's surface in the swath SW is thus detected. The radar device is configured such that it is arranged for both transmitting radar pulses and receiving corresponding radar echoes. The radio echoes detected by the radar device depend on the form and character of the earth's surface and enable, with a known subsequent signal processing, the calculation of SAR images of the earth's surface.

(8) Before performing the subsequent signal processing, the detected radar echoes are represented, after an analog-digital conversion, by so-called SAR raw data. These raw data are data samples which contain the amplitude and the phase of the sampled radio echoes. The raw data are arranged in a two-dimensional matrix, wherein one dimension of the matrix corresponds to the respective emitted radar pulse (represented by a pulse number) and the other dimension of the matrix represents a time delay, which represents the time period which a sampled radar echo requires to propagate from the radar device to the earth's surface and back to radar device 1. In other words, the so-called slant range R is represented by this time period, said slant range corresponding to the distance between the radar device and the scattering point of the radar echo on the earth's surface. The distance can thus be equated with a direction from which the sampled radar echo passes from the earth's surface GR to radar device 1.

(9) In the scenario of FIG. 1, the reflection of a radar echo at scattering point P on the earth's surface is indicated by way of example. The slant range of this radar echo is denoted by Ro. The slant range is in a geometrical relationship with the so-called ground range (ground distance), which is denoted in FIG. 1 by y and represents the distance between nadir path NT and the corresponding scattering point. The value of a slant range R can thus be unambiguously converted into the value of a corresponding ground range y. Nadir path NT is the vertical projection of radar path RT onto the earth's surface GR. Position N on path NT is the so-called nadir, i.e. the point on the earth's surface with the shortest distance to the radar device. A radar echo from nadir N is referred to as a nadir echo.

(10) As already mentioned, the SAR raw data undergo a subsequent signal processing. Depending on the embodiment, this signal processing can already take place in the satellite, wherein the processed information is then transmitted to a ground station on the earth's surface. Alternatively, the SAR raw data may be transmitted to a ground station without post-processing, the post-processing being carried out in the ground station so as to obtain SAR images from the raw data. The post-processing comprises two filter operations, which are carried out along the range direction R or y and along the azimuth direction x. The filter operation along the range direction is often referred to as a range compression and the filter operation along the azimuth direction is referred to as an azimuth compression. With these operations, a focusing of the SAR raw data is achieved, in order to calculate SAR images.

(11) As part of the range compression, a signal-matched filter is used for the data samples in the range direction of the radar echoes of each radar pulse. In this filter, the data samples along the range direction and therefore along the time axis are convolved with function h(t), which results from the complex-conjugated and time-inverted waveform s(t) of the corresponding radar pulse, i.e. the following holds:
h(t)=s*(t)

(12) By means of this range compression, the SNR ratio (SNR=Signal to Noise Ratio) of the signal is maximized and the scattering point is located at the correct range position.

(13) Besides the range compression, the aforementioned azimuth compression is carried out. Like the range compression, the azimuth compression is known per se and thus will not be explained in detail. The azimuth compression is based on the finding that scattering points on the earth's surface with different azimuth angles relative to radar path RT comprise frequency shifts in the radar echoes on account of the Doppler effect. By analyzing of the Doppler frequency spectrum, a scattering point on the earth's surface can be located in the azimuth direction. The azimuth compression thus leads to a sharpening of the data in the azimuth direction.

(14) The SAR data acquisition has the problem that a significant radar echo is also received from nadir N outside the swath SW. Hence, a radar echo for a specific range position is usually superimposed by the nadir echo of a subsequently transmitted radar pulse, since the nadir echo, due to the shorter travel path, is detected by the radar device at the same point in time as the considered radar echo. A nadir echo occurs for the range position Ro shown in FIG. 1 when the following condition is met:

(15) $\begin{matrix} R_{0} = h + k \frac{c_{0}}{2 PRF} & (1) \end{matrix}$
k denotes an integer value and PRF corresponds to the pulse repetition frequency with which the radar pulses are emitted. h is the height (already mentioned) of the radar device above the earth's surface GR and c.sub.0 corresponds to the speed of light.

(16) Although radar device 1 is intended to concentrate the energy in the swath SW, the nadir echo cannot be prevented due to the special properties of the scattering process (specular reflection). Under certain circumstances, the nadir echo may be stronger than the radar echo on which it is superimposed. The nadir echo can thus significantly degrade the quality of a focused SAR image. The nadir echo appears in the SAR image as a bright stripe at the position on the earth's surface which corresponds to the range position of the radar echo which is disturbed by the nadir echo.

(17) Conventionally, nadir echoes are avoided in SAR images by matching the pulse repetition frequency and the position or width of the detected swath on the earth's surface to one another, so that the nadir echoes do not occur inside the swath. As a result, however, the width of the swath is limited and furthermore the performance of the SAR system with regard to ambiguities is degraded. In order to improve the performance, it may be necessary to increase the complexity of the system, e.g. by increasing the size of the radar antennas.

(18) Within the scope of the present disclosure, a new approach is used in the signal processing of SAR raw data, in order reliably to remove the nadir echo from these data. This new approach is explained below based on FIG. 2. The starting point of the method of FIG. 2 is the originally detected SAR raw data, which are denoted by RD. As already mentioned above, the SAR raw data contain, for each radar pulse, data samples for a plurality of range positions. The raw data for a respective radar pulse are denoted in FIG. 2 by RD.sub.i, wherein i indicates the corresponding radar pulse. As an example, the raw data are represented by an image which contains two objects in the form of a house and a tree, wherein both objects are blurred in the SAR raw data, as indicated by blotches represented dotted. The raw data also contain a nadir echo, which is indicated by a broad vertical stripe at the range position from which the radar echo is received on which the nadir echo is superimposed. This stripe is indicated in FIG. 1 by a square pattern.

(19) To implement the method according to the present disclosure, it has been ensured during the acquisition of the SAR raw data that radar pulses are emitted in repeating sequences of waveforms orthogonal to one another, i.e. each pair of waveforms within a sequence has a scalar product of zero. In the embodiment described here, the waveforms within a sequence are cyclically shifted chirps, wherein a chirp is a linearly frequency-modulated waveform with a predetermined pulse duration. A corresponding chirp s.sub.i(t) is given by the following equation:

(20) $\begin{matrix} s_{i} (t) = {\begin{matrix} e^{j \frac{B}{T} {(t - t_{i} - T (.Math. \frac{t + \frac{T}{2} - t_{i}}{T} .Math.))}^{2},} & - \frac{T}{2} t \frac{T}{2} \\ 0, & otherwise \end{matrix} & (2) \end{matrix}$

(21) T denotes the pulse duration of the chirp, B is the waveform bandwidth (i.e. the bandwidth of the frequencies contained therein) and i represents the i-th pulse within a sequence. The quantity t.sub.i specifies the magnitude of the cyclic shift, wherein t.sub.i is defined as follows:

(22) $\begin{matrix} t_{i} = \frac{i (i + 1)}{2 B} - T .Math. \frac{i (i + 1) + BT}{2 BT} .Math., i = 0, .Math., 2 BT - 1 & (3) \end{matrix}$

(23) A sequence thus contains 2BT pulses, which differ from one another by cyclic shifts. A cyclic shift means that the pulse is shifted in the direction of the time axis and the portion of the pulse which exceeds the limit of the pulse duration at T/2 is added at the front end of the pulse duration at T/2.

(24) The cyclic shift of the pulses just described for achieving orthogonal waveforms is indicated in diagram DI of FIG. 3. Sub-diagram DI1 shows the real part of above pulse s.sub.0(T), whereas diagram DI2 illustrates the imaginary part of this pulse. In contrast, diagram DI3 shows the real part of the pulse for an index i, which corresponds to a shift of 0.3 s, whereas diagram DI4 illustrates the corresponding imaginary part. It can clearly be seen from diagrams DI1 to DI4 that both the real parts and the imaginary parts of the pulse waveforms are shifted within the pulse duration along the time axis. By using the above orthogonal waveforms, it is ensured that the nadir echo in the processing of FIG. 2 has a waveform which is orthogonal to the waveform of the pulse which arrives simultaneously with the nadir echo.

(25) In the signal processing of FIG. 2, a first focusing FC1 of the SAR raw data is first carried out. As part of this first focusing, a first range compression RC1 is carried out for the respective pulses. The range-compressed data are further compressed by a first azimuth compression AC1 in the azimuth direction. Range compression RC1 and azimuth compression AC1 and analogously also compressions RC2 and AC2 as well as corresponding decompressions RDE and ADE are represented in FIG. 2 as separate steps for the sake of greater clarity. These compressions and correspondingly also the decompressions are usually nested into one another, in order to correct further effects such as the range cell migration known per se.

(26) As part of the first range compression RC1, corresponding SAR raw data RD.sub.i of each radar pulse undergo a first filtering FI1. In such a filtering, a temporal convolution is carried out, which is calculated in the frequency domain in the embodiment described here. To do so, the raw data RD.sub.i first undergo an FFT transformation (FFT=Fast Fourier Transformation) and the signal thus contained is then multiplied by the Fourier transformation of a pulse response. In contrast with a conventional range compression, this pulse response is based not on the i-th pulse just processed, but on the radar pulse that belongs to the nadir echo. In the shown scenario, this radar pulse has index (i+k), i.e. it is the k-th radar pulse after the i-th radar pulse.

(27) Accordingly, a multiplication by function H.sub.i+k*(f) takes place during filtering FI1, which function is the Fourier transform of function h.sub.i+k=s.sub.i+k*(t), wherein s.sub.i+k(t) corresponds to the cyclically shifted chirp of the (i+k)-th radar pulse. An IFFT transformation (IFFT=Inverse Fast Fourier Transformation) is then carried out. After filtering all radar pulses, the first azimuth compression AC1 follows, which does not differ from a conventional azimuth compression according to the prior art. Thus, this compression is not described in detail. A particular feature, however, is that this azimuth compression is not matched to the image distance of the useful signal, but rather to the distance of the nadir pulse. First data DA1 are finally obtained in the form of an SAR image, in which the nadir echo is focused without the useful signal being lost. In other words, the data DA1 now contain a nadir echo in the form of a narrow vertical line, whereas the useful data remain almost completely intact in the form of blotches in the image. The useful data are in addition blurred, which is indicated by different hatching of the blotches compared to raw data RD.

(28) The data DA1 then undergo a signal suppression SU. Since the position of the radar echo disturbed by the nadir echo is known from the geometry of the SAR measurement (see equation (1)), the signal is suppressed in a dedicated manner only at the location of the disturbed radar pulse in data DA1 in the embodiment described here. Consequently, the data DA2 are obtained in the form of an SAR image, in which the nadir echo is removed. The samples of the data DA2 for a radar pulse with index i are denoted by DA2.sub.i.

(29) The data DA2 then undergo a defocusing DF known per se, which comprises a range decompression RDE and an azimuth decompression ADE. As part of the range decompression, an inverse filtering FI2 is carried out. This takes place again in the frequency domain. Hence, an FFT transformation of the samples DA2.sub.i is first carried out, followed by a multiplication by the reciprocal of above function H.sub.i+k(f) as well as by a subsequent IFFT transformation. Besides the range decompression RDE, the azimuth decompression ADE is applied to the data in a manner known per se. Since this azimuth decompression takes place with methods known per se, it is not described in detail here.

(30) After the defocusing DF, modified raw data RD are obtained, which apart from a minimal degradation correspond to original raw data RD and from which the nadir echo is removed. These modified raw data RD then undergo a second focusing FC2, in which however a filtering F13 is now used as part of range compression RC2, which filtering is matched to the respective radar pulses. As part of this filtering, the modified raw data RD.sub.i of a respective radar pulse again undergo an FFT transformation, wherein a multiplication by function H.sub.i*(f) then takes place, which is the Fourier transform of function h.sub.i(t)=s.sub.i*(t), wherein s.sub.i(t) corresponds to the chirp of the i-th radar pulse. Finally, the signal is again transformed with an IFFT transformation into the time space. Apart from second range compression RC2, a second azimuth compression AC2 is carried out in a manner known per se.

(31) As a result of the second focusing FC2, third data DA3 are finally obtained, which correspond to a focused SAR image, which no longer contains the nadir echo. As can be seen from FIG. 2, both objects in the form of a house and a tree can now be seen in the image, without the image containing a disturbing stripe that is caused by the nadir echo.

(32) The embodiment of the present disclosure as described above has a number of advantages. In particular, as part of the processing of SAR raw data, a first focusing on the nadir echo is carried out, whereby the filtering is matched to the waveform of the pulse whose radar echo corresponds to the nadir echo. In this way, the nadir echo is focused and can be removed in the data thus obtained by a suitable signal suppression. By means of a subsequent defocusing, modified raw data are obtained without a nadir echo, which raw data can then be converted with a second focusing, now matched to the useful signal, into the SAR image, which then no longer contains a nadir echo. It is therefore no longer necessary to observe any limitations with regard to the pulse repetition frequency of the radar pulses and the nadir echo can essentially be completely removed. By using orthogonal waveforms, it is at the same time ensured that the information from the useful data remains intact with a negligibly small degradation.

(33) The above embodiment of the present disclosure has been explained based on a suppression of a nadir echo. Nonetheless, the present disclosure can also be used to remove so-called range ambiguities from the SAR raw data. Range ambiguities occur when radar echoes are received from different range positions within the swath detected from the earth's surface. For a range position, a range ambiguity can occur at one or also at a plurality of different range positions. In order to suppress such range ambiguities for a predetermined range position (with associated i-th radar pulse), the steps of FIG. 2 up to the acquisition of the modified radar echo RD are carried out iteratively for each range ambiguity, which is then represented by an interference pulse with the index (i+k) (wherein k may possibly also be negative). The raw data RD are equated with the modified raw data RD when processing of a new range ambiguity. The above steps can also be repeated iteratively for a plurality of range positions, at which respective ambiguities exist at other range positions. In contrast with the signal suppression for the nadir echo described above, the suppression of range ambiguities is based on a threshold value, wherein all signals which lie above the threshold value within data DA1 are cut off at the threshold value. Alternatively, the signals are set to zero only at the positions which exceed the threshold value.

(34) Although the operations or methods described above are presented in a particular order and/or as being performed by an example component, it should be understood that the ordering of the actions and the components performing the actions may be varied, depending on the implementation. In addition, aspects of any one of the methods described above can be combined with aspects of any other of the methods.

(35) The above detailed description set forth above in connection with the appended drawings describes examples and does not represent the only examples that may be implemented or that are within the scope of the claims. The term example, when used in this description, means serving as an example, instance, or illustration, and not preferred or advantageous over other examples. The detailed description includes specific details for the purpose of providing an understanding of the described techniques. These techniques, however, may be practiced without these specific details. In some instances, well-known structures and apparatuses are shown in block diagram form in order to avoid obscuring the concepts of the described examples.

(36) Information and signals may be represented using any of a variety of different technologies and techniques. For example, data, instructions, commands, information, signals, bits, symbols, and chips that may be referenced throughout the above description may be represented by voltages, currents, electromagnetic waves, magnetic fields or particles, optical fields or particles, computer-executable code or instructions stored on a computer-readable medium, or any combination thereof.

(37) The various illustrative blocks and components described in connection with the disclosure herein may be implemented or performed with a specially-programmed device, such as but not limited to a processor, a digital signal processor (DSP), an ASIC, a FPGA or other programmable logic device, a discrete gate or transistor logic, a discrete hardware component, or any combination thereof designed to perform the functions described herein. A specially-programmed processor may be a microprocessor, but in the alternative, the processor may be any conventional processor, controller, microcontroller, or state machine. A specially-programmed processor may also be implemented as a combination of computing devices, e.g., a combination of a DSP and a microprocessor, multiple microprocessors, one or more microprocessors in conjunction with a DSP core, or any other such configuration.

(38) The functions described herein may be implemented in hardware, software executed by a processor, firmware, or any combination thereof. If implemented in software executed by a processor, the functions may be stored on or transmitted over as one or more instructions or code on a non-transitory computer-readable medium. Other examples and implementations are within the scope and spirit of the disclosure and appended claims. For example, due to the nature of software, functions described above can be implemented using software executed by a specially programmed processor, hardware, firmware, hardwiring, or combinations of any of these. Features implementing functions may also be physically located at various positions, including being distributed such that portions of functions are implemented at different physical locations. Also, as used herein, including in the claims, or as used in a list of items prefaced by at least one of indicates a disjunctive list such that, for example, a list of at least one of A, B, or C means A or B or C or AB or AC or BC or ABC (i.e., A and B and C).

(39) Computer-readable media includes both computer storage media and communication media including any medium that facilitates transfer of a computer program from one place to another. A storage medium may be any available medium that can be accessed by a general purpose or special purpose computer. By way of example, and not limitation, computer-readable media can comprise RAM, ROM, EEPROM, CD-ROM or other optical disk storage, magnetic disk storage or other magnetic storage devices, or any other medium that can be used to carry or store desired program code means in the form of instructions or data structures and that can be accessed by a general-purpose or special-purpose computer, or a general-purpose or special-purpose processor. Also, any connection is properly termed a computer-readable medium. For example, if the software is transmitted from a website, server, or other remote source using a coaxial cable, fiber optic cable, twisted pair, digital subscriber line (DSL), or wireless technologies such as infrared, radio, and microwave, then the coaxial cable, fiber optic cable, twisted pair, DSL, or wireless technologies such as infrared, radio, and microwave are included in the definition of medium. Disk and disc, as used herein, include compact disc (CD), laser disc, optical disc, digital versatile disc (DVD), floppy disk and Blu-ray disc where disks usually reproduce data magnetically, while discs reproduce data optically with lasers. Combinations of the above are also included within the scope of computer-readable media.

(40) The previous description of the disclosure is provided to enable a person skilled in the art to make or use the disclosure. Various modifications to the disclosure will be readily apparent to those skilled in the art, and the common principles defined herein may be applied to other variations without departing from the spirit or scope of the disclosure. Furthermore, although elements of the described aspects and/or embodiments may be described or claimed in the singular, the plural is contemplated unless limitation to the singular is explicitly stated. Additionally, all or a portion of any aspect and/or embodiment may be utilized with all or a portion of any other aspect and/or embodiment, unless stated otherwise. Thus, the disclosure is not to be limited to the examples and designs described herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Method and apparatus for computer-assisted processing of SAR raw data

Assignee

Inventors

Cpc classification

Classification Explorer

G01S7/2813

PHYSICS

Classification Explorer

G01S13/90

PHYSICS

Classification Explorer

G01S13/9011

PHYSICS

Classification Explorer

G01S7/295

PHYSICS

International classification

Classification Explorer

G01S13/90

PHYSICS

Classification Explorer

G01S7/295

PHYSICS

Classification Explorer

G01S7/28

PHYSICS

Abstract

Claims

Description