Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Study of the Buried Basin C-H, Based on the Multi-Source Remote Sensing Data
Previous Article in Journal
Evaluating Trade-Off and Synergies of Ecosystem Services Values of a Representative Resources-Based Urban Ecosystem: A Coupled Modeling Framework Applied to Panzhihua City, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design of f-SCAN Acquisition Mode for Synthetic Aperture Radar

by
Pietro Guccione
1,2,*,
Daniele Mapelli
1,
Davide Giudici
1 and
Adriano Rosario Persico
1
1
Aresys SRL, 20132 Milan, Italy
2
Department of Electrical Information Engineering, Politecnico di Bari, 70126 Bari, Italy
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(20), 5283; https://doi.org/10.3390/rs14205283
Submission received: 31 July 2022 / Revised: 18 October 2022 / Accepted: 19 October 2022 / Published: 21 October 2022

Abstract

:
This paper presents the design and processing of the SAR acquisition technique named frequency scanning (f-SCAN), aimed to obtain high sensitivity to targets with low backscattering and to improve the signal-to-noise ratio (SNR) in wide-swath systems. The f-SCAN is an interesting alternative to the scanning on receive method (SCORE), which needs multiple phase centres achieved using the digital beam forming (DBF) technique. f-SCAN requires less hardware complexity than SCORE; at the same time, it improves the sidelobes and ambiguities’ suppression. The elements used in f-SCAN to generate the pencil beam are the true time delay lines (TTDLs) and the phase shifters (PSs). The general methodology to design an f-SCAN spaceborne SAR high-resolution wide-swath (HRWS) system is introduced; emphasis is put on the mathematical definition of the timing parameters and on a novel method of using TTDLs to achieve the full spanning of wide swaths. The processing of f-SCAN data is also considered: we introduce a novel algorithm to limit the data volume and to guarantee an almost invariant slant range impulse response function (IRF) by removing spectral distortions. Eventually, new definitions, specific for f-SCAN, of the well-known SAR performance parameters, are provided. Simulation results and performances are presented. The advantages and disadvantages with respect to SCORE are discussed using the design of a real case system.

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 m2 resolution in the 100 km of the swath of the ERS-1 satellite (launched in 1991) to the 5 × 20 m2 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, N c h 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 N c h 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]:
ϕ = 2 π f c sin θ · Δ y
where f is the frequency of the signal, Δ y 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, x ( t ) = s i n ( 2 π f 0 t ) , applying a delay or a phase produces the same effect: x ( t t 0 ) = s i n ( 2 π f 0 ( t t 0 ) ) and x ( t t 0 ) = s i n ( 2 π f 0 t + ϕ ) , where ϕ = 2 π f 0 t 0 . 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]:
ϕ = 2 π f Δ τ 2 π f c sin θ · Δ y
If the previous relation is written as a function of the scan direction θ [17]:
sin θ = c Δ y Δ τ ϕ 2 π f
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 θ 0 and width Δ θ = θ f a r , s w a t h θ n e a r , s w a t h 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:
  • L e l is the elevation size (i.e., height) of the whole array antenna.
  • N is the number of elements along the elevation.
  • N T D L is the number of TTDLs. Then, the N elements from the whole antenna are grouped into N T D L sub-arrays of N / N T D L elements each.
According to the first point of the strategy, defined at the end of the previous section, the phase shift ϕ p s was designed to provide an illumination towards the swath centre θ 0 at the carrier frequency f c . The relation between ϕ p s and θ 0 , already given in Equation (1), is here repeated [16]:
ϕ p s = 2 π f c c sin θ 0 · Δ y
Due to the finite bandwidth of the signal, B c h , the corresponding frequency dispersion is given by inverting Equation (4) with the chirp frequency edges:
Δ θ p s = sin 1 f c f c B c h / 2 · sin θ 0 sin 1 f c f c + B c h / 2 · sin θ 0
Usually, the angle interval Δ θ p s 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]:
F θ = sin π L e l λ sin θ sin π L e l / N T D L λ sin θ θ g r l o b e , 1 = sin 1 λ c L e l / N T D L λ c L e l / N T D L
where we recall that L e l is the antenna height and N T D L is the number of available delay lines.
The frequency dispersion at the k-th grating lobe is then:
Δ θ g r l o b e , k = sin 1 k · λ max L e l / N T D L sin 1 k · λ min L e l / N T D L k · Δ λ L e l / N T D L
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):
k r o u n d Δ θ Δ θ p s Δ θ g r l o b e , 1
Eventually, the delay for each line is derived by combining Equations (3), (6), and (8):
Δ τ = Δ Y c k · λ c L e l / N T D L
In Equations (9), Δ Y 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]:
Θ e l ( ϕ ) = Θ e m b ( ϕ ) m = 0 K 1 exp j 2 π f c Δ τ · m l = 0 N / K 1 exp j 2 π λ ( N K m + l ) Δ Y k · sin ϕ · exp j ( N K m + l ) ϕ p s
In the previous equation, K = N T D L , Θ e m b ( ϕ ) is the common embedded pattern from each radiating element, while λ = c / f 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., N T D L decreases), then:
  • The angular position of the grating lobes ( θ g r l o b e , 1 ) 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 T c h and total bandwidth B c h . 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: k c h = B c h / T c h . 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 k f s c a n and duration T f s c a n . Its width, B 0 , 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 B c h 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 k c h and the steering rate k f s c a n 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 R 1 , 2 , the slant ranges at the near and far range, the corresponding delays are: τ 1 , 2 = 2 R 1 , 2 / c . The geometric sampling window start time ( S W S T g e o ) depends on τ 1 (it is the fractional part of τ 1 once the integer numbers of P R I = 1 / P R F are removed), and the geometric sampling window end time ( S W E T g e o ) depends on τ 2 , leading to the definition of the geometric sampling window length ( S W L g e o ):
S W L g e o = S W E T g e o S W S T g e o
To allow the correct focusing of the very-far-range target, the geometric SWL ( S W L g e o , Equation (11)) of a conventional instrument must be extended by the duration of one chirp:
S W L i n s t r = S W L g e o + T c h
Thanks to the filtering effect of the f-SCAN beam (see Figure 4), the sampling window length of the f-SCAN system ( S W L f s c a n ) is shorter than what is stated in Equation (12):
S W S T f s c a n = S W S T g e o + B c h B k c h S W E T f s c a n = S W E T g e o + T c h B c h B k c h S W L f s c a n = S W L i n s t r 2 B c h B k c h
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 T i n t 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 T i n t is the integration time, which corresponds to the resolution band B:
T i n t = B k c h
Eventually, the steering rate k f s c a n 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 k f s c a n is:
k f s c a n = B c h B T f s c a n
where the scanning time T f s c a n is:
T f s c a n = S W L f s c a n T i n t
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: S W L f s c a n × f s , where f s > B c h 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 B c h 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 B 0 in Figure 4, rather than by B. B 0 depends not only on the resolution bandwidth B, but also on both the chirp rate k c h and steering rate k f s c a n . From the geometric considerations on Figure 4 and from Equation (14), we can write:
T i n t = B 0 B | k f s c a n | T i n t = B | k c h |
Then, by setting equal the right sides of Equation (17), we obtain the relation between B 0 and B:
B 0 = | k f s c a n | + | k c h | | k c h | · B
We define as the shrinking factor of the f-SCAN acquisition the quantity:
k s h r i n k = | k c h | | k f s c a n | + | k c h |
Equation (19) states that the target, because of the steering rate k f s c a n 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, B 0 is always greater than the band B necessary to fulfil the requirement of the geometric resolution. Eventually, we can take some margin α on B 0 to have the minimum sampling frequency f s :
f s = α · B 0 < < α · B c h
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 f s . 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 f s 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 N m o s a i k copies of the spectrum are stacked. Note that the mosaicking operation is equivalent, in the time domain, to the zero interleaving of N m o s a i k samples. According to Figure 5, the number of spectral supports mosaicked is:
N m o s a i k = c e i l B c h B 0
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 k f s c a n :
h d e r a m p t = exp j π k f s c a n · t 2
with t [ S W L f s c a n / 2 , S W L f s c a n / 2 ] . 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 N m o s a i k . Please note that: (i) the nominal bandwidth of the LPF is B 0 , 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 S W L f s c a n given by Equation (13) to S W L i n s t r given by Equation (12). This is obtained by padding the data with N 0 samples at the beginning and at the end, where N 0 , according to Equation (13), is defined as:
N 0 = B c h B B c h / T c h · f s , e x t
Please note that, after frequency spectrum restoration, f s , e x t = N m o s a i k · f s 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 ( t a c q ) with the time of the range-focused data ( t f o c ). 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, t f o c ). Hence, we can derive the mathematical relation between t f o c and t a c q by defining the quantity f a (highlighted in magenta in Figure 7) as:
f a = t a c q · | k f s c a n | + B 0 2 f a = B c h 2 t a c q t f o c · | k c h |
Then, by making the right sides of Equation (24) equal, we obtain the relation between t a c q and t f o c :
t a c q = | k c h | | k c h | + | k f s c a n | · t f o c + 1 2 B c h B 0 | k c h | + | k f s c a n |
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 t f o c . The central frequency f r g of the target at time t f o c is computed taking advantage of Equation (25) and k f s c a n as:
f r g = | k f s c a n | t a c q = | k c h | | k f s c a n | | k c h | + | k f s c a n | · t f o c + B c h B 0 2 | k f s c a n | | k c h | + | k f s c a n |
Thus, from the first term of Equation (26), we obtain the expression of the rate k f o c , which is the frequency modulation rate of the range-compressed data and, indeed, corresponds to the slope of the yellowish shape in Figure 7:
k f o c = | k c h | | k f s c a n | | k c h | + | k f s c a n |
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 k f o c , where k f o c 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 k c h 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.

5. FSCAN Performances

In this section, the typical SAR L1 performance figures are theoretically extended to the case of f-SCAN. Since f-SCAN is a method related to the elevation direction, here, we focused our attention only on the elevation-dependent parameters, i.e., the noise equivalent sigma zero (NESZ) and range-distributed target-to-ambiguity ratio (RgDTAR).

5.1. NESZ and RgDTAR Models

The traditional equation for the noise sensitivity (NESZ) [30] in the case of an f-SCAN system is modified as follows:
N E S Z i = 4 4 π 3 R θ i 3 · k T B F · V s s i n θ i · L S C · L A z · L P T X · d c e f f · D e l , 2 w , i θ i · λ i 3 · c N E S Z = Δ θ / 2 Δ θ / 2 N E S Z i · δ θ 0 θ i d θ i
In Equation (28), R θ i is the slant range at the given off-nadir angle θ i ; T is the temperature of noise; k is the Boltzmann’s constant; F is the noise figure of the receiver; V s is the orbital velocity of the satellite; P T X is the peak transmitted power; d c e f f is the effective transmission duty cycle (expressed as a fraction of one), which accounts only for the part of the chirp effectively used for range compression, i.e., d c e f f = d c · B B c h ; D e l , 2 w , i θ i is the two-way antenna directivity at the given angle θ i ; λ i is the wavelength corresponding to the instrument frequency for the i-th direction; c is the light speed. The range bandwidth B is the one that drives the resolution, according to the definition in Section 2.2.
As concerns the losses, there are the traditional contributions such as the integration loss, L A z , and all the other well-known terms summarized in L, which encompasses, e.g., TX and RX instrument losses, propagation losses, etc.
The f-SCAN system, however, experiences another kind of loss, which accounts for the effect of range spectral shaping due to the scanning pencil beam. This loss, which resembles the concept of the pulse extension loss (PEL) in SCORE [31], can be named the pulse loss, and it is addressed by L S C in Equation (28). It can be found by the following expression:
L S C = B / 2 B / 2 | W r g f | 2 d f B / 2 B / 2 | W r g f | 2 · D e l , 2 w f d f
where the term W r g f represents the contribution of the antenna whitening and windowing, if applied.
For the RgDTAR [6], the conventional formulation can be extended to f-SCAN with the same approach used for NESZ in Equation (28):
R g D T A R i = n = N N D e l , 2 w , i θ R i + n · c 2 · P R F · σ 0 θ R i + n · c 2 · P R F s i n θ i n c R i + n · c 2 · P R F · R i + n · c 2 · P R F 3 D e l , 2 w , i θ R i · σ 0 θ R i s i n θ i n c R i · R i 3 R g D T A R = Δ θ / 2 Δ θ / 2 R g D T A R i · δ θ 0 θ i d θ i
In the previous Equation (30), σ 0 θ i n c is the reflectivity at a given incidence angle. The other terms have the same meaning as in Equation (28). In particular, the terms R i + n · c / 2 · P R F are the range ambiguous positions.

5.2. Comparison with No-DBF and SCORE SAR

In order to quantify the capability of the f-SCAN system, here, we compared the system defined in Section 4 with two equivalent systems: (i) a traditional SAR acquiring in stripmap mode; (ii) a SAR using digital channels in elevation to obtain SCORE. In the two cases, the chirp bandwidth was set to B = 304 MHz (to achieve the same range resolution), and the swath extension was as in the f-SCAN simulation, i.e., about 44 km. Moreover, in the case of SCORE, eight digital channels were used. In Figure 12, the comparison between the NESZ and the RgDTAR of the three systems is provided. We can clearly see that the f-SCAN outperformed the other two systems both on NESZ and the RgDTAR.

5.3. Discussion on f-SCAN’s Advantages and Drawbacks

The main advantages of f-SCAN with respect to SCORE are:
  • A simpler onboard implementation. Only one receiver channel, rather than N, is required. Moreover, there is no need to have an onboard DBF recombination matrix to combine the N channels before downloading the data to the ground as the scanning of the pencil beam is achieved in the analogue domain. That means less mass (because of less hardware), volume, and power consumption.
  • The chance of using the scanning high-directivity beam in both transmit and receive. Hence, the NESZ of an f-SCAN system is expected to be lower than for a SCORE system (see Figure 12).
On the other hand, f-SCAN requires:
  • A more delicate design due to the close relation existing among directivity, the direction of pointing, and chirp frequency.
  • A wider pulse bandwidth and, in general, a higher duty cycle, which result in higher average onboard power compared to the usual systems and to SCORE.
Finally, also f-SCAN suffers from the same impairments as SCORE, i.e., the sensitivity to timing variation due to topography. A proper redesign of the timing would be necessary in that case, as it is being analysed in the progression of DBF missions such as ROSE-L and Sentinel-1NG.

6. Conclusions

In the paper, we addressed the f-SCAN technique, which is an alternative to SCORE to achieve almost constant high directivity throughout a wide swath. f-SCAN implements the scanning beam in the analogue domain by using a proper combination of phase shifters and true time delays. The role of phase shifters and TTDLs was discussed, and a novel application was presented, aimed at extending the angle interval swept by the phase shifter only. We showed how to reduce the data volume by an effective method that involves a pre-processing step such as in TOPS SAR. The relation between the fast time of the raw data and the fast time of the range-focused data was explained, which is not trivial. Eventually, new definitions for both NESZ and RgDATR were provided. The performance of an f-SCAN system was finally compared to the corresponding performances of a SCORE system and a conventional SAR, proving the potential benefits of the former.

Author Contributions

Conceptualization: P.G. and D.M.; Methodology: P.G., D.M. and A.R.P.; Software: P.G., D.M. and A.R.P.; Validation and overview: D.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Authors ensure that the data shared are in accordance with consent provided by participants on the use of confidential data.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In Section 2.2, we addressed the problem of designing the phase shifters and the delay lines for f-SCAN application. In particular, as concerns the dispersion contribution coming from the grating lobes, we approximated the factor k according to Equation (8). Here, we want to propose a method to have a better estimation of the parameter k. First of all, we point out two main aspects:
  • The angular interval to be scanned thanks to the frequency dispersion is usually larger than the swath width Δ θ as, according to Figure 8, left panel, the aim is to have the necessary bandwidth B, properly weighted by the elevation pencil beam, also for the targets at the very edges of the swath (both near and far range).
  • The scanning rate of the antenna pencil beam is indeed a further degree of freedom in the design of an f-SCAN instrument. Furthermore, k f s c a n , defined in Equation (15), does not depend on the antenna parameters (i.e., phase shifters and TTDLs). Then, a proper antenna design shall aim at maximizing the energy collected along the reddish shape in Figure 4.
To find out the best value k, we start by searching for the new angular interval that satisfies the first observation in the previous paragraph. We can write the following system of equations:
θ n e a r , b e a m = θ n e a r , s w a t h B 2 1 k f s c a n · Δ θ b e a m T c h θ f a r , b e a m = θ f a r , s w a t h + B 2 1 k f s c a n · Δ θ b e a m T c h
In the previous Equations (A1), Δ θ b e a m = θ f a r , b e a m θ n e a r , b e a m is the angular interval swept by the pencil beam, and then, the angle θ n e a r , b e a m is off-nadir angle, corresponding to the maximally steered (toward the nadir) pencil beam. On the contrary, with θ n e a r , s w a t h , we refer to the off-nadir angle corresponding to the near range of the acquired swath, which is different (in general: major) with respect to θ n e a r , b e a m . The ratio Δ θ b e a m T c h takes into account the angle spanned with the frequency of the chirp during the pulse time T c h . Obviously, similar notation stands also for the far-range angles. By solving Equation (A1) for the two unknowns θ n e a r , b e a m and θ f a r , b e a m , the solution is:
θ n e a r , b e a m = θ n e a r , s w a t h B k f s c a n 1 T c h · θ n e a r , s w a t h + θ f a r , s w a t h 2 1 B k f s c a n 1 T c h θ f a r , b e a m = θ f a r , s w a t h B k f s c a n 1 T c h · θ n e a r , s w a t h + θ f a r , s w a t h 2 1 B k f s c a n 1 T c h
Then, we can compute the optimal value of k by taking advantage of Equations (A2) and (6) (and using a more precise version of Equation (8)). Eventually, we can write:
k = a r g m i n k N { | Δ θ b e a m Δ θ p s sin 1 k · λ max L e l / N T D L sin 1 k · λ min L e l / N T D L | }
where the previous equation must be numerically solved with the constraint that k must be an integer.
To conclude, we point out that Equation (8) is derived from Equation (A3) by approximating the sin 1 · function with its argument. Obviously, the higher the value of k is, the less accurate the approximation is.

References

  1. Moreira, A.; Prats-Iraola, P.; Younis, M.; Krieger, G.; Hajnsek, I.; Papathanassiou, K.P. A tutorial on synthetic aperture radar. IEEE Geosci. Remote Sens. Mag. 2013, 1, 6–43. [Google Scholar] [CrossRef] [Green Version]
  2. Younis, M.; Fischer, C.; Wiesbeck, W. Digital beamforming in SAR systems. IEEE Trans. Geosci. Remote Sens. 2003, 41, 1735–1739. [Google Scholar] [CrossRef]
  3. Krieger, G.; Gebert, N.; Moreira, A. Unambiguous SAR signal reconstruction from nonuniform displaced phase center sampling. IEEE Geosci. Remote Sens. Lett. 2004, 1, 260–264. [Google Scholar] [CrossRef] [Green Version]
  4. Suess, M.; Grafmueller, B.; Zahn, R. A novel high resolution, wide swath SAR system. Scanning the Present and Resolving the Future. In Proceedings of the IEEE 2001 International Geoscience and Remote Sensing Symposium (Cat. No.01CH37217), Sydney, Australia, 9–13 July 2001; Volume 3, pp. 1013–1015. [Google Scholar] [CrossRef]
  5. Freeman, A.; Johnson, W.; Huneycutt, B.; Jordan, R.; Hensley, S.; Siqueira, P.; Curlander, J. The “Myth” of the minimum SAR antenna area constraint. IEEE Trans. Geosci. Remote Sens. 2000, 38, 320–324. [Google Scholar] [CrossRef] [Green Version]
  6. John, C. Curlander, R.N.M. Synthetic Aperture Radar: Systems and Signal Processing; Wiley: Hoboken, NJ, USA, 1992; p. 672. [Google Scholar]
  7. Davidson, M.W.J.; Furnell, R. ROSE-L: Copernicus L-Band Sar Mission. In Proceedings of the 2021 IEEE International Geoscience and Remote Sensing Symposium IGARSS, Brussels, Belgium, 11–16 July 2021; pp. 872–873. [Google Scholar]
  8. Geudtner, D.; Tossaint, M.; Davidson, M.; Torres, R. Copernicus Sentinel-1 Next Generation Mission. In Proceedings of the 2021 IEEE International Geoscience and Remote Sensing Symposium IGARSS, Brussels, Belgium, 11–16 July 2021; pp. 874–876. [Google Scholar] [CrossRef]
  9. Kim, J.H.; Younis, M.; Prats-Iraola, P.; Gabele, M.; Krieger, G. First Spaceborne Demonstration of Digital Beamforming for Azimuth Ambiguity Suppression. IEEE Trans. Geosci. Remote Sens. 2013, 51, 579–590. [Google Scholar] [CrossRef] [Green Version]
  10. Gebert, N.; Krieger, G.; Moreira, A. Multichannel Azimuth Processing in ScanSAR and TOPS Mode Operation. IEEE Trans. Geosci. Remote Sens. 2010, 48, 2994–3008. [Google Scholar] [CrossRef] [Green Version]
  11. Giudici, D.; Guccione, P.; Manzoni, M.; Guarnieri, A.M.; Rocca, F. Compact and Free-Floating Satellite MIMO SAR Formations. IEEE Trans. Geosci. Remote Sens. 2022, 60, 1–12. [Google Scholar] [CrossRef]
  12. Guccione, P.; Monti Guarnieri, A.; Rocca, F.; Giudici, D.; Gebert, N. Along-Track Multistatic Synthetic Aperture Radar Formations of Minisatellites. Remote Sens. 2020, 12, 124. [Google Scholar] [CrossRef]
  13. Tan, W.; Xu, W.; Huang, P.; Huang, Z.; Qi, Y.; Han, K. Investigation of Azimuth Multichannel Reconstruction for Moving Targets in High Resolution Wide Swath SAR. Sensors 2017, 17, 1270. [Google Scholar] [CrossRef] [Green Version]
  14. Freeman, A.; Krieger, G.; Rosen, P.; Younis, M.; Johnson, W.T.K.; Huber, S.; Jordan, R.; Moreira, A. SweepSAR: Beam-forming on receive using a reflector-phased array feed combination for spaceborne SAR. In Proceedings of the 2009 IEEE Radar Conference, Pasadena, CA, USA, 4–8 May 2009; pp. 1–9. [Google Scholar] [CrossRef]
  15. Gebert, N.; Krieger, G.; Moreira, A. Digital Beamforming on Receive: Techniques and Optimization Strategies for High-Resolution Wide-Swath SAR Imaging. IEEE Trans. Aerosp. Electron. Syst. 2009, 45, 564–592. [Google Scholar] [CrossRef] [Green Version]
  16. Mailloux, R. Phased Array Antenna Handbook, 3rd ed.; Artech House: Norwood, MA, USA, 2017. [Google Scholar]
  17. Roemer, C. Introduction to a new wide area sar mode using the F-SCAN principle. In Proceedings of the 2017 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Fort Worth, TX, USA, 23–28 July 2017; pp. 3844–3847. [Google Scholar] [CrossRef]
  18. Roemer, C.; Gierlich, R.; Marquez-Martinez, J.; Notter, M. Frequency Scanning applied to Wide Area SAR Imaging. In Proceedings of the 12th European Conference on Synthetic Aperture Radar (EUSAR 2018), Aachen, Germany, 4–7 June 2018; pp. 1–5. [Google Scholar]
  19. Zhao, Q.; Zhang, Y.; Wang, W.; Liu, K.; Deng, Y.; Zhang, H.; Wang, Y.; Zhou, Y.; Wang, R. On the Frequency Dispersion in DBF SAR and Digital Scalloped Beamforming. IEEE Trans. Geosci. Remote Sens. 2020, 58, 3619–3632. [Google Scholar] [CrossRef]
  20. Chen, Z.; Zhang, Z.; Wang, W.; Zhao, Q.; Wen, Y.; Meng, X. Performance Demonstration of Dispersive SCORE: Digital Scalloped Beamforming With X-Band and C-Band DBF-SARs. IEEE Geosci. Remote Sens. Lett. 2022, 19, 1–5. [Google Scholar] [CrossRef]
  21. Xu, W.; Yu, Q.; Fang, C.; Huang, P.; Tan, W.; Qi, Y. Onboard Digital Beamformer with Multi-Frequency and Multi-Group Time Delays for High-Resolution Wide-Swath SAR. Remote Sens. 2021, 13, 4354. [Google Scholar] [CrossRef]
  22. Bordoni, F.; Krieger, G.; Younis, M. Multifrequency Subpulse SAR: Exploiting Chirp Bandwidth for an Increased Coverage. IEEE Geosci. Remote Sens. Lett. 2019, 16, 40–44. [Google Scholar] [CrossRef] [Green Version]
  23. Alvarez, Y.; García González, C.; Vazquez-Antuna, C.; Hoeye, S.; Las-Heras, F. Frequency Scanning Based Radar System. Prog. Electromagn. Res. 2012, 132, 275–296. [Google Scholar] [CrossRef] [Green Version]
  24. Safaai-Jazi, A.; Haupt, R.L.; Nayeri, P. Pulse Dispersion in Phased Arrays. Int. J. Antennas Propag. 2017, 2017, 275–296. [Google Scholar] [CrossRef]
  25. Scheiber, R.; Martone, M.; Gollin, N. Chirp Selection and Data Compression for Spaceborne Wide-Swath SAR in FScan-Mode. In Proceedings of the 13th European Conference on Synthetic Aperture Radar (EUSAR 2021), Online, 29 March–1 April 2021; pp. 1–6. [Google Scholar]
  26. Nan, L.; Gai, G.; Shiyang, T.; Linrang, Z. Signal Modeling and Analysis for Elevation Frequency Scanning HRWS SAR. IEEE Trans. Geosci. Remote Sens. 2020, 58, 6434–6450. [Google Scholar] [CrossRef]
  27. Gatelli, F.; Monti Guamieri, A.; Parizzi, F.; Pasquali, P.; Prati, C.; Rocca, F. The wavenumber shift in SAR interferometry. IEEE Trans. Geosci. Remote Sens. 1994, 32, 855–865. [Google Scholar] [CrossRef] [Green Version]
  28. De Zan, F.; Monti Guarnieri, A. TOPSAR: Terrain Observation by Progressive Scans. IEEE Trans. Geosci. Remote Sens. 2006, 44, 2352–2360. [Google Scholar] [CrossRef]
  29. Cafforio, C.; Prati, C.; Rocca, F. SAR data focusing using seismic migration techniques. IEEE Trans. Aerosp. Electron. Syst. 1991, 27, 194–207. [Google Scholar] [CrossRef]
  30. Calabrese, D.; Episcopo, R. Derivation of the SAR Noise Equivalent Sigma Nought. In Proceedings of the 10th European Conference on Synthetic Aperture Radar (EUSAR 2014), Berlin, Germany, 3–5 June 2014; pp. 1–4. [Google Scholar]
  31. Younis, M.; Rommel, T.; Bordoni, F.; Krieger, G.; Moreira, A. On the Pulse Extension Loss in Digital Beamforming SAR. IEEE Geosci. Remote Sens. Lett. 2015, 12, 1436–1440. [Google Scholar] [CrossRef]
Figure 1. Antenna patterns of a frequency scanning array using phase shifters only. The figure was generated using the parameters presented in Section 4.
Figure 1. Antenna patterns of a frequency scanning array using phase shifters only. The figure was generated using the parameters presented in Section 4.
Remotesensing 14 05283 g001
Figure 2. Angular axis stretching effect, from TTDLs only, due to wide-pulse bandwidth effect. Left panel: k = 0 (see Equation (7)). Right panel: k = 5 (see Equation (7)). The figures were generated using the parameters presented in Section 4.
Figure 2. Angular axis stretching effect, from TTDLs only, due to wide-pulse bandwidth effect. Left panel: k = 0 (see Equation (7)). Right panel: k = 5 (see Equation (7)). The figures were generated using the parameters presented in Section 4.
Remotesensing 14 05283 g002
Figure 3. A number of pencil beams at the right swath edge are shown. Solid lines: 8 TTDLs used; dashed lines: 4 TTDLs used. The colour code of the beams addresses the relative frequency: from green to red, the frequency increases. The figure was generated using the parameters presented in Section 4.
Figure 3. A number of pencil beams at the right swath edge are shown. Solid lines: 8 TTDLs used; dashed lines: 4 TTDLs used. The colour code of the beams addresses the relative frequency: from green to red, the frequency increases. The figure was generated using the parameters presented in Section 4.
Remotesensing 14 05283 g003
Figure 4. Fast time–frequency diagram of raw SAR data. Quantities in the figure are described in Equations (11)–(16).
Figure 4. Fast time–frequency diagram of raw SAR data. Quantities in the figure are described in Equations (11)–(16).
Remotesensing 14 05283 g004
Figure 5. Fast time–frequency diagram after SWL and f s reduction.
Figure 5. Fast time–frequency diagram after SWL and f s reduction.
Remotesensing 14 05283 g005
Figure 6. Block diagram of the unfolding procedure.
Figure 6. Block diagram of the unfolding procedure.
Remotesensing 14 05283 g006
Figure 7. Fast time–frequency diagram after the restoration of the full spectral support and geometrical receiving window length interval.
Figure 7. Fast time–frequency diagram after the restoration of the full spectral support and geometrical receiving window length interval.
Remotesensing 14 05283 g007
Figure 8. Each target experiences a different bandwidth. (Left panel) The spectral support of each simulated target in the swath. For reference, the nominal chirp spectrum (in black) is also drawn. (Right panel) Time–frequency diagram of the simulated point targets.
Figure 8. Each target experiences a different bandwidth. (Left panel) The spectral support of each simulated target in the swath. For reference, the nominal chirp spectrum (in black) is also drawn. (Right panel) Time–frequency diagram of the simulated point targets.
Remotesensing 14 05283 g008
Figure 9. Time–frequency diagram of the simulated point targets acquired using the subsampled range frequency.
Figure 9. Time–frequency diagram of the simulated point targets acquired using the subsampled range frequency.
Remotesensing 14 05283 g009
Figure 10. (Left panel) Range line with the 11 simulated point target after pre-processing and range compression. Peaks’ heights were normalized to 1. (Right panel) Zoom on the central point target aimed to show its IRF sinc-like shape.
Figure 10. (Left panel) Range line with the 11 simulated point target after pre-processing and range compression. Peaks’ heights were normalized to 1. (Right panel) Zoom on the central point target aimed to show its IRF sinc-like shape.
Remotesensing 14 05283 g010
Figure 11. Top: ground range resolution; bottom, left-axis: PSLR; right-axis: ISLR.
Figure 11. Top: ground range resolution; bottom, left-axis: PSLR; right-axis: ISLR.
Remotesensing 14 05283 g011
Figure 12. Compared performances for an f-SCAN (blue dotted line), SCORE (red dashed line), and standard (solid yellow line) SAR system. Left panel: NESZ; right panel: RgDTAR.
Figure 12. Compared performances for an f-SCAN (blue dotted line), SCORE (red dashed line), and standard (solid yellow line) SAR system. Left panel: NESZ; right panel: RgDTAR.
Remotesensing 14 05283 g012
Table 1. System parameters.
Table 1. System parameters.
ParameterValueUoM
Central frequency9.8GHz
Available bandwidth1.2GHz
Antenna height1.5m
No. phase centre64
No. time delay line8
Mechanical boresight30deg
PRF2560Hz
Duty cycle15%
Table 2. Geometric parameters.
Table 2. Geometric parameters.
ParameterValueUoM
Incidence angle interval[21.35–25.95]deg
Satellite height510km
Off-nadir angle interval[19.70–23.90]deg
Slant range extension17.77km
Ground range extension44.28km
Ground range resolution1.2m
Table 3. System timing.
Table 3. System timing.
ParameterSymbolValueUoM
Geometric SWL S W L g e o 118.56 μ s
Instrument SWL S W L i n s t r 177.15 μ s
Chirp duration T c h 58.59 μ s
Resolution time T i n t 14.84 μ s
Scanning time T f s c a n 74.81 μ s
Scanning SWL S W L f s c a n 89.65 μ s
Resolution bandwidthB304MHz
Instantaneous bandwidth B 0 481.80MHz
Chirp rate k c h −20.48MHz/ μ s
fSCAN rate k f s c a n 11.98MHz/ μ s
Shrink factor α 0.631
True time delay Δ τ 0.410ns
Phase shifter ϕ p s −39.34deg
Table 4. Pulse timing.
Table 4. Pulse timing.
ParameterValue [ μ s]
TX pulse start (ref)0
TX pulse end58.59
RX pulse start160.72
RX pulse end191.77
Start of next pulse (PRI)390.62
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Guccione, P.; Mapelli, D.; Giudici, D.; Persico, A.R. Design of f-SCAN Acquisition Mode for Synthetic Aperture Radar. Remote Sens. 2022, 14, 5283. https://doi.org/10.3390/rs14205283

AMA Style

Guccione P, Mapelli D, Giudici D, Persico AR. Design of f-SCAN Acquisition Mode for Synthetic Aperture Radar. Remote Sensing. 2022; 14(20):5283. https://doi.org/10.3390/rs14205283

Chicago/Turabian Style

Guccione, Pietro, Daniele Mapelli, Davide Giudici, and Adriano Rosario Persico. 2022. "Design of f-SCAN Acquisition Mode for Synthetic Aperture Radar" Remote Sensing 14, no. 20: 5283. https://doi.org/10.3390/rs14205283

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop