1. Introduction
Synthetic aperture radar (SAR) comprises active microwave antennas mounted on airborne or spaceborne platforms that can achieve high-resolution images of the illuminated area regardless of an external source and weather conditions. They have been used for more than forty years in different space missions to monitor natural and anthropogenic phenomena: earthquakes, landslides, forest growth, desertification, glacier and sea-ice state, urban development, agriculture extension, deforestation, and many other applications [
1].
With the increasing number of applications, two major requirements arose from SAR users: higher-resolution images and, at the same time, a wider swath [
2,
3,
4]. Indeed, we passed from the 30 × 25 m
2 resolution in the 100 km of the swath of the ERS-1 satellite (launched in 1991) to the 5 × 20 m
2 resolution on the 250 km of swath in the still-operative Sentinel-1, with an improvement of 6- and 1.25-times in range and azimuth resolutions and 2.5-times in the swath width.
The contrast between the two requirements lies in the need to reduce the pulse repetition frequency (PRF) to allow the signal from a wider swath to be received, with the need to increase the PRF to allow a wider azimuth bandwidth (and so, a higher resolution) to be achieved without ambiguities [
5,
6].
In new-generation systems such as ROSE-L and Sentinel-1 Next Generation [
7,
8], a high azimuth resolution and low ambiguities are achieved by the use of multiple channels, i.e., multiple sampling during acquisition, and on-ground reconstruction of the full azimuth spectrum in the
N-times PRF non-ambiguous interval. This is the well-known multiple aperture processing scheme (MAPS) technique, widely discussed in [
2,
9,
10,
11,
12,
13].
A high range of resolutionsis commonly achieved by increasing the range bandwidth, a resource usually not expensive in SAR systems, apart from the increase of captured noise. However, imaging a wide swath poses another problem: often, the antenna beam in elevation is too narrow to illuminate the complete swath. This brings a reduced directivity typically at swath borders and, so, the degradation of sensitivity.
Among the solutions recently proposed to overcome this limitation, we recall SweepSAR [
14] and SCORE [
15]. Both techniques illuminate the entire swath upon transmission, whilst, upon receiving, a high-directivity pencil beam scans the swath, tracking the locus of the return echo. SweepSAR is based on an analogue beam steering network; SCORE is instead a digital beam forming (DBF) method based on multiple receivers in elevation, each with its own digitized receiving chain. In SCORE, the wide transmit beam is generated by the full antenna height by a proper excitation of the amplitudes and phases of sub-array elements of the phased array antenna [
16]. In reception,
broader beams are used by grouping subsets of sub-arrays in elevation. These beams are onboard composed to obtain an equivalent narrow beam pointing in the current position of the receiving echo. The notable element of this technique is that the composition is performed in the digital domain, so that the narrow beam can be steered at the time of the receiving echo, i.e., in the order of microseconds. SCORE then allows the directivity to remain close to its maximum value throughout the swath length. This way, the sensitivity is optimized to its highest achievable value.
A recent upcoming alternative to SCORE is represented by the frequency scanning principle. The frequency scanning principle was firstly presented in [
17,
18] (Airbus owns also a patent on this method). f-SCAN is based on the idea of obtaining the same advantages of the SCORE technique, but using fewer demanding electronics. The SCORE technique requires: (i) the independent acquisition and digitization of
channels and (ii) an onboard recombining network whose coefficients must change according to the fast time sampling frequency. In f-SCAN, the same result is achieved by using phase shifters (PSs) and true time delay lines (TTDLs) on the phased array antenna. Since the beam composition is completely analogous, there is no need to maintain multiple channels and, then, multiple samplers, ADCs, and so on.
The opposite of SCORE, the pencil beam is generated in both TX and RX and always using the whole antenna. The pulse transmission is performed as in a conventional SAR system, but the frequency dispersion is fruitfully exploited to gather different parts of the transmitted chirp (which indeed correspond to different frequencies) at different ground ranges. The design of the frequency rate of the pulse (i.e., an “up” or “down” chirp) is performed in such a way that, for the chirp’s frequency transmitted at first, the pencil beam must point at the far range; conversely, for the last transmitted frequencies, the pencil beam must point at near range. The same principle is valid also for the received echoes. A consequence is that the received signal is shorter than what is expected from geometry, due to the directional filter operated by the TX and RX beams.
Thanks to the narrow beam used for scanning, both the transmit and receive beams have higher directivity, which means a better noise equivalent sigma zero (NESZ) and better range-distributed target-to-ambiguity ratio (RgDTAR). As will be shown, the main drawback of this method is that f-SCAN requires a larger bandwidth in transmission and a higher duty cycle.
In the recent literature, the principle of frequency scanning and frequency dispersion have been addressed. In [
19,
20], the concept of the conventional SCORE was extended to wide-pulse systems by proposing digital scalloped beam forming. This method has been revealed to be an effective method to soften the influence of frequency dispersion (which indeed is still interpreted as an unwanted effect) in DBF SAR. A similar result was achieved also in [
21]. In [
22], a novel SAR operational mode named the multi-frequency sub-pulse (MFSP) was presented: the wide-pulse bandwidth was exploited to increase the imaged swath extension without the emergence of range ambiguities. References [
23,
24] gave a detailed description of the implementation of the frequency scanning principle with a particular focus on the antenna, but the analysis of the corresponding SAR image capability was missing. Again, in the papers [
17,
18] by Roemer, where the f-SCAN principle was applied to SAR for the first time, there was not yet an exhaustive characterization of the constraints to be considered while designing a system. Moreover, the evaluation of SAR performance was only quickly introduced. In the more recent paper by Scheiber et al. [
25], the focus was on the transmit pulse chirp length, which was interpreted as a trade-off parameter: the consequences on the echo window length and the benefits of dedicated onboard processing for data volume reduction were discussed. Eventually, in [
26], a further version of the f-SCAN SAR, named elevation frequency scanning synthetic aperture radar (EF-SCAN SAR), was presented. In this work, the analytic space–time signal model of an intrapulse frequency scanning array was derived. Based on this, an imaging algorithm for EF-SCAN SAR was developed.
In this paper, a detailed general methodology to design an f-SCAN acquisition system is presented, giving more details than in the existing literature in the previous paragraph. The design was applied in the case of an exemplary demanding X-band spaceborne SAR acquisition mode. The processing of f-SCAN data was addressed by presenting a novel (to the authors’ knowledge) algorithm that aimed to guarantee an invariant IRF throughout the swath by removing spectral distortions. The definitions of NESZ and RgDTAR were accordingly updated. Eventually, a numerical simulation was implemented to estimate the performance figures: IRF resolution, peak and integrated sidelobe ratios (PSLR, ISLR), NESZ, and RgDTAR.
The paper is organized as follows. In
Section 2, the frequency dispersion concept and the f-SCAN system and its design are explained. In
Section 3, the f-SCAN processing is described, while in
Section 4, the simulation results of a case study are shown and discussed. In
Section 5 the theoretical performances of the f-SCAN method are defined and compared to a reference system. Eventually, conclusions in
Section 6 close the paper.
2. f-SCAN System and Design
In the present section, the description of the f-SCAN principle and the main design parameters and constraints are defined. The section is organized as follows:
Section 2.1 introduces the frequency dispersion principle;
Section 2.2 describes how to realize f-SCAN, with particular attention to the TTDL; in
Section 2.3, the definition of the elevation pattern is provided; in
Section 2.4, the timing of the instrument is characterized; in
Section 2.5, the conditions for data sub-sampling are derived.
2.1. The Frequency Dispersion Principle
A common method to achieve the antenna pattern shape and proper steering is by using complex coefficients for the excitation of the radiating elements of an array, in both TX and RX. The embedded patterns generated by the radiating elements are usually a function of the elements’ characteristics (geometry, size, frequency), but are quite large due to their limited size. The scan direction
of a linear array is related to a phase increment
introduced in the excitation coefficients used for the different elements by the following relation [
16]:
where
f is the frequency of the signal,
the distance between the elements of the array, and
c the e.m. speed.
For wide-band emitted signals, however, the same phase shift
generates different scan directions of the array factor since the phase shift and direction coincide only for a single frequency. In practice, for a monochromatic signal,
, applying a delay or a phase produces the same effect:
and
, where
. For a non-monochromatic signal, the application of a constant phase shift will produce different delays at different frequencies. This effect is known as frequency dispersion [
16]. In
Figure 1, an example of frequency shift due to the use of phase shifters only at different frequencies is shown. The modulation effect of the directivity is due to the embedded pattern of the single element.
In radar technology, a slight scanning effect of frequency is a well-known (and unwanted) problem for wide-band signals where phase shifters are used as an approximation of the (true) delay lines [
17].
The frequency dispersion effect can be used as an advantage in the f-SCAN phase arrays. The basic idea is to introduce in the phase increment of the array contribution (i.e., in Equation (
1) a delay line of amount
[
17]:
If the previous relation is written as a function of the scan direction
[
17]:
The scan direction results in a function of two contributions: (i) a frequency-independent term, which is commanded by a true delay increment, and (ii) a frequency-dependent term, commanded by the phase shifter.
In f-SCAN, the two degrees of freedom, i.e., the delay and the phase shifter, are used together to obtain the desired frequency dispersion. In detail:
The phase shifter is used to steer the array factor toward the desired angle, i.e., the off-nadir centre of the swath for a SAR system. The resulting frequency dispersion of the peak is computed accordingly.
The delay lines produce a non-dispersive shift of the peak, but also a stretching of the pattern (because of the bandwidth of the pulse), which actually results in a frequency dispersion on the grating lobe. This unwanted effect is used to complete the frequency dispersion caused by the phase shifters so as to cover the whole angular interval of the swath.
In the following section, the combined effect of the two is detailed.
2.2. Design of Phase Shifters and Delay Lines in f-SCAN
Let us design a SAR acquisition mode in which a swath with central off-nadir angle and width must be fully scanned in both TX and RX by a pencil beam. Let us also assume a planar array antenna with the following characteristics:
is the elevation size (i.e., height) of the whole array antenna.
N is the number of elements along the elevation.
is the number of TTDLs. Then, the N elements from the whole antenna are grouped into sub-arrays of elements each.
According to the first point of the strategy, defined at the end of the previous section, the phase shift
was designed to provide an illumination towards the swath centre
at the carrier frequency
. The relation between
and
, already given in Equation (
1), is here repeated [
16]:
Due to the finite bandwidth of the signal,
, the corresponding frequency dispersion is given by inverting Equation (
4) with the chirp frequency edges:
Usually, the angle interval does not coincide with the swath angular extension , which is dictated by other design considerations (if it coincides or is greater, this means that there is no need for further dispersion, and then, the TTDLs are useless).
For this reason, according to Step 2 of the strategy defined in
Section 2.1, a set of true time delay lines was designed to achieve the full coverage of the swath. The TTDLs do not create mispointing on the main lobe, but create a frequency dispersion on the grating lobes because of the stretching of the angle axis (see
Figure 2, left panel). The position of the first grating lobe depends on the wavelength
. It is derived from the expression of the periodic sinc representing the array factor of the array antenna [
16]:
where we recall that
is the antenna height and
is the number of available delay lines.
The frequency dispersion at the
k-th grating lobe is then:
Please note that the previous approximation is valid until the quantities in the inverse sine are little. The frequency dispersion is then larger for farther grating lobes (i.e., the higher
k, the larger the dispersion), as shown in
Figure 2, right panel.
The angular interval of the swath not covered by the frequency dispersion due to the phase shifters (see Equation (
5)) is covered by the stretching effect due to the use of delay lines (see Equation (
7)). A good approximation of the integer number of grating lobes to be shifted is (for the exact relation, please refer to
Appendix A):
Eventually, the delay for each line is derived by combining Equations (
3), (
6), and (
8):
In Equations (
9),
is the distance between the phase centre of the group of radiating elements, which share the same delay line.
2.3. Generation of the Elevation Pattern
The elevation pattern of the phased array antenna and due to phase shifters and delay lines is [
16]:
In the previous equation, , is the common embedded pattern from each radiating element, while is the variable wavelength.
The number of delay lines is a critical design parameter and may also represent a constraint coming from the hardware. If the number of available delay lines is reduced (i.e., decreases), then:
The angular position of the grating lobes (
) comes closer to the main lobe, as comes from Equation (
6), and also, their dispersion becomes smaller. This means that the integer number of the grating lobes must increase to achieve the same angular dispersion required in Equation (
8). Increasing
k means a higher value of the delay, which impacts its hardware implementation.
With fewer TTDLs, the sub-array composed pattern is narrower and the corresponding attenuation around the grating lobe is stronger. This may jeopardize the performances (e.g., NESZ) at the edge of the swath.
In
Figure 3, there is a detail of the patterns in two cases, using 8 or using 4 delay lines. In the first case, the second grating lobe was used, while the fourth grating lobe was used in the second case. Even the amount of delay doubled, passing from 0.205 ns to 0.409 ns for the specific example, which shall be detailed in
Section 4.
2.4. Relation between Chirp and Timing
While the chirp is being transmitted, the pencil beam scans the swath so as to illuminate mainly a given portion of the underneath area. This means that each target in the swath shall be (mainly) illuminated only by a part of the chirp. This useful part of the chirp is driven by the desired range resolution of the system and is shaped, in the frequency domain, by the scanning pencil beam of f-SCAN. The principle can be effectively described using a time–frequency diagram (T-F), as the one commonly used for the slow-time Doppler relation in azimuth, with the difference that, now, it is applied to the fast time and range frequency.
To describe the f-SCAN acquisition in such diagrams, we refer to
Figure 4. We assumed a chirp pulse of duration
and total bandwidth
. Indeed, these parameters are often fixed by the SAR system electronics and mode design (e.g., PRF). The chirp rate is the ratio between the total bandwidth and the chirp duration:
. In the T-F diagram, the chirp signal is represented by the green lines with a negative slope (a down chirp was assumed).
The scanning beam is instead illustrated by the reddish strip with positive slope and duration . Its width, , is the beam width in the frequency domain.
On the other hand, the geometric range resolution is a user requirement that determines the band
B. The relation between
and
B is illustrated in
Figure 4;
B is the vertical projection of the intersection between the chirp line and the steering strip.
If the chirp rate and the steering rate are chosen with opposite signs, the length of the receiving window can be made significantly shorter than the one defined by the geometry. In practice, it is possible to discard the parts of the signal already filtered out by the pencil beam, bringing an effective reduction of the receiving length and of the data volume.
For a given reference geometry, say
, the slant ranges at the near and far range, the corresponding delays are:
. The geometric sampling window start time (
) depends on
(it is the fractional part of
once the integer numbers of
are removed), and the geometric sampling window end time (
) depends on
, leading to the definition of the geometric sampling window length (
):
To allow the correct focusing of the very-far-range target, the geometric SWL (
, Equation (
11)) of a conventional instrument must be extended by the duration of one chirp:
Thanks to the filtering effect of the f-SCAN beam (see
Figure 4), the sampling window length of the f-SCAN system (
) is shorter than what is stated in Equation (
12):
After Equation (
13), the typical one-to-one correspondence between the raw data fast time and slant range distance is lost, i.e., the fast time of the raw data is
compressed with respect to the physical delay of the targets. The advantage is that it is easier to fit the sampling window in the timing diagram. On the other hand, it is necessary to
decompress the data when performing the range focusing. This aspect is addressed in
Section 3.
Please note that, according to
Figure 4 and Equation (
13), a portion
of the full chirp length, enough to retrieve the needed resolution, is always ensured, even for the targets at the very near and far ranges of the swath. The time
is the integration time, which corresponds to the resolution band
B:
Eventually, the steering rate
can be evaluated. Here, it is defined in Hz/s (as the chirp rate) since we were interested in preserving the range bandwidth
B, even in the reduced sampling window length of f-SCAN. From the geometry in
Figure 4, the steering rate
is:
where the scanning time
is:
A remarkable aspect when designing an f-SCAN system is the need to have an optimal weighting of the range-variant central frequency (namely: the red dashed line in
Figure 4) by the scanning pencil beam, so as to gather the maximum energy. In this perspective, the scanning beam is a further degree of freedom, whose design was already addressed in
Section 2.2. Further insight on this point is detailed in
Appendix A.
Equation (
15) states that the SAR signal has a range-dependent central frequency. This peculiarity of f-SCAN mode deserves attention in the design of an interferometric SAR since the spectral overlap is necessary to perform interferometry over distributed targets [
27].
2.5. Data Sub-Sampling and Shrinking Factor
Looking at
Figure 4 and Equation (
13), the 2D support of the f-SCAN system is given by the product:
, where
is the range sampling frequency. Nevertheless, it is evident that most of the domain of the data is empty, i.e., only the diagonal reddish shape contains the useful signal. Therefore, there is room for optimizing the data sampling of the received signal. We can fruitfully exploit the filtering effect of the pencil beam by allowing aliasing from those frequencies attenuated by the pencil beam itself.
Furthermore, the received SAR raw data must be sampled at a frequency: (i) significantly lower than the one required by the full band
of the transmitted chirp, but (ii) still high enough to preserve the information in the resolution band
B. Indeed, the lower bound of the sampling frequency is given by the parameter
in
Figure 4, rather than by
B.
depends not only on the resolution bandwidth
B, but also on both the chirp rate
and steering rate
. From the geometric considerations on
Figure 4 and from Equation (
14), we can write:
Then, by setting equal the right sides of Equation (
17), we obtain the relation between
and
B:
We define as the shrinking factor of the f-SCAN acquisition the quantity:
Equation (
19) states that the target, because of the steering rate
of the pencil beam, experiences an antenna sharper than the physical width of the pencil beam itself. This effect is very similar to the one observed by TOPS [
28] in azimuth. Moreover, Equation (
19) is important because it clarifies that the minimum sampling frequency of the SAR data must be defined considering the inverse of this shrinking factor. In fact,
is always greater than the band
B necessary to fulfil the requirement of the geometric resolution. Eventually, we can take some margin
on
to have the minimum sampling frequency
:
The T-F diagram of the SAR data using the reduced sampling frequency in Equation (
20) is in
Figure 5. The reduction of the 2D support comes from limiting both the
SWL and the sampling frequency
. Compared to this, in
Figure 5, there are the light blue rectangle (i.e., the extended original domain) and the blue rectangle in the centre (i.e., the new reduced domain).
The drawback of using a reduced sampling frequency is the creation of some aliasing, which would lead to “ghost” targets (which is a novel contribution to the overall ambiguity level, not to be confused with the canonical range ambiguities measured by RgDTAR). Nevertheless, the level of these ambiguous targets can be kept under control by properly choosing the sampling frequency and, in the case of necessity, giving a proper shape to the pencil beam, e.g., using tapering on the amplitude RX excitation coefficients. It is worth noticing that such a phenomenon is not new for SAR data as it happens along the azimuth every time the sampling frequency (i.e., PRF) is lower than the Nyquist condition. Therefore, in the instrument design stage, it would be a good choice to have a level of the range ghost comparable to (or slightly lower than) the azimuth ambiguities.
3. f-SCAN Processing
In the current section, the description of the processing of the f-SCAN data is provided. The algorithm hereinafter presented shall be applied to each received sampling window (i.e., over the fast time) at each PRI. As illustrated in
Figure 5, the data are characterized by a fast time-dependent spectrum: targets at different ranges experience different central frequencies. This is the consequence of combining a wide-pulse bandwidth with a narrow pencil beam, which scans the swath using the highest directivity. Furthermore, according to the choice described in
Section 2.5, the raw data spectrum is folded in frequency.
The processing of the f-SCAN data consists of: (i) restoring the full time–frequency support of the data; (ii) performing the typical SAR range compression.
The processing steps needed to restore the original frequency support of the data are shown in the block diagram of
Figure 6 and hereinafter described. We remark that this algorithm shall be applied over fast time at each PRI.
First, the Fourier transform of the signal is taken and
copies of the spectrum are stacked. Note that the mosaicking operation is equivalent, in the time domain, to the zero interleaving of
samples. According to
Figure 5, the number of spectral supports mosaicked is:
The second step is a time-variant demodulation of the mosaicked data. This operation is known as deramping and is performed by multiplying in the time domain the SAR data times a chirp. Referring to
Figure 4, the modulation rate of the deramping function is the steering rate
:
with
. The third step is the application of a low-pass filter (LPF) to select only one spectral replica (the one in the baseband) among the
. Please note that: (i) the nominal bandwidth of the LPF is
, according to Equation (
18), and (ii) together with the LPF, it is possible to apply a whitening filter (to compensate any spectral distortion) and/or a spectral windowing to give the spectrum the desired shape. This is commonly used in SAR processing to regularize the shape of the IRF The last step to accomplish the frequency unfolding is to reramp the data by multiplying it, in the time domain, by the conjugate of Equation (
22).
Once the frequency unfolding of the data is completed, the data time support is extended from
given by Equation (
13) to
given by Equation (
12). This is obtained by padding the data with
samples at the beginning and at the end, where
, according to Equation (
13), is defined as:
Please note that, after frequency spectrum restoration, is the new sampling frequency.
At this stage, the extended support of
Figure 4 has been restored from the compact support of
Figure 5, and the traditional range focusing procedure can be applied, i.e., the correlation of the signal with a copy of the transmitted chirp.
Hereinafter, a mathematical description of the effects of range focusing on the data is explained. The aim was to find an equation that puts into relation the acquisition time of the raw data (
) with the time of the range-focused data (
). We can refer to
Figure 7, where the reddish shape highlights the support of the raw data (after the unfolding procedure), whilst the yellowish shape is for the support of the range-focused data. Intuitively, the relation between the two shapes is the following: the linear dependency between time and frequency, typical of the chirp signal, is removed by the range focusing processing, and all the frequencies in the band (
B) collapse in one single time instant (i.e., the time of range-focused data,
). Hence, we can derive the mathematical relation between
and
by defining the quantity
(highlighted in magenta in
Figure 7) as:
Then, by making the right sides of Equation (
24) equal, we obtain the relation between
and
:
where the first term takes into account the stretching of the time axis, while the second term is a constant offset.
Different from conventional SAR, for f-SCAN acquisition mode, the spectra of targets at different ranges are disjoint and not fully overlapped. This is evident from
Figure 7 where the yellowish shape has a central frequency that varies with the time
. The central frequency
of the target at time
is computed taking advantage of Equation (
25) and
as:
Thus, from the first term of Equation (
26), we obtain the expression of the rate
, which is the frequency modulation rate of the range-compressed data and, indeed, corresponds to the slope of the yellowish shape in
Figure 7:
Equation (
27) can be conveniently used to simplify the azimuth focusing of the f-SCAN data with respect to the processing proposed in [
26]: the advantage is using a conventional (and already-existing) azimuth focusing processor. After having performed range focusing according to the algorithm previously described, the following steps are needed:
Deramping each azimuth line from the range-compressed data according to
, where
is from Equation (
27);
Applying a standard azimuth focusing processor [
29];
Reramping each azimuth line with its natural modulation rate from Equation (
27).
It is worth noting that, for the generation of standard range-compressed data and single-look complex data, it is not necessary to perform the time-variant demodulation along the range (deramping) using Equation (
27) as the extended spectral support allows preserving the physical frequency at which the target has been acquired. Moreover, such residual modulation does not affect the amplitude of the single-look complex, but is still present in its phase. Therefore, all the characteristics of the target response are still preserved in the single-look complex data. On the contrary, for higher-level products (e.g., detected data), it could be useful to perform the deramping along the range, using the rate in Equation (
27) to re-align the spectra of the targets along the baseband. Then, a final sub-sampling can be performed, reducing the volume of the data.
As a final remark, it is worth pointing out that, here, the most intuitive, but not optimized processing of the f-SCAN data is described. Indeed, taking advantage of the similarity between f-SCAN and the TOPS system, it is possible to define an optimized processor such as the one in [
28], which directly transforms the data from the compact support of
Figure 5 to those in
Figure 7.
4. Exemplary System Design
In the following, the design of a realistic system is shown and the time-frequency diagrams in different stages are explained.
The system is composed of an X-band antenna acquiring a stripmap swath of about 40 km on the ground. The details of the system and geometrical parameters are given in
Table 1 and
Table 2, respectively. From the user’s point of view, the incidence angle interval and the satellite height are sufficient to define a (simplified) geometry. The remaining geometric parameters in the table (slant range, swath extension, and so on) were retrieved using geometrical relations and the local spherical Earth assumption.
By following the relations in
Section 2, it is possible to obtain the timing design, as in
Figure 4. Timing values are reported in
Table 3. Please note that the chirp slope
is negative. The system timing, i.e., the sequence of pulse transmission and receiving intervals, together with the true time delay and phase shifter values, is finally given in
Table 4.
Successively, a MATLAB simulator was developed and set up with 11 ideal point targets aligned at the same azimuth and equally spaced in the off-nadir interval, from 19.9 deg to 23.68 deg (border effects were excluded). Each target was illuminated by a portion of the complete chirp due to being set in a given geometric position. Due to the antenna pattern shape (the embedded pattern of each radiator element was supposed to be a sinc shape), the frequency interval of each target was also pattern-shaped, as illustrated in
Figure 8, left panel. The colour code in the figure is proportional to the frequency interval: hot colours for high frequencies; cold colours for low frequencies. Since each target generates a signal response at a different time and with different spectral support, the time–frequency diagram resulted in a sequence of lozenge supports, as shown in
Figure 8, right panel.
A subsampling factor equal to 3 was applied during the acquisition, so that the sampling frequency, initially set to 1.5-times the chirp bandwidth, was set to 600 MHz. Due to spectrum folding, aliasing occurs. However, the spectrograms of the different targets generate different aliasing patterns. This means that, for some of the subsampled spectrograms, the folded tails of the lozenge may overlap on the main part; in other spectrograms, they may not. This superposition, a function of the chirp and f-SCAN slopes, and the chosen sub-sampling factor make the aliasing unavoidable. This is the cost to pay to keep the range sampling frequency and the data volume low during the acquisition.
Raw data have spectral support, as reported in
Figure 9. Before the range compression, they were pre-processed according to the algorithm proposed in
Section 3 and illustrated in
Figure 6. The steps are: (i) spectral mosaicking (by time domain zero interleaving) to raise the frequency; (ii) deramping; (iii) low-pass filtering and whitening; (iv) reramping. The (optional) whitening step was introduced to remove the antenna pattern shape of the target spectral support shown in
Figure 8. Of course, whitening will raise the sidelobes, making the main lobe narrower as in traditional SAR processing, and this may be followed by a proper windowing. Conversely, whitening and windowing give the advantage of achieving an invariant IRF wherever the position of the target in the swath.
The range line after pre-processing can be compressed as in a traditional SAR system, i.e., correlating the range line with the reference pulse replica. The compressed range line is shown in
Figure 10. As expected, the effect of subsampling generated spurious peaks (i.e., “ghosts”) whose level was, anyway, acceptable, being 30 dB lower than the peak. Due to the whitening (no windowing was applied in the simulation), the peak had the usual sinc shape (see
Figure 10, right panel).
Quality measures of the IRF, i.e., ground range resolution, peak-to-sidelobe ratio (PSLR), and integrated sidelobe ratio (ISLR), are finally shown in
Figure 11. As can be seen, the ground range resolution is within the user requirements (i.e., 1.2 m at the near range), and the PSLR and ISLR are very close to the values achieved from an ideal IRF.