PRF Selection in Formation-Flying SAR: Experimental Verification on Sentinel-1 Monostatic Repeat-Pass Data

: Formation-flying synthetic aperture radar (FF-SAR) enables new working modes and can achieve very high performance through a series of very compact, low-weight, satellite platforms thanks to passive operations of conveniently distributed formation-flying receivers. System timing is a crucial aspect of FF-SAR. The manuscript presents a novel approach to pulse repetition frequency (PRF) selection in order to obtain a uniform distribution of samples at given platform positions. A digital beamforming algorithm is applied on a stack of monostatic repeat-pass images collected by the Sentinel-1 system to test the validity of the PRF selection method. Processed images were thus properly selected to achieve the best merit index measuring the quality of samples distribution. The results show that: (a) the image resulting from beamforming features better azimuth ambiguity-to-signal ratio and (b) the proposed approach for PRF selection allows one to individuate a subset of the available images leading to uniform distribution of samples which can be used to support FF-SAR processing.


Introduction
A formation-flying synthetic aperture radar (FF-SAR) is a spaceborne distributed synthetic aperture radar [1][2][3] in which the signal emitted by the transmitter and scattered from the area of interest is not collected by a single receiver but by many conveniently distributed formation-flying receivers [4]. FF-SAR is a semiactive multistatic system [5] composed by a single active illuminator and a number of receiving-only passive systems. Bistatic raw data are thus collected by each receiver: these data can be coherently combined to generate either a single, higher performance image [6][7][8][9][10] or to implement single-pass interferometric SAR techniques [5,11]. This paper investigates the first topic and specifically analyzes the role played by the along-track separations among the receivers. An FF-SAR in which N receivers fly with an along-track baseline is a sensor receiving N pulses [7,12] or samples at azimuth-displaced positions for each transmitted pulse. The factor N can be interpreted as the redundancy, or the maximum number of degrees of freedom [13] of the distributed system when compared to a monostatic SAR. Again, such a redundancy can be exploited in different ways [13][14][15][16][17], from Ground Moving Target Indication (GMTI) to High-Resolution Wide-Swath (HRWS) imaging, pending the application of suitable reconstruction algorithms [7][8][9]13] and the availability of a sufficient number of receivers.
It is evident that the concept shows similarities with the along-track multichannel SAR [7,13], but the same important differences arise. Specifically, along-track separations among the receivers of an FF-SAR are longer and change either along the orbit or during different phases of the mission. Moreover, cross-track and vertical baselines are established due to relative orbit dynamics and for reducing the collision risk [18,19]. System timing is thus more intricate than for a multichannel SAR [20], and proper selection of the pulse repetition frequency (PRF) is mandatory. This work firstly contributes to that topic presenting an approach to PRF selection.
A key feature of FF-SAR is the capability to synthesize an image with a higher PRF with respect to the actual PRF of the transmitter. Different techniques for image enhancement result from this capability like coherent resolution improvement, swath width widening, and reduction of the ambiguity levels. PRF can be increased when the collected samples are not overlapped [12,13] and relatively distant to each other. The position of the samples within the synthetic aperture depends on the along-track separations among the receivers but also on the actual transmitter PRF. Indeed, due to errors in the initial conditions and orbit perturbations, satellite separations will drift from the nominal values. Unfortunately, slight variations in the along-track baselines have a strong impact on the distribution of samples, thus asking for a very tight relative orbit control (i.e., at sub-wavelength level). Proper PRF selection is a way to achieve the desired distribution of samples while relaxing the requirement of relative orbit control. Specifically, based on the current relative separations, a suitable PRF value can be computed and used by the transmitter. The approach allows the focus to be moved from tight orbit control to accurate relative navigation, which is easier to be achieved, especially for small satellites.
The synthesis of an enhanced image from the bistatic raw data collected by each receiver is obtained by the application of a reconstruction or beamforming algorithm. The term reconstruction refers to the capability to recover an unambiguous Doppler spectrum from aliased azimuth signals (i.e., signals sampled at a frequency lower than Nyquist rate). Depending on the residual number of degrees of freedom, more general processing methods [13] can be introduced to obtain other features than the mere reconstruction, for example, preservation of the Signal-to-Noise Ratio (SNR). The term beamforming is instead related to the possibility of considering the physical array realized by the relative positions among the receivers and of applying beamforming strategies to that array to counteract the effect of aliasing (e.g., to perform null-steering toward the direction of azimuth ambiguities). As noted in [7], these two approaches are, to a given extent, equivalent because they represent two alternative ways to formulate the same problem. An example of reconstruction algorithm is that presented in [12] and further analyzed in [7]. The algorithm was specifically developed for a multichannel, single-platform SAR and then tailored to a multiplatform, distributed SAR by introducing a corrective factor to locate the samples properly [12]. The antenna patterns were not considered in [12] as they were assumed to be equal. Signal processing was generalized in [13] by accounting for the more realistic case in which i) antenna patterns can differ and ii) the actual system redundancy is smaller than the number of available receivers. Experimental demonstrations of reconstruction performance of the different methods were also carried on quad-channel Radarsat-2 data [13]. When an FF-SAR is considered, with both along-track and cross-track/vertical baselines, the reconstruction algorithm of [7] is no longer valid because further phase and frequency contributions are generated in that case. In practice, the approach of [7] could be applied for receivers characterized by very limited (tens of meters), along-track-only separations. For coping with longer baseline formations, the idea has been proposed [8,9] to apply the reconstruction algorithm after range compression and, more importantly, following a phase compensation step in which, for instance, the topographic phase, generated by cross-track/vertical baselines, is removed. In addition, a generalized version of the reconstruction algorithm is applied to deal with the case of longer alongtrack baselines to compensate for Doppler centroid and Doppler rate variations among different receivers. The approach, already tested on TanDEM-X data, is able to deal with much longer baselines. Alternatively, the impact of range cell migrations among different receivers was addressed in [17] performing signal reconstruction in the wavenumber domain.
Because of the unavailability of bistatic raw, or range compressed, data collected by more than two formation-flying receivers, a beamforming algorithm is considered in this paper to verify the proposed approach for PRF selection. The algorithm, based on [6], works in the space-time domain and processes focused bistatic images exploiting the principle of convolution between the abovementioned physical array, realized by the positions of the receivers, and the bistatic images [10]. The added value of this choice is that the verification can be performed even on multimonostatic data (i.e., collected in a repeat-pass scenario). In detail, from the point of view of the spatial distribution of receivers, a standard repeat-pass interferometric data stack is equivalent (up to a factor of two in the phase terms) to a single-pass FF-SAR configuration, so the same processing algorithm can be applied. However, if from one side this geometric equivalence is verified, from the other side a repeatpass data stack shows significant time decorrelation which makes the interpretation of the obtained results difficult. For this reason, the analysis is conducted in this paper over selected portions of the images showing high SNR and enough stability throughout the data stack, such as a strong manmade structure close to the shoreline and the relevant azimuth ambiguity which is a very stable feature, easily identifiable over the darker sea background.
The paper is organized as follows. Section 2 discusses constraints and requirements of PRF selection for FF-SAR in different working modes and presents the PRF selection method. Section 3 details the considered beamforming approach. Then, Sections 4 and 5 apply the beamforming procedure to images gathered by the Sentinel-1 mission.

PRF Selection for FF-SAR
Proper PRF selection is crucial to enable FF-SAR operation. If one considers the Doppler bandwidth, BD, of the scene illuminated by the transmitter, and an FF-SAR including N receivers, two main different cases must be discussed: Equation (1) implies that the signal collected by each receiver undersamples BD, leading to signal aliasing and azimuth ambiguities. In such a case, depending on the considered satellite formation, the reconstruction algorithm presented in [12] can be used to combine N received signals for generating a single azimuth signal sampled at N·PRF enabling the suppression of azimuth ambiguities. In an ideal, noise-free, and clutter-free scenario, the result of these algorithms is not affected by PRF, as long as samples of different receivers do not coincide. However, as discussed in [12,13,16,21], reconstruction of azimuth signals from a periodic nonuniform sampling can amplify, notably, either side-lobe clutter or noise level. For reducing such an effect, the distribution of samples must be as uniform as possible.
According to [7,12], once having set the PRF (e.g., ⋅ ≅ D N PRF B ), Equation (3) can be used to calculate the exact relative positions among the receivers able to realize a uniform distribution of samples.
where Δ j x is the along-track separation between the transmitter and the j th receiver, s v is platform velocity, and j k is an integer number that is, in general, different for each receiver.
The application of these models to an FF-SAR involves tight constraints on the control of relative positions among the formation-flying satellites [15]. However, the problem can be formulated in an inverse fashion. Specifically, based on satellite orbit dynamics, satellites can be expected to fly with relatively irregular separations. This means that PRF can be used as the only available (and to a certain extent) free parameter that must be tuned or adjusted to guarantee that collected samples show the highest degree of uniformity which is achievable by the actual spatial distribution of receivers. The basic assumption for making this approach feasible is the accurate knowledge of relative positions among the receivers.
An algorithm to perform PRF selection can be defined as follows. One can assume a linear FF-SAR formation with N satellites, flying at constant velocity vs, where the satellite number N is the leading satellite of the formation and the satellite 1 is the last one in the along-track direction. One can also assume that satellite 1 is also the transmitting one. At a given time instant, t = t0, relative platform positions are known, so the locations of samples x10, x20, …xN0, at this time instant are known too. With specific reference to receiver N, it is important to focus on the sample collected at t = t0 (i.e., xN0) and on the next collected sample by the same receiver. This sample is collected at t = t0 + 1/PRF and it is located at xN0 + vs/PRF, for any given value of PRF. The goal of PRF selection is to find the repetition frequency able to fill, almost uniformly, the distance from xN0 to xN0 + vs/PRF, using samples collected by the other receivers at different time epochs, that is, corresponding to different transmitted pulses (see Figure 1). If the generic receiver j is considered, at a time epoch, t = t0 + κj/PRF, the sample collected by this receiver will be located at The integer κj, which is different from kj in Equation (3), is to be evaluated. The sample collected by receiver j after κ j pulse repetition interval (PRI) from t0, is located at: Based on Equation (4) one can write so, the unknown integer parameter can be estimated as indicates truncation to integer. If one computes κ j for each receiver, the sequence of samples and the relevant distribution can be computed for any value of PRF. It is worth noting that the sequence of samples can change with the varying PRF. For instance, assuming N = 3, two sequences of samples are possible, namely, 3-1-2-3 and 3-2-1-3.
As an example, the case of three receivers is considered as illustrated in Figure 1, which shows the positions of samples (upper line) and platform positions as a function of time (lower lines).  x 20 x 30 The spacing among the receivers is not uniform, though it is not very irregular as expected in an FF-SAR spaceborne mission. In detail, the transmitter platform is also receiving, hence Tx and Rx1 positions coincide, whereas the receivers 2 and 3 are dislocated. In Figure 1, different time instants have been shown. The first is after 1/PRF from t0. Indeed, a distance of vs/PRF exists between sample 1 at t0 (shown in red) and the one gathered after 1/PRF (shown in violet). The second time instant (t = t0 + κ2/PRF) is when sample 2 is between x30 and sample 3 gathered after 1/PRF (sample 3 in violet). Similarly, the third time (t = t0+ κ1/PRF) is when sample 1 is between x30 and sample 3 gathered after 1/PRF (sample 3 in violet).
The best PRF, that is, the PRF increasing the level of uniformity of samples, can be thus individuated considering the distance n d between two consecutive samples within the interval of Equation (4) and introducing the J-index as which measures the deviation of n d from the case of uniform distribution of samples. According to the definition of J, the index ranges from 0 to 1 with the minimum value achieved for uniform distribution and the maximum value obtained for the case in which all the samples are overlapped.
Assuming the parameters of Table 1, the trend of J as a function of PRF is presented in Figure  2a. Specifically, PRF is increased up to 5% of the nominal value of 3628.4 Hz. J achieves both maxima and minima within the simulated range of PRF. Figure 2b shows the distribution of samples corresponding to some remarkable points along the curve of the index J. Starting from the first point A, both sample 2 and sample 1 drift leftward so the level of nonuniformity increases. J-index increases as well. When sample 2 reaches sample 3, the condition of point B is obtained. From this case, J-index continues to increase with the increasing PRF, up to the condition in which both 2 and 1 are very close to 3, thus representing a local maximum for J. If PRF is further increased, J decreases because 1 and 2 drift leftwards but with different drifting rates with respect to PRF. The result is that samples tend to become more and more uniform, up to the condition of a local minimum of J in which sample 1 is at one-third of the distance of vs/PRF from the left sample 3 and sample 2 is very close to a distance of one-third of vs/PRF from the right sample 3 (this corresponds to point E in Figure 2b). In this case, samples are very close to the condition of uniform distribution. Any PRF increase from this condition leads to an increase of J-index, which reaches a local maximum (point F) when sample 1 is overlapped to sample 3. Finally, the behavior of points G, H, and I can be explained using a similar argument as for points E, B, and C, respectively. In conclusion, based on the presented analysis, the best PRF should be selected as that corresponding to point E or to point G. Specifically, from the point of view of spatial sampling, those points are very similar.  As mentioned, the approach is feasible under the assumption of accurate knowledge of relative positions among the receivers. To gain further insight into this aspect, J-index is analyzed assuming an error of 10 cm (1 δ) for the relative positions among the receiver. This error matches the onboard realtime performance demonstrated by the PRISMA mission [22]. The selected PRF values, around the condition of the point E of Figure 2, for 100 realizations of the relative positioning errors, are shown in Figure 3. It is evident that when relative positioning errors are taken into account, the selected PRF is subject to very limited variations (i.e., order of 5 Hz) which do not affect the quality of the reconstructed signals. The proposed approach to PRF selection is valid for the case of Equation (1), that is, when FF-SAR is used to realize a higher PRF than that actually used. However, the approach can find application also in the case of Equation (2), when signals collected by the receivers are not aliased. Such a working condition can be implemented to perform GMTI. Specifically, it is well known that a very useful condition for successful GMTI application is that for which samples collected by different receivers are overlapped in space. This corresponds to the so-called Displaced Phase Center Antenna (DPCA) condition [14]. Under the DPCA condition, ground clutter for the stationary scene is common to all the collected signals so it can be canceled, thus allowing one to isolate moving targets and to estimate their velocity. As a consequence, for GMTI, in addition to Equation (2), PRF selection can be performed choosing the value of PRF that maximizes the index J for the current satellite configuration.
The concept of anti-DPCA condition is discussed in [13] with reference to a single-platform multichannel system. Specifically, the anti-DPCA condition is verified when samples from different receivers are perfectly interleaved: several well-distributed, isolated PRF values realize the anti-DPCA condition for a single-platform multichannel system. When an FF-SAR is considered, characterized by irregular and large physical separations among the receivers, neither DPCA nor anti-DPCA can be theoretically verified. Specifically, if one adjusts a single parameter (i.e., the PRF), it is not possible to obtain perfectly interleaved (or overlapping) samples. This is evident from Figure 4, in which the J-index is estimated for a multichannel (Figure 4a) and for an FF-SAR configuration (Figure 4c). In Figure 4b, it has been shown that the DPCA condition (point D) is perfectly realized for a given PRF value and the anti-DPCA condition (point C) occurs for a PRF value leading to a zero J-index. On the contrary, the trend of the J-index trend for an FF-SAR is more complex and many PRF values place the configuration near to DPCA and anti-DPCA condition. Thus, the proposed method measures the degree of uniformity of samples distribution and allows one to identify conditions that can resemble the anti-DPCA (or the DPCA) condition.  FF-SAR used to increase the PRF of the transmitter by a factor given by the total number of receivers.
Depending on system requirements and constraints, an FF-SAR with N receivers can be designed (i) to increase the PRF by a factor of K < N and (ii) to use residual degrees of freedom to satisfy requirements concerning further signal parameters, like SNR, resolution, Peak-to-Side-Lobe Ratio (PSLR), Point Target Ambiguity Ratio (PTAR), and Signal-to-Ambiguity-plus-Noise Ratio (SANR). The reconstruction algorithm of [7] cannot be applied in this case and more general processing methods must be used. The problem is addressed in [13] thoroughly. When the PRF of the transmitter has to be enhanced by a factor of K, it is still very important for K out of N samples to be nearly equispaced. This is, for instance, the case of the generalized matrix inversion method [13] and of the projection method [16] in which, even if all the N received signals are processed by the reconstruction algorithm, only K samples, realizing the most uniform samples distribution, must be selected to define the kernel of the algorithm. With reference to these more general processing methods, the proposed approach to PRF selection is still useful. Namely, one can compute the J-index for a range of PRF values which is slightly higher than 1/K-th of the Nyquist rate (e.g., 5%-10%) and considering all the possible configurations of K receivers out of the available N ones. The minimum J-index thus corresponds to the best degree of uniformity achieved by K out of N receivers (Section 4 will come back to this topic).

Beamforming Algorithm for FF-SAR
Beamforming is a processing approach that can be used to coherently combine the data collected by each receiver of an FF-SAR into a higher performance image or product.
An example of digital beamforming is the one used to suppress azimuth ambiguities generated when the working PRF is lower than the Doppler bandwidth of the illuminated scene. This principle is herein applied to a formation-flying satellite geometry that is accounting for velocity variations among the receivers and nonlinear motion of the platforms, together with Earth surface curvature.
Considering an FF-SAR formation of N receivers, and following bistatic SAR focusing for each Tx/Rx couple, N focused bistatic images, In with n = 1, …, N, are available. Without loss of generality, it is assumed herein that (i) the transmitter is also a receiver, (ii) bistatic images have been coregistered and georeferenced to a common image grid, and (iii) the first forward and backward azimuth ambiguities have to be suppressed. The procedure can be easily generalized to the case of more ambiguities to be suppressed (e.g., K < N). It is worth noting that the presented beamforming algorithm, following [14] and differently from [7,12,13,16,17], works in the space-time domain and not in the frequency domain. The algorithm is far from being efficient from the point of view of the computation burden since bistatic SAR focusing must be applied, separately, to each raw bistatic SAR image (i.e., N times). Herein, based on the unavailability of real FF-SAR data, the goal is not to develop or present a novel algorithm but to individuate an approach to test the developed PRF selection method on focused, either bistatic or multimonostatic, images.
A flowchart of the beamforming strategy is presented in Figure 5. Starting from the image grid, for each image pixel with coordinates ( ) where ( ) Tx R = − P P P and rE is the local Earth radius. Since herein the first ambiguity is only of interest, the location of the first forward P , A back and backward P , A forw ambiguities can be computed with m equal, respectively, to −1 and 1 in Equation (9). More details about the computation of the location of the azimuth ambiguities P A are available in Appendix A. The algorithm can be easily generalized to accommodate more accurate models of the Earth surface.
Using the computed P , A back and P , A forw , beamforming weights have to be defined. This can be done by setting the following conditions to the pattern of the sparse antenna formed by the location of the receivers at the slow time instant corresponding to the current pixel: the pattern must point toward P ; nulls must be placed in the directions of P , A back and P , A forw .
The conditions can be formulated as  a P P P (11) and the same equation can be used to compute the steering vector of the ambiguities, replacing P with the position of the ambiguity. Considering the phase of the current pixel as the reference phase to define the pattern, the steering vectors can be rewritten as It is evident that very accurate knowledge (i.e., at sub-wavelength scale) of the relative positions among the receivers is required to provide the relevant sparse antenna with the desired features. This accuracy, now, can be only achieved in postprocessing, so the beamforming algorithm is assumed to be applied on the ground processing bistatic data downloaded from each receiver.
The problem of Equation (9) Therefore, the weight can be computed using the concept of pseudoinverse of a nonsquare matrix where the superscript "H" indicates the Hermitian matrix. The weights of Equation (15) can be used to calculate a "partial" image , where Nc and Nr are respectively the number of columns and rows of the focused bistatic images.

Experimental Results
The developed approach to PRF selection in FF-SAR was tested on an interferometric data stack of Sentinel-1 images. The stack included 20 single-look complex images, acquired in VV polarization and in descending orbit phase ( Table 2). The considered acquisition mode was Interferometric Wide Swath (IW) [23], which acquired three subswaths of the imaged terrain at 5 m azimuth and 20 m range resolution. The pixel spacing of the processed data was 2.33 m in azimuth and 13.96 m in range. For simplicity, two bursts of the first subswath of each image was considered for further processing steps, characterized by about 34° incidence angle. Thus debursting was applied using SNAP toolbox [24] to merge the bursts into a single, continuous image. Images were coregistered in SNAP by backgeocoding and using precision orbit files. The master image was then selected as the one minimizing the dispersion of the perpendicular baseline (see Table 2). Figure 6 shows the master image over the tested area. The image covers the port of Naples, Italy, and the surrounding coastal area. A portion of the airport can be clearly identified. As expected, a man-made infrastructure close to the shoreline generated an evident azimuth ambiguity, which is imaged as a bright feature over the darker sea background. The beamforming algorithm discussed in Section 3 can be applied to the considered monostatic repeat-pass data stack. Namely, the sole difference with respect to FF-SAR operation is that the bistatic range must be replaced by the two-way monostatic path. As noted in Section 1, from the point of view of the algorithm, there is no further difference between FF-SAR processing and the monostatic repeat-pass one. Obviously, important differences arise from the point of view of the achieved results, which are discussed below.
Indeed, the available data stack is a good testbench for the developed PRF selection method. Specifically, the selected data stack is characterized by a single PRF value (which cannot be changed) (i.e., the same PRF is used in all the 20 considered satellite passes). Nonetheless, different combinations of the available images lead to different values of the J-index depending on the distribution of the corresponding samples. To study the behavior of the J-index, the locations of samples of each image corresponding to the first line time are analyzed and the κ j integer is computed for each platform. Conceptually, assuming a desired increase of PRF value by a factor of K < N, the J-index is computed for each possible combination without repetitions of all the N images in sets of K elements. As explained in Section 2, the goal is to select the best combination of K satellite passes for each value of K (i.e., the K images leading to the most uniform distribution of samples). Hence, the trend of J as a function of K is analyzed: Equation (8) 18) where N is the number of available images. Finally, the parameter K is chosen as the one minimizing JK as shown in Figure 7. Based on the analysis of J-index, and for each value of K, images corresponding to uniform samples are readily selected. The maximum degree of uniformity for the collected samples was obtained for K = 7. This means that, with reference to the PRF of the data stack, 7 out of 20 available images led to the most uniform distribution of samples. The beamforming algorithm of Section 3 was then applied to this case. It is clear that, because of decorrelation among the images collected in a repeat-pass fashion, the overall image quality of the computed final image was very low. So, to make the results less affected by time decorrelation, the analysis of the final image was limited to the high SNR, and stable azimuth ambiguity of Figure 6 and the relevant ground target. To appreciate the effect of the beamforming process, Figure 8 shows the signal intensity (in dB) along the iso-range line of the ground target generating the ambiguity. The strong reduction of the ambiguous signal after beamforming is evident. In order to quantify such an effect, the azimuth ambiguity-to-signal ratio (AASR) was evaluated. The AASR can be computed as ( ) 10 10 log where azimuth signal power is intended as the power in the azimuth main lobe and the power of the i-th azimuth ambiguity is the power in the lobe centered around the ambiguity. The summation in the numerator of Equation (19) should be evaluated from −∞ to +∞, however, only the first terms in the series (i.e., i close to ±1) significantly contribute to the value of the numerator.
The results show an AASR of -20.3 dB after the beamforming application with respect to an AASR of -4.5 dB of the master image.

Discussion
The experimental results presented in Section 4 showed that beamforming can be successfully applied to the subset of the available images characterized by the best J-index. To further point out the role played by the J-index, other possible subsets are herein analyzed. Figure 9a  x v PRF + . In the figure, the following notation is used: slave#| j κ to indicate, in addition to the slave number, also the relevant pulse number. It is evident that the distribution of samples for this configuration when all the images are used, referred to as Configuration A, is far from being uniform. If one indicates as Configuration B the one leading to the best value of the J-index, including just 7 out of 19 slaves, the corresponding distribution ( Figure 9b) is much more uniform. Finally, one can consider another configuration of 7 images, indicated as Configuration C (Figure 9c) (e.g., showing very limited baselines among some passes). Figure 9 indicates J-index values for all the selected configurations.
To further compare the selected configurations, the corresponding array patterns resulting from the beamforming are shown in Figure 10. Better null-steering capabilities are obtained when the Jindex is minimized. The array pattern in this case shows that the applied null reaches a much lower value and the application of the null affects a larger image portion around the ambiguity (see Configuration B in Figure 10). The computed values of AASR further confirm the validity of the proposed approach.
Concluding this section, it is important to remark that the analysis of the J-index, and more in general, the proper selection of the PRF, is of general validity, whereas only a simplified approach to FF-SAR beamforming/reconstruction has been tested because of the limitations of the available multimonostatic interferometric data stack. The considered beamforming algorithm only tends to suppress the ambiguities without accounting for the noise or other imaging parameters. Image quality for scenes with uniform backscattering or distributed scatterers is thus poor. Future activities are planned to test more general reconstruction algorithms such as those proposed in [13,16,17].

Conclusions
The paper analyzed the problem of PRF selection in a formation-flying SAR. A solution to compute a suitable PRF value was presented for any assigned formation of irregularly distributed The approximate solution can be refined by iteration. In the present work, convergence is assumed to be reached when . This is typically achieved in just 3-4 iterations even if P0 is tens of kilometers away from the true pixel location.
The gradient in Equation (A3) can be analytically derived x y z P P P P f P P P P P P P P P P P P P P P P P f P P P P P P P P P P P P P P P P P f P P P P P P P P