<figure><img src="image_1.png" /></figure>

Article

# A Comprehensive Comparison of Far-Field and Near-Field Imaging Radiometry in Synthetic Aperture Interferometry

Eric Anterrieu <fnref n="1" />*, Louise Yu <fnref n="2" /> and Nicolas Jeannin <fnref n="3" />

**Abstract:** Synthetic aperture interferometry (SAI) is a signal processing technique that mixes the signals collected by pairs of elementary antennas to obtain high-resolution images with the aid of a computer. This note aims at studying the effects of the distance between the synthetic aperture interferometer and an observed scene with respect to the size of the antenna array onto the imaging capabilities of the instrument. Far-field conditions and near-field ones are compared from an algebraic perspective with the aid of simulations conducted at microwave frequencies with the Microwave Imaging Radiometer by Aperture Synthesis (MIRAS) onboard the Soil Moisture and Ocean Salinity (SMOS) mission. Although in both cases the signals kept by pairs of elementary antennas are cross-correlated to obtain complex visibilities, there are several differences that deserve attention at the modeling level, as well as at the imaging one. These particularities are clearly identified, and they are all taken into account in this study: near-field imaging is investigated with spherical waves, without neglecting any terms, whereas far-field imaging approximation is considered with plane waves according to the Van–Citter Zernike theorem. From an algebraic point of view, although the corresponding modeling matrices are both rank-deficient, we explain why the singular value distributions of these matrices are different. It is also shown how the angular synthesized point-spread function of the antenna array, whose shape varies with the distance to the instrument, can be helpful for estimating the boundary between the Fresnel region and the Fraunhofer one. Finally, whatever the region concerned by the aperture synthesis operation, it is shown that the imaging capabilities and the performances in the near-field and far-field regions are almost the same, provided the appropriate modeling matrix is taken into account.

**Keywords:** imaging radiometry; antenna array; aperture synthesis; near-field; far-field

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Anterrieu, E.; Yu, L.; Jeannin, N. A Comprehensive Comparison of Far-Field and Near-Field Imaging Radiometry in Synthetic Aperture Interferometry.
*Remote Sens.* **2024**, *16*, 3584. https://doi.org/10.3390/rs16193584

Academic Editor: Massimiliano Pieraccini

Received: 5 August 2024

Revised: 17 September 2024

Accepted: 20 September 2024

Published: 26 September 2024

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2024 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

Regardless of the immediate environment of a receiving or radiating aerial, the space surrounding an antenna is usually subdivided into two regions: the near-field and the far-field regions [1], also known as the Fresnel and the Fraunhofer regions, respectively. While there is no clear boundary between these regions, a common and acceptable criterion is that $2D^2/\lambda$ represents a safe far-zone distance, where D is the maximum dimension of an antenna aperture and $\lambda$ is the operating wavelength [2]. At such a Fraunhofer distance, the maximum path-length difference between the contribution from the edge of the aperture and that from the center corresponds to the fraction $\lambda/16$ or to a phase $\pi/8$ radians, as very well illustrated by Figure 1.13 in [3]. However, as suggested by Equation (9) in [4], this criterion might not be sufficiently strict for aperture synthesis imaging applications. This is what is observed with the study reported in this note, which aims at comparing the well-known far-field conditions to the near-field ones in detail and to estimate if this criterion, which has been historically established for elementary antennas outside any imaging context, has to be revisited, or not, for antenna arrays when they are used for synthesis imaging with interferometric data.

<footnote n="1">¹ Centre d'Études Spatiales de la BIOsphère (CESBIO), 13 Avenue Colonel Roche, 31400 Toulouse, France</footnote>

<footnote n="2">² Centre National d'Études Spatiales (CNES), 18 Avenue Edouard Belin, 31400 Toulouse, France</footnote>

<footnote n="3">³ Airbus Defence and Space (ADS), 31 Rue des Cosmonautes, 31400 Toulouse, France</footnote>

<footnote>* Correspondence: eric.anterrieu@cnrs.fr</footnote>
Such a comparison between near-field and far-field imaging has already been addressed in synthetic aperture imaging radiometry. However, among the differences that can be found between these two situations, one or the other, or several particularities, are very often neglected. This is not the case in this study, where near-field imaging is considered with spherical waves and without neglecting any term in the modeling equation, whatever the price to pay for taking them all into account in the numerical processing conducted with the aid of a computer. On the other hand, far-field imaging approximation is considered with plane waves according to the Van Citter-Zernike theorem [5,6]. The question of the boundary between these two regions arises in synthetic aperture interferometry when complex visibilities of a scene observed at a given distance from the antenna array are inverted in order to retrieve the brightness temperature distribution of the synthesized field of view. Indeed, as soon as these complex visibilities have been obtained in the near-field region, they have to be inverted with the appropriate resolving matrix, namely the inverse of the modeling matrix of the instrument in near-field imaging conditions, not that in the far-field approximation; otherwise, errors might be expected in the retrieved map.

Regarding synthesis imaging operations with an array of elementary antennas, it is often suggested in the literature to convert near-field visibilities to far-field ones and then to invert the latter with a standard far-field resolving matrix, thus making the calculation and the inversion of the near-field imaging matrix unnecessary. However, the conditions of validity and reliability of this kind of pre-processing are not always mentioned. Although it obviously has many advantages in terms of computational effort, it is well suited only for unresolved point sources, but it is definitively not for extended ones. As argued in this note, these approaches are therefore useless in many cases. This provides clear support for the approach adopted in this study, which makes no approximation in the modeling and no compromise in the processing as far as near-field conditions are concerned. It aims at inverting near-field complex visibilities with the appropriate resolving matrix obtained from the near-field modeling matrix without neglecting any term in the corresponding equation and whatever the computational effort to calculate it and to invert it.

This note is devoted to the comparison between near-field and far-field imaging radiometry in synthetic aperture interferometry from a modeling point of view, from an algebraic point of view, and from an imaging point of view with the aid of numerical simulations. Section 2 gives an overview of the main differences between far-field and near-field synthesis imaging with an antenna array, as illustrated by a pair of elementary antennas with geometrical considerations. The modeling of a synthetic aperture interferometric radiometer is discussed in Section 3, where the integral relation between the brightness temperature distribution $T$ of an observed scene and the complex visibility $V$ for a pair of elementary antennas $A_p$ and $A_q$ is established for near-field imaging conditions when considering spherical waves as well as for far-field ones within the frame of the approximation with plane waves, as it is assumed in the Van Citter-Zernike theorem. At this level of the modeling equations, three differences between these two regions are clearly identified, and their potential impacts on the synthesis imaging are reported. The regularized reconstruction procedure that will be implemented for data processing is summarized in Section 4. To support the comparison, numerical simulations conducted at microwave frequencies with the Microwave Imaging Radiometer by Aperture Synthesis (MIRAS) onboard the Soil Moisture and Ocean Salinity (SMOS) mission are presented in Section 5. (In order to establish some consistency between this note and a previous one addressing a similar comparison under the same algebraic angle but between two observational paradigms in the far-field region only [7], the same array of elementary antennas have been used, although it is well known that MIRAS does not operate in the near-field region.) Although the two imaging conditions use the same interferometric paradigm and share the same goal, there are several differences that deserve attention while performing the comparison. From an algebraic point of view, this is the case of the corresponding modeling matrices, which are both rank-deficient with the same effective rank but with different distributions of their singular values. The impacts of the three differences identified in Section 3 are illustrated
and quantified. Complex visibilities of point sources are simulated and inverted, both in the near-field region and in the far-field one. It is shown that the angular resolution [8] does not depend on the region concerned by synthesis imaging. On the contrary, as soon as near-field visibilities are inverted with the far-field resolving matrix, it turns out that the shape of the angular synthesized point-spread function varies with the distance between the point-source and the phase center of the antenna array. A similar analysis is conducted with extended sources to quantify the so-called reconstruction floor error [9]. Finally, the radiometric sensitivity [10] is addressed, and simulations do not reveal any significant difference in the propagation of random Gaussian noise from the complex visibilities to the retrieved brightness temperature maps. For the sake of transparency and pedagogy, this note is completed with an annex, where the complex cross-correlation between two quasi-monochromatic electromagnetic fields is established for spherical waves. Although this study does not claim to define a boundary between the near-field region and the far-field one for synthetic aperture imaging radiometers, it is shown that for the antenna array of SMOS, the common criterion $2D^2/\lambda$ is under-estimating the extent of the near-field, or conversely over-estimating that of the far-field. Indeed, the safe far zone would probably be around $10D^2/\lambda$ for this antenna array when it is used for synthesis imaging with interferometric data.

## 2. Far-Field vs. Near-Field: A Geometrical Comparison

Since the very first successful experiment of interferometry in radio astronomy conducted after World War II by Ryle and Vonberg [11] with a two-element radio interferometer for studying the Sun at 175 MHz, synthetic aperture interferometry (SAI) techniques have evolved well, but the observing principle has not changed [12,13].

Shown in Figure 1 is the geometry of a one-dimensional SAI radiometer with two elementary antennas $A_p$ and $A_q$ observing a source under far-field conditions, as well as under near-field ones. The two antennas point in the same (Nadir) direction and are separated by the baseline vector $b_{pq}$. In the far-field region, owing to the distance from the source, the incident wavefronts are assimilable to plane waves coming from the same direction $\theta$ for both elementary antennas. As a consequence, the differential path $r_q - r_p$ is equal to $b_{pq} \sin \theta$, where $b_{pq}$ is the magnitude of the baseline vector $b_{pq}$. On the contrary, in the near-field region, incident wavefronts are spherical waves coming from different directions, $\theta_p$ for $A_p$ and $\theta_q$ for $A_q$, which are no longer equal to $\theta$. The differential path $r_q - r_p$ is here equal to $b_{pq} \sin \frac{1}{2}(\theta_p + \theta_q) / \cos \frac{1}{2}(\theta_p - \theta_q)$ which tends to $b_{pq} \sin \theta$ as soon as $\theta_p$ and $\theta_q$ both tend towards $\theta$ whenever far-field conditions are approached, as expected.

As a consequence of this very short preliminary comparison, SAI imaging radiometry might be more complicated under near-field conditions than under far-field ones because, for every baseline $b_{pq}$, it requires a precise calculation of the vectors $r_p(\theta, \phi)$ and $r_q(\theta, \phi)$ between the elementary antennas and any radiation coming from a direction $(\theta, \phi)$. It is also necessary to have an accurate estimation of the local incoming directions $(\theta_p, \phi_p)$ and $(\theta_q, \phi_q)$ to evaluate the appropriate attenuation by the elementary voltage patterns of $A_p$ and $A_q$.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Geometry of a one-dimensional SAI radiometer with two elementary antennas $A_p$ and $A_q$ and with incident wavefronts (in red) from a source: far-field conditions with plane waves (**left**) and near-field ones with spherical waves (**right**). In both cases, $A_p$ and $A_q$ are pointing in the same direction; here, the Nadir is illustrated by the main beam of the elementary power patterns (in blue), $r_p$ and $r_q$ are location vectors with respect to the source, and the two antennas are separated by the baseline vector $b_{pq}$.</figcaption></figure>

### 3. Instrument Modeling

According to Equation (20) page 573 in [12] and Equation (15.6) page 770 in [13], when neglecting decorrelation effects that can be reduced with spectral sub-banding [14] if narrow-band conditions are not satisfied [15], the complex visibility $V$ measured with a pair of elementary antennas $A_p$ and $A_q$ separated by the baseline vector $b_{pq}$, like those in Figure 1, is given by the integral relation established in the Appendix A at the end of this note in spherical coordinates $\theta$ and $\phi$:

$$V(\mathbf{b}_{pq}) = \iint_{source} \frac{j k_o (r_p(\theta, \phi) - r_q(\theta, \phi))}{\sqrt{\Omega_p} \sqrt{\Omega_q}} T(\theta, \phi) \frac{e}{r_p(\theta, \phi) r_q(\theta, \phi)} dS, \quad (1)$$

where $\mathcal{F}_p$ and $\mathcal{F}_q$ are the voltage patterns of the two elementary antennas with equivalent solid angles $\Omega_p$ and $\Omega_q$, $T$ is the brightness temperature distribution of an observed scene, $k_o \equiv 2\pi/\lambda_o$ is the central angular wavenumber, and $\lambda_o$ is the central wavelength of observation. The angular coordinates $(\theta_p, \phi_p)$ and $(\theta_q, \phi_q)$ are local spherical angles (the colatitude and the azimuth, respectively) from which the elementary brightness temperature contribution $T(\theta, \phi)dS$ is seen by the antennas $A_p$ and $A_q$, as illustrated in Figure 1. They are completed by the two distances $r_p$ and $r_q$ between the two antennas and the source element $dS$ responsible for that elementary contribution to the complex visibility. On introducing the solid angle $d\Omega \equiv dS/r^2$, we therefore have the general expression for $V$:

$$V(\mathbf{b}_{pq}) = \iint_{source} \frac{j k_o (r_p(\theta, \phi) - r_q(\theta, \phi))}{\sqrt{\Omega_p} \sqrt{\Omega_q}} T(\theta, \phi) \frac{r^2(\theta, \phi)}{r_p(\theta, \phi) r_q(\theta, \phi)} e^{j k_o (r_p(\theta, \phi) - r_q(\theta, \phi))} d\Omega, \quad (2)$$
where $r(\theta, \phi)$ is the distance from the origin of the reference frame of the antenna array (generally taken at the phase center of the array) to the source element $dS \equiv r^2 d\Omega$ responsible for the elementary contribution $T(\theta, \phi)dS$, as illustrated by Figure 1 in [4], Figure 2 in [16], Figure 3 in [17] and Equation (2) in [16] and Equation (2) in [18].

On introducing the angular position variable $\xi$ whose components $\xi_1 = \sin \theta \cos \phi$ and $\xi_2 = \sin \theta \sin \phi$ are direction cosines in the reference frame of the antenna array, the solid angle $d\Omega \equiv dS/r^2 \equiv \sin \theta d\theta d\phi$ turns out to be equal to $d\xi/\sqrt{1 - \|\xi\|^2}$ when accounting for the Jacobian of the transformation from spherical coordinates $(\theta, \phi)$ to direction cosines $(\xi_1, \xi_2)$, so Equation (2) now reads

$$V(\mathbf{b}_{pq}) = \iint_{\|\xi\| \le 1} \frac{\mathcal{F}_p(\xi_p) \mathcal{F}_q^*(\xi_q)}{\sqrt{\Omega_p} \sqrt{\Omega_q}} T(\xi) \frac{r^2(\xi)}{r_p(\xi) r_q(\xi)} e^{-2j\pi \frac{r_q(\xi) - r_p(\xi)}{\lambda_o}} \frac{d\xi}{\sqrt{1 - \|\xi\|^2}} \quad (3)$$

It is often read in the literature that expression (3) is valid only in near-field imaging conditions, but it is exactly the opposite: this expression is the most general one. Indeed, so far, no assumption has been made on the region of the electromagnetic field to be concerned so that it is valid in the near-field region as well as in the far-field one. However, from a numerical point of view, it is admitted that Expression (3) is not the easiest one to implement. This is why it is often used in the literature with approximations or assumptions to make it simpler.

Referring back to Figure 1, let us now assume that the source under observation is far enough from the antenna array so that $r_p \equiv r_q \equiv r$, $\theta_p \equiv \theta_q \equiv \theta$ and $\phi_p \equiv \phi_q \equiv \phi$ (or equivalently $\xi_p \equiv \xi_q \equiv \xi$). The differential path $r_q - r_p$ is therefore nothing but the scalar product $\mathbf{b}_{pq} \cdot \xi$ and as a consequence of these assumptions, expression (3) reduces to

$$\tilde{V}(\mathbf{b}_{pq}) = \iint_{\|\xi\| \le 1} \frac{\mathcal{F}_p(\xi) \mathcal{F}_q^*(\xi)}{\sqrt{\Omega_p} \sqrt{\Omega_q}} T(\xi) e^{-2j\pi \frac{\mathbf{b}_{pq} \cdot \xi}{\lambda_o}} \frac{d\xi}{\sqrt{1 - \|\xi\|^2}} \quad (4)$$

Expression (4) is the far-field expression of the complex visibility, also known in the literature as the Van Cittert–Zernike theorem [5,6]. The reasons for deriving such an approximation of the general expression (3) under some assumptions listed in page 774 of [13] are due to the observing conditions experienced in radio astronomy, which might be different from those encountered in Earth observation. This is particularly true with ground or airborne instruments that may sometimes need to be used in the Fresnel area, whereas spaceborne ones operate naturally in the Fraunhofer region. Moreover, contrary to (3), Expression (4) is very simple to implement as it is easier to know the baseline $\mathbf{b}_{pq}$ with some accuracy than the distances $r_p$ and $r_q$ for every viewing direction $\xi$.

When closely inspecting Expressions (3) and (4), it turns out that near-field conditions have three impacts on the modeling equation with respect to far-field approximation:

* The very first dissemblance between (3) and (4) is found in the phase term of the exponential where the difference $r_q(\xi) - r_p(\xi)$ is approximated by the scalar product $\mathbf{b}_{pq} \cdot \xi$ in the far-field approximation;

* The second difference is the significant presence of the fraction $r^2(\xi)/r_p(\xi)r_q(\xi)$ in (3), whereas it is absent in (4) since it tends toward 1 in the far-field approximation;

* The third and last dissemblance is found in the antenna voltage patterns $\mathcal{F}_p$ and $\mathcal{F}_q$ where local directions $\xi_p$ and $\xi_q$ of the local spherical angles $(\theta_p, \phi_p)$ and $(\theta_q, \phi_q)$ are assumed to be equal to $\xi$ in the far-field approximation.

As a consequence of the last item of this list, it is important to note that with an array of identical elementary antennas with the same voltage pattern $\mathcal{F}$, if synthesis imaging under far-field conditions can take into account this property to simplify the imaging process, it is no longer possible under near-field conditions. Indeed, for every
viewing direction $(\theta, \phi)$, the attenuation of the brightness temperature $T(\theta, \phi)$ by the voltage pattern $\mathcal{F}$ varies from one elementary antenna to another because the local incoming directions $(\theta_p, \phi_p)$ and $(\theta_q, \phi_q)$ are no longer equal. In the words of Zhang et al. [19] “that means for different baselines the source distribution is different, thus the Fourier transform cannot be used for the inversion.” From a modeling point of view, it is as if each antenna had a different diagram from the other one, according to the identities $\mathcal{F}_p(\theta, \phi) \equiv \mathcal{F}(\theta_p, \phi_p)$ and $\mathcal{F}_q(\theta, \phi) \equiv \mathcal{F}(\theta_q, \phi_q)$. It would therefore not be surprising to see the effect of this disparity of the elementary voltage patterns [20] on the properties of the modeling matrix of the instrument, as well as on the qualities of the inversions.

The previous effect cannot be seen in [21] nor in [22] as a “simplify visibility function ignoring the antenna voltage pattern” is used. Taking into account the antenna voltage pattern at the modeling level is not a guarantee of being able to capture this effect. This is the case in the recent studies reported in [23,24], where Equation (4) in the former and Equation (1) in the latter are both suffering from the same assumption that is not justified nor documented by the authors: the local incoming directions $\xi_p$ and $\xi_q$ corresponding to the local spherical angles $(\theta_p, \phi_p)$ and $(\theta_q, \phi_q)$ are assumed to be identical (and equal to $\xi$). This is in contradiction with expression (2), which is often referred to as the near-field complex visibility in the literature, like in Equation (1) of [16], which is very frequently cited in this imaging domain. Moreover, it is clear from Figure 1 that it cannot be the case when one claims to address near-field imaging. What is true with local spherical angles must remain true with direction cosines, as is the case in the present study with expressions (2) and (3). Finally, unlike the assumption made in [17,23] where the product of the distances between the two antennas and the source element dS is approximated by $r^2$ and therefore the corresponding fraction is neglected in the expression of the complex visibilities, the present note aims at using the complete expression (3) without neglecting any of the three terms listed previously and without introducing any reductive hypothesis, as is often the case when expending the difference $r_q(\xi) - r_p(\xi)$ in Taylor series [25] like in [19,24,26].

The usefulness of approximations or reductive assumptions in the near-field modeling is questionable, especially since there is no reason not to consider all the terms found in the general expression (3). This is also true for antenna arrays with an aperture so wide [27] that some baselines are in the near-field region while the others are in the far-field one, according to the criteria provided by the Fraunhofer distance [2]. As a consequence, this study focuses only on near-field and far-field synthesis imaging as illustrated by the modeling Equations (3) and (4) and by the corresponding linear operators $\mathbf{G}$ and $\tilde{\mathbf{G}}$ introduced in the next section.

## 4. Regularized Inversion

According to a recent revisit of (4), accounting for the mutual coupling between elementary antennas, the brightness temperature distribution $T(\xi)$ has to be substituted with the difference $T(\xi) - T_r$, where $T_r$ is the physical temperature of the receivers [28]. In any case, this constant term does not affect the zero-spacing baselines, and it can be removed by measuring a flat target [29]. As a consequence, without any loss of generality, after the discretization of the integral found in (4) over an appropriate sampling grid in the direction of the cosine domain [30], the relationship between the complex visibilities and the brightness temperature distribution of the scene under observation can be written in algebraic form:

$$ \tilde{V} = \tilde{\mathbf{G}} T \qquad (5) $$

where $\tilde{\mathbf{G}}$ is the linear modeling matrix of the imaging radiometer using SAI in the far-field approximation. Likewise, the general relationship (3) can be written in a very similar algebraic form:

$$ V = \mathbf{G} T \qquad (6) $$

where $\mathbf{G}$ is the general linear modeling matrix of the imaging radiometer using SAI, contrary to what is implied in the title of [31] as the difference $r_q(\xi) - r_p(\xi)$ is there
expended in Taylor series, thus resulting in an approximated linear modeling matrix, not a general one.

It is now well established that these modeling matrices are rank-deficient [32] as a consequence of a number of unknowns (the values of the brightness temperature distribution over the sampling grid) being larger than the number of data (the complex visibilities). In both paradigms, the effective rank of the modeling matrix is equal to the number of Fourier frequencies in the experimental frequency coverage $\mathcal{H}$ of the antenna array. Indeed, such imaging radiometers are band-limiting devices as the extent of the baselines $b_{pq}$ is limited to $\mathcal{H}$. As a consequence, the corresponding inverse problems have to be regularized in order to provide unique solutions to the linear systems (5) and (6). For obvious reasons of comparison, in this study, the same regularization is implemented in either case. Otherwise, the differences that might be observed could be attributed to the regularization itself and not to the modeling matrix (and, finally, not to the concerned electromagnetic region). Among the many regularization methods that can be found in the literature, the minimum-norm one is widely used in SAI [33]:

$$ \tilde{T}_r = \min_T \|T\|^2 \text{ s.t. } \tilde{G}T = \tilde{V}. \qquad (7) $$

The extension to the general case does not raise any difficulty:

$$ T_r = \min_T \|T\|^2 \text{ s.t. } GT = V. \qquad (8) $$

Numerical implementations are available in many algorithmic libraries and with many programming languages. In both cases, the reconstructed map is obtained through the computation of the pseudo-inverse [34] of the corresponding modeling matrix:

$$ \tilde{T}_r = \tilde{G}^+ \tilde{V}, \qquad (9) $$

$$ T_r = G^+ V. \qquad (10) $$

In either case, the pseudo-inverse is computed with the aid of a truncated singular value decomposition [35], where only the largest significant singular values are kept during the inversion.

Finally, in order to filter out the Gibbs effects [36] caused by the sharp cut-off in the Fourier domain due to the limited extent of the frequency coverage $\mathcal{H}$, these maps have to be damped by an appropriate apodization (or windowing) function $W$ [37]:

$$ \tilde{T}_r = W \star \tilde{T}_r, \qquad (11) $$

$$ T_r = W \star T_r, \qquad (12) $$

where $\star$ is the convolution operator. These maps have to be compared to

$$ T = W \star T, \qquad (13) $$

which is the “ideal” temperature map to be reconstructed and apodized with the same window $W$. Indeed, they cannot be compared to $T$ because they are not at the same angular resolution.

## 5. Numerical Simulations and Comparison

The numerical simulations discussed in this section have been performed with the antenna array of the imaging radiometer MIRAS [38] onboard the SMOS [39] mission. The reader will surely notice that this work and the conclusions drawn from the results presented in this section have evidently no impact on the operational processing of SMOS data as MIRAS is always operating in the far-field region and never in the near-field one. This antenna array has been chosen for feeding these numerical simulations just because it is
a concrete and tangible one with available voltage patterns measured for every elementary antenna [40], as well as for establishing some consistency with a previous note [7], as already mentioned in the Section 1. It should also be noted that the study reported in this note is missing real interferometric data from both electromagnetic regions with the same instrument to make the comparison complete. Such a possibility is currently not available, but it will be in the near future, as the French space agency (CNES) and the CESBIO are together developing an L-band antenna array equipped with 32 elementary antennas and with a longest baseline *D* of about 1 m intended to be used in a wide range of distances *h* (typically from 5 m to 40 m) suspended from a crane over an experienced and controlled field, as was performed few years ago with MIRAS breadboard [41]. Although the main goal of this forthcoming antenna array will not be to compare far-field and near-field imaging conditions but rather to confront synthetic aperture interferometry and digital beam forming [7], the flexibility of the instrument will allow such a comparison.

The SMOS mission of the European Space Agency (ESA) is devoted to Soil Moisture and Ocean Salinity remote sensing from space [39]. It was launched in November 2009, and since then, it has been operational and has provided accurate L-band brightness temperature maps [42] with a spatial resolution between ~25 km and ~60 km, depending on the position in the field of view. These maps are used for retrieving surface soil moisture (SM) [43], as well as sea surface salinity (SSS) [44]. The very principle of the imaging radiometer MIRAS is that of aperture synthesis, a technique initially developed for radio astronomy where the signals collected by each pair of interferometric arrays of elementary antennas are cross-correlated [13]. It is therefore not surprising that the antenna array of MIRAS was inspired by the Very Large Array (VLA) [45] for the geometry of the array and by the Cosmic Background Imager (CBI) [46] for the regular spacing between the antennas. As shown in Figure 2, it is a Y-shaped array equipped with 69 equi-spaced elementary antennas. As a consequence of the longest baseline *D* = 6.75 m and of the central wavelength of the instrument λ = 21.21 cm, the Fraunhofer distance 2D²/λ is here about 430 m [2]. Unless it is explicitly precised, in the remainder of this note, this antenna array will be placed in the near-field region at a distance *h* = 100 m from a source. However, unlike SMOS, which is pointing at the Earth with a constant tilt angle of 32.5° in the orbital plane [38,39], here, this offset is reduced to zero so that the number of pixels of the sampling grid in the synthesized field domain belonging to the Earth surface is constant over a wide range of values for *h*. This is essential for conducting studies where this distance is a key parameter; otherwise, the differences that might be observed could be attributed to this varying number of pixels and not to the modeling equations themselves. Referring back to the geometry of Figure 1 and to illustrate it concretely, shown in Figure 3 are the distances *r*<sub>p</sub>(*ξ*) from every antenna *A*<sub>p</sub> of this array to that source in the Nadir direction *ξ* = (0,0), as well as the corresponding local angles θ<sub>p</sub>(*ξ*). For this particular direction pointed at by every antenna, the distance to the source can be increased here by 8 cm, that is to say, larger than λ/3, which is not negligible. Likewise, whereas in far-field conditions, the colatitude angle corresponding to that direction is by definition equal to 0° for every antenna, at a distance *h* = 100 m, it varies from one antenna to the other and can be as large as 2.3°. As a consequence, these conditions are clearly not those of far-field imaging.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Location and numbering of the 69 elementary antennas of MIRAS in the reference frame of the instrument. The spacing between the elements along the three arms is equal to 0.875λ = 18.56 cm with λ = 21.21 cm (the central wavelength of the instrument). The longest baselines *D* between the antennas at the ends of the three arms are equal to 21√3 0.875λ = 6.75 m.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Distances $r_p(\xi)$ between a source located in the near-field region at a distance $h = 100$ m in the Nadir direction $\xi = (0,0)$ of the SMOS array and each of the 69 elementary antennas $A_p$ (**left**), together with the corresponding local colatitude angles $\theta_p(\xi)$ (**right**). The antenna index $p$ refers to the numbering of the elementary antennas of MIRAS shown in Figure 2.</figcaption></figure>

5.1. *Singular Values*

Referring back to (3) and (4), in both cases, the number of complex visibilities is
equal to the number of baselines **b**<sub>pq</sub>. For SMOS, it is equal to 69 × (69 − 1) = 4692. In
addition to these interferometric measurements, 3 out of the 69 receivers are also operat-
ing as reference radiometers to provide precise measurements of the average brightness
temperature of the observed scene [47]. The total number of measurements is therefore
equal to $4692 + 3 = 4695$. However, when accounting for the redundant baselines along the three arms of the Y-shaped array, the number of Fourier frequencies in the corresponding star-shaped frequency coverage $\mathcal{H}$ is only equal to $2790 + 1 = 2791$, including the zero frequency.

The distributions of the 4695 singular values of the modeling matrices $G$ and $\tilde{G}$ are shown in Figure 4. According to the floating point relative accuracy of MATLAB ($\epsilon \approx 2 \cdot 10^{-16}$), the 1904 smallest singular values of $G$ could be considered equal to zero [48]. The number of remaining strictly positive singular values is therefore equal to 4695 – 1904 = 2791, which is also the number of Fourier frequencies in the frequency coverage $\mathcal{H}$. This property remains true in the distribution of the singular values of the modeling matrix $G$; however, the gap between the two groups of singular values is much smaller: here less than one order of magnitude instead of 13 in the far-field conditions.

<figure><img src="image_7.png" /></figure>

**Figure 4.** Distributions of the 4695 singular values of the modeling matrices $G$ (red) and $\tilde{G}$ (green) of the SMOS array in X polarization when it is placed at a distance $h = 100$ m from a source. Whatever the electromagnetic region concerned by the modeling Equations (5) and (6), two groups of singular values separated by a well-determined gap are observed. In every case, the first group is composed of the 2791 largest singular values, and the second one is composed of the 1904 remaining singular values ($\epsilon \approx 2 \cdot 10^{-16}$ is the floating point relative accuracy of MATLAB). However, this gap between the two groups is much larger in the far-field case (green) than in the near-field one (red). The same behavior is observed in Y polarization.

As a consequence, when computing the pseudo-inverse with the aid of the truncated singular value decomposition, the impact of the singular values discarded prior to inversion might not be the same and a larger floor error might be observed in the near-field reconstructions compared to that in the far-field ones. Finally, the propagation of random noise is governed by the conditioning of the modeling matrix, which is the ratio between the largest and the smallest singular values kept in the inversion. It turns out that here, this ratio is about 11.2 for $\tilde{G}$ and 12.6 for $G$. Although this ratio is an upper bound for the actual amplification factor of random noise, one can expect very similar behaviors of the two matrices.

## 5.2. Far-Field vs. Near-Field

Referring back to Section 3, near-field conditions have three impacts on the modeling equation when comparing expressions (3) and (4). Shown in Figures 5 and 6 are examples of the spatial distribution of the error made when the difference $r_q(\xi) - r_p(\xi)$ is approximated by the dot product $b_{pq} \cdot \xi$. If the shape of this error is quite simple for pairs of elementary antennas located on the same arm, this is no longer the case when the two antennas are located on two different arms. With regards to the magnitude of this error, it is a not negligible fraction of the wavelength as high as 7.5 cm compared to the 3 mm accuracy required for the baselines coordinates at the central wavelength of SMOS $\lambda = 21.21$ cm. As
a consequence, such an approximation should have on strong impact on the brightness temperature retrievals.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Spatial distribution of the difference $\mathbf{b}_{pq} \cdot \mathbf{\zeta} - (r_q(\mathbf{\zeta}) - r_p(\mathbf{\zeta}))$ for three pairs of elementary antennas: $p = 3$ and $q = 23$ (**left**), $p = 26$ and $q = 46$ (**middle**), and $p = 49$ and $q = 69$ (**right**). Here, the array is placed at a distance $h = 100$ m from a source. The difference could be as high as 7.5 cm, that is to say, about $\lambda/3$ with $\lambda = 21.21$ cm (the central wavelength of SMOS antenna array).</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Spatial distribution of the difference $\mathbf{b}_{pq} \cdot \mathbf{\zeta} - (r_q(\mathbf{\zeta}) - r_p(\mathbf{\zeta}))$ for the three pairs of elementary antennas corresponding to the longest baselines: $p = 23$ and $q = 46$ (**left**), $p = 23$ and $q = 69$ (**middle**), and $p = 46$ and $q = 69$ (**right**). Here, the array is placed at a distance $h = 100$ m from a source. The difference could be as high as 2.5 cm, that is to say, about $\lambda/9$ with $\lambda = 21.21$ cm (the central wavelength of SMOS antenna array).</figcaption></figure>

Likewise, shown in Figure 7 are examples of the spatial distribution of the multi-plicative real-valued fraction $r^2(\mathbf{\zeta})/r_p(\mathbf{\zeta})r_q(\mathbf{\zeta})$ of (3). In the far-field approximation (4), this fraction is assumed to be equal to 1. In the near-field case shown here, the deviation from 1 could be as high as 4% for the shortest baselines at the extremities of the arms, but it decreases as soon as the antennas are closer to the phase center of the array. However, this is a not negligible fraction as the brightness temperature seen by the antenna array for these baselines from the corresponding directions is directly affected by this deviation. Although this term does not affect the phase of the complex visibilities, the impact on the magnitude is strong enough to affirm that it should not be neglected in near-field synthesis imaging.
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Spatial distribution of the fraction $r^2(\xi)/r_p(\xi)r_q(\xi)$ for the three pairs of elementary antennas corresponding to the shortest baseline at the extremities of the three arms: $p = 22$ and $q = 23$ (**left**), $p = 45$ and $q = 46$ (**middle**), and $p = 68$ and $q = 69$ (**right**). Here, the array is placed at a distance $h = 100$ m from a source. The deviation from 1 could be as high as 4%.</figcaption></figure>

Finally, shown in Figure 8 are examples of the spatial distribution of the magnitude and the phase of the difference $\mathcal{F}_p(\xi)\mathcal{F}_q^*(\xi) - \mathcal{F}_p(\xi_p)\mathcal{F}_q^*(\xi_q)$. The error introduced when ignoring local viewing directions $\xi_p$ and $\xi_q$, assuming $\xi = \xi_p = \xi_q$ in that part of the modeling equation involving the voltage pattern of the elementary antennascould be as high as 0.4 dB and 4°. Here again, this is large enough to affirm that local viewing directions at the level of the antenna's voltage pattern should be taken into account in near-field synthesis imaging.

<figure><img src="image_11.png" /><figcaption>**Figure 8.** Spatial distribution of the magnitude (**bottom**) and the phase (**top**) of the difference $\mathcal{F}_p(\xi)\mathcal{F}_q^*(\xi) - \mathcal{F}_p(\xi_p)\mathcal{F}_q^*(\xi_q)$ for the three pairs of elementary antennas corresponding to the shortest baseline at the extremities of the three arms: $p = 22$ and $q = 23$ (**left**), $p = 45$ and $q = 46$ (**middle**), and $p = 68$ and $q = 69$ (**right**). Here, the array is placed at a distance $h = 100$ m from a source. The difference could be as high as 0.4 dB and 4°.</figcaption></figure>
Keeping in mind that the objective of this study is to determine when the far-field approximation can be used and when near-field synthesis imaging has to be performed, shown in Figures 9–11 are the variations of the largest value of the previous spatial distributions with the observing distance *h* from the phase center of the antenna array. With regards to the difference **b**<sub>pq</sub>·ξ = (*r*<sub>q</sub>(ξ) − *r*<sub>p</sub>(ξ)), it turns out that for the Franhauser distance 2D<sup>2</sup>/λ ≈ 430 m, the largest value is here about λ/16, as expected, according to the common criterion for separating the Fresnel region from the Fraunhofer one. With respect to the fraction r<sup>2</sup>(ξ)/r<sub>p</sub>(ξ)r<sub>q</sub>(ξ), the largest deviation from the value 1 for the same Franhauser distance is slightly less than 1%. Finally, regarding the impact of the local viewing directions ξ<sub>p</sub> and ξ<sub>q</sub> onto the antenna voltage patterns F<sub>p</sub> and F<sub>q</sub>, the largest error introduced at the amplitude level is less than 0.1 dB for the same Franhauser distance.

<figure><img src="image_12.png" /><figcaption>**Figure 9.** Comparison between the difference $r_q(\xi) - r_p(\xi)$ and the dot product $\mathbf{b}_{pq} \cdot \xi$ (the far-field approximation) as a function of the distance between the source and the phase center of the antenna array. For a distance equal to the Fraunhofer distance $2D^2/\lambda$ with $D = 6.75$ m and $\lambda = 21.21$ cm (the longest baseline and the central wavelength of SMOS antenna array), the largest discrepancy between the two quantities is about $\lambda/16$.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 10.** Comparison between the ratio $r^2(\xi)/r_q(\xi)r_p(\xi)$ and 1 (the far-field approximation) as a function of the distance between the source and the antenna array. For a distance equal to the Fraunhofer distance $2D^2/\lambda$ with $D = 6.75$ m and $\lambda = 21.21$ cm (the longest baseline and the central wavelength of SMOS antenna array), the largest deviation from 1 is less than 1%.</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>**Figure 11.** Comparison between the cross-power patterns $\mathcal{F}_p(\xi_p)\mathcal{F}_q^*(\xi_q)$ and $\mathcal{F}_p(\xi)\mathcal{F}_q^*(\xi)$ (the far-field approximation) as a function of the distance between the source and the antenna array. For a distance equal to the Fraunhofer distance $2D^2/\lambda$ with $D = 6.75$ m and $\lambda = 21.21$ cm (the longest baseline and the central wavelength of SMOS antenna array), the largest discrepancy between the two quantities is less than 0.1 dB.</figcaption></figure>

As a consequence of this preliminary analysis, none of the three terms identified in Section 3 can be seen as a clear indicator of the boundary between the near-field region and the far-field one. Although the approximation of $r_q(\xi) - r_p(\xi)$ by $b_{pq} \cdot \xi$ leads to an error very close to that associated with the Fraunhofer distance, according to Equation (9) in [4], the threshold of $\lambda/16$ is questionable for imaging purposes with an antenna array: so why would a phase error of $\pi/8$ at complex visibilities level be acceptable? Referring back to the principle of synthesis imaging, the appropriate threshold has to be found at the level of the imaging process in the comparison made between the reconstructed maps $\tilde{T}_r$ and $T_r$ with the ideal one $\mathcal{T}$. This is precisely what is carried out in the next subsections, first with a point source, then with an extended one.

### 5.3. Point Source

Complex visibilities $V$ and $\tilde{V}$ have been simulated according to (3) and (4) for the same brightness temperature distribution $T(\xi)$ of a point-source, which is equal to 0 everywhere except in the boresight direction of the antenna array. They are shown with polar plots on Figure 12. Whatever the pair of elementary antennas $A_p$ and $A_q$, in such a special case of brightness temperature distribution, the far-field visibilities $\tilde{V}(b_{pq})$ are real-valued, as expected according to (4) as soon as the only contribution to the integral is for $\xi = (0, 0)$. This is no longer the case with the near-field visibilities $V(b_{pq})$ simulated when the antenna array is placed at a distance $h = 100$ m from the point source. Although the instrument is still the same, the phase of these complex visibilities varies in the range $\pm 130^\circ$, as a consequence of the effects of the geometry illustrated in Figure 1. Furthermore, simulations conducted with (3) for different values of $h$ have shown that this range of phases is larger as the distance between the antenna array and the point-source gets closer. On the contrary, it tends toward zero as this distance increases and reaches the far-field region, as expected and as shown in Figure 13. For the Fraunhofer distance of $2D^2/\lambda = 430$ m, it turns out that this range of phases is about $\pm 30^\circ$. Referring back to the far-field approximation, such a phase translates into a baseline error as large as $\lambda/12 \simeq 18$ mm. This is 6 times larger than the 3 mm baseline uncertainty of SMOS, which corresponds to a phase error of the order of $5^\circ$ that is reached here only from a distance $h$ greater than 2 Km, that is to say, about $10D^2/\lambda$. As a consequence, there is no doubt that inverting the near-field complex visibilities $V$ with the pseudo inverse $\tilde{G}^+$ of the far-field modeling operator will lead to large errors and artifacts not only for $h = 100$ m, which is evidently in the near-field region, but also for any value of $h \lesssim 10D^2/\lambda$, including of course for the Fraunhofer distance $2D^2/\lambda$.
<figure><img src="image_15.png" /><figcaption>**Figure 12.** Complex visibilities of an unresolved point-source when the array is in the far-field region (**left**) and when it is placed at a distance $h = 100$ m from the source (**right**). Attention should focus on the phase of the far-field visibilities $\tilde{V}$, which is equal to zero for every baseline $b_{pq}$ whereas that of the near-field ones $V$ varies in the range $\pm 130^\circ$ with an average value of $0^\circ$, depending on the pair of elementary antennas $A_p$ and $A_q$.</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>**Figure 13.** The range of the phase of the complex visibilities $V$ of an unresolved point-source as a function of the distance $h$ between the source and the antenna array. Attention should focus on the decrease towards $0^\circ$ as $h$ increases. However, for the Fraunhofer distance $2D^2/\lambda = 430$ m, this range is still about $\pm 30^\circ$. It remains below $\pm 5^\circ$ as soon as $h$ exceeds 2 Km $\simeq 10D^2/\lambda$.</figcaption></figure>

The previous complex visibilities $\tilde{V}$ and $V$ have been inverted according to (9) and (10). Shown in Figure 14 is the reconstructed point-spread function of SMOS when the antenna array is placed at a distance $h = 100$ m from the source in the near-field region, and the visibilities $V$ shown in red in Figure 12 are inverted with the pseudo-inverse of the near-field modeling operator $G$ or with that of the far-field one $\tilde{G}$. Both have to be compared to the point-spread function obtained when the array is in the far-field region and the corresponding visibilities $\tilde{V}$ shown in green in Figure 12 are inverted with $\tilde{G}^+$. In every case, a rectangular windowing function $W$ has been used so that the reconstructed point-spread function is not apodized and the angular resolution, as estimated by the full-width at half-maximum of the main lobe, is the finest one that could be obtained from this antenna array.

As soon as the complex visibilities are inverted with the appropriate operator, the one corresponding to the concerned electromagnetic region, the retrieved point-spread functions exhibit the expected imaging qualities. More precisely, $G^+V$ in the near-field case and $\tilde{G}^+\tilde{V}$ in the far-field one are perfectly identical, and they evidence the predicted values for many figures of merit (full-width at half-maximum of the main lobe, level of the side lobes, etc.). On the contrary, when inverting near-field visibilities with a far-field
operator, it turns out that $\tilde{G}^+V$ is completely different from the expected point-spread function. Clearly, the far-field imaging operator is responsible for large aberrations in the near-field point-spread function. This is exactly what happens when far-field imaging conditions are assumed to be met when they are not. As the energy transported by the electromagnetic field is the same because it does not depend on the distance but only on the brightness temperature of an observed scene, the norm of the retrieved maps is also the same. However, it is clear that the spatial distribution of this energy is distorted so that the shape of the point-spread function has significantly changed. It is no longer that of a main lobe surrounded by six-fold patterns with an alternating succession of positive and negative side lobes, but it is now similar to the shape of the antenna array itself. Finally, as illustrated by Figure 15, this shape also depends on the location of the point source in the field of view synthesized by the instrument.

<figure><img src="image_17.png" /><figcaption>**Figure 14.** Reconstructed point-spread function when the array is in the far-field region and the visibilities $\tilde{V}$ are inverted with the far-field modeling operator $\tilde{G}^+$ (**left**), when the array is placed at a distance $h = 100$ m from the source at the boresight $\zeta = (0,0)$ and the visibilities $V$ are inverted with the far-field modeling operator $\tilde{G}^+$ (**middle**) or with the near-field one $G^+$ (**right**). The color scale is the same for all images, and contour levels are drawn for -0.10, 0.10, 0.25 and 0.50.</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>**Figure 15.** Reconstructed point-spread function when the array is in the near-field region at a distance $h = 100$ m from the boresight point and the near-field visibilities $V$ are inverted with the far-field modeling operator $\tilde{G}^+$. The reconstruction varies with the position $\zeta$ of the source in the field of view: from left to right, $\zeta = (-0.5, 0), \zeta = (0, 0)$ and $\zeta = (0.5, 0)$. The color scale is the same for all images, contour levels are drawn for -0.10, 0.10, 0.25 and 0.50.</figcaption></figure>
Shown in Figure 16 are reconstructed point-spread functions of SMOS when the antenna array is placed at a distance $h = 2D^2/\lambda$, $3D^2/\lambda$, $4D^2/\lambda$, $5D^2/\lambda$, $10D^2/\lambda$ and $20D^2/\lambda$ from a source and when the complex visibilities $V$ are inverted with the pseudo-inverse of the far-field modeling operator $\tilde{G}$. As predicted earlier in this subsection when looking at the range of the phase of the near-field visibilities $V$ reported in Figure 13, large errors and artifacts are observed here at the level of the reconstructed point-spread functions for any value of $h \lesssim 10D^2/\lambda$. Although the Mean Squared Error (MSE) has many attractive features, it is not well suited for quantifying the sameness, or on the contrary the difference, between images [49]. The Structural SIMilarity (SSIM) index [50] is a well-known method for measuring the similarity, or the dissimilarity, between two images. This index can be viewed as a quality measure of one of the images being compared, provided the other image is considered of perfect quality. Here, the reference image is the reconstructed point-spread function shown in Figure 14 when the antenna array is in the far-field region and the visibilities $\tilde{V}$ are inverted with the far-field modeling operator $\tilde{G}^+$. It is compared to every reconstructed point-spread function, like those shown in Figure 16, obtained when the array is placed at a distance of $h$ from the source and when the corresponding visibilities $V$ are nevertheless inverted with the far-field modeling operator $\tilde{G}^+$. Shown in Figure 17 are the variations in the SSIM index obtained for this comparison. For the Fraunhofer distance $2D^2/\lambda = 430$ m, this index is only about 0.82, but it is larger than 0.99 for any value of $h \gtrsim 10D^2/\lambda$.

<figure><img src="image_19.png" /><figcaption>**Figure 16.** Reconstructed point-spread function when the array is placed at a distance $h = 2D^2/\lambda$, $3D^2/\lambda$, $4D^2/\lambda$, $5D^2/\lambda$, $10D^2/\lambda$, and $20D^2/\lambda$ from a source at boresight $\xi = (0,0)$, with $D = 6.75$ m and $\lambda = 21.21$ cm (the longest baseline and the central wavelength of SMOS antenna array). In every case, the near-field visibilities $V$ are inverted with the far-field modeling operator $\tilde{G}^+$. The color scale is the same for all images, and contour levels are drawn for $-0.10, 0.10, 0.25$ and $0.50$.</figcaption></figure>
<figure><img src="image_20.png" /><figcaption>**Figure 17.** Structural SIMilarity (SSIM) index between the far-field point-spread function and the reconstructed point-spread functions obtained when the array is placed at a distance of *h* from the source and when the corresponding visibilities are inverted with the far-field modeling operator. Attention should focus on the increase towards 1 as *h* increases. For the Fraunhofer distance $2D^2/\lambda = 430$ m, this index is only about 0.82. It is larger than 0.95 as soon as *h* exceeds 1 Km $\approx 5D^2/\lambda$ and larger than 0.99 for any $h \gtrsim 2$ Km $\approx 10D^2/\lambda$.</figcaption></figure>

Finally, it is frequently read in the literature (Equation (5) in [16], Equation (9) in [19] or Equation (2) in [51], for example) that near-field complex visibilities $V$ can be corrected from their distortion to obtain far-field visibilities $\tilde{V}$ by simply subtracting the near-field phase and adding the far-field one for every baseline $\mathbf{b}_{pq}$ according to

$$V(\mathbf{b}_{pq}) e^{+2j\pi \frac{r_q(\xi) - r_p(\xi)}{\lambda_o}} e^{-2j\pi \frac{\mathbf{b}_{pq} \cdot \xi}{\lambda_o}} \approx \tilde{V}(\mathbf{b}_{pq}). \quad (14)$$

However, referring back to Equations (3) and (4), only a point source can be concerned by such a correction because these modeling equations are integral relations so that the complex visibility for any antenna pair $A_p$ and $A_q$ does not depend on the variable $\xi$. It should be read between the lines in [16,19] that $\xi$ is here nothing but the viewing direction from a single point source; more precisely, the brightness temperature distribution $T(\xi)$ in (3) and (4) is assumed here to be null everywhere except for that particular value of $\xi$. To be complete with this correction, it is also assumed that the second and the third impact out of the three listed in Section 3 are neglected. This is exactly what should be understood in the words of [16] "in the case that only path delay phase errors are taken into account... if we take into account a single punctual source." Such an approach is therefore useless since observed scenes are not always reduced to a single point source. Moreover, applying this phase correction to an extended source or to a collection of point sources is simply not possible, unless the same phase correction is applied for every viewing direction in the field of view, which would assume that it does not change from one direction to another. This is exactly what is suggested in [16], where "near-field visibilities from an extended source can be partially corrected just by focussing the near-field visibilities to the pixel at boresight." This is not always true as illustrated here by Figure 5, where, for three pairs of elementary antennas, the difference between $\mathbf{b}_{pq} \cdot \xi$ and $r_q(\xi) - r_p(\xi)$ can be as large as $\pm \lambda/3$ between the boresight direction and other viewing directions in the field of view. Thus, for the corresponding baselines, if the phase correction (14) is correctly introduced at boresight, an error as large as $\pm 2\pi/3$ is not corrected in the complex visibility for the contributions from other directions. The same remark holds for Figure 6, where, for three other pairs of elementary antennas, this phase error is as large as $\pm 2\pi/9$, which is still not negligible. As a consequence, owing to the role played by the phase of the complex visibilities in imaging synthesis operations, large errors should be expected in the regions, or for the
point sources, far from the boresight direction. This is exactly what is shown in Figure 18, where the brightness temperature distribution of three point sources located at $\xi = (0,0)$ and at $\tilde{\xi} = (\pm0.5,0)$ is observed when the antenna array is in the near-field region at a distance $h = 100$ m from the first source at boresight (the two other sources are therefore at a distance $h/\cos 30^\circ \simeq 115$ m from the phase center). The near field visibilities $V$ are corrected according to (14) with $\xi = (0,0)$, and the corresponding “almost” far-field visibilities $\tilde{V}$ thus obtained are inverted with the far-field modeling operator $\tilde{G}^+$. As expected, the reconstructed point-spread function at boresight is identical to the expected one already shown in Figure 14. On the contrary, as soon as the point source moves away from the very center of the field of view, the reconstructed point-spread function is still suffering from strong aberrations, although the situation is obviously slightly better than that shown in Figure 15 without the phase correction (14). Accounting for the previous considerations and referring back to Section 4, it is recommended to invert near-field complex visibilities $V$ with the appropriate imaging operator $G^+$ rather than trying to convert them in far-field ones $\tilde{V}$ to be inverted with $\tilde{G}^+$, contrary to what is suggested in [51] where it is claimed that “after near-field correction, far-field image retrieval tools can be used without any additional modification.”

<figure><img src="image_21.png" /><figcaption>**Figure 18.** Reconstructed point-spread function when the array is in the near-field region at a distance $h = 100$ m from the boresight point and the near-field visibilities $V$ are corrected from their distortion to obtain far-field visibilities, which are inverted with the far-field modeling operator $\tilde{G}^+$. The reconstruction varies with the position $\tilde{\xi}$ of the source in the field of view: from left to right, $\tilde{\xi} = (-0.5,0)$, $\tilde{\xi} = (0,0)$ and $\tilde{\xi} = (0.5,0)$. The color scale is the same for all images, contour levels are drawn for $-0.10, 0.10, 0.25$ and $0.50$.</figcaption></figure>

### 5.4. *Extended Source*

Near-field imaging is no exception to field aliasing; as a consequence, unlike what is often found in the literature, it has to be studied and compared to what is observed and very well known when far-field imaging conditions are met. Here again, complex visibilities $V$ and $\tilde{V}$ have been simulated according to (3) and (4) for the same brightness temperature distribution $T(\tilde{\xi})$ of an extended source with a constant temperature. Both have been inverted with the appropriate operator: $G^+$ for $V$ according to (10) and $\tilde{G}^+$ for $\tilde{V}$ according to (9).

Two cases have been considered, depending on the extent of the source. The first one corresponds to a situation that has been widely simulated with a circular [16] or a square [24] flat target that never extends beyond the field of view synthesized by the antenna array. The brightness temperature distribution $T(\tilde{\xi})$ of an extended source set to 100 K in the field of view synthesized by the antenna array of Figure 2 and to 0 K outside this hexagon is shown in Figure 19, together with the corresponding complex visibilities $V$ and $\tilde{V}$.
<figure><img src="image_22.png" /><figcaption>**Figure 19.** Brightness temperature distribution $T(\xi)$ of an extended source that is set to 100 K in the field of view synthesized by the antenna array of Figure 2 and to 0 K outside this hexagon. Shown in green are the corresponding complex visibilities $\tilde{V}$ when the antenna array is in the far-field region, while those in red are the complex visibilities $V$ when the phase center of the array is placed at a distance $h = 100$ m from the source in the boresight direction.</figcaption></figure>

Shown in Figure 20 are the brightness temperature maps $\tilde{T}_r$ and $T_r$ retrieved after inverting the complex visibilities $\tilde{V}$ and $V$ of Figure 19 according to (9) and (10). In every case, a Blackman windowing function $W$ has been used to dampen the reconstructed maps $\tilde{T}_r$ and $T_r$. As the brightness temperature $T$ of the observed scene is here set to 0 K outside the field of view synthesized by the antenna array, field aliasing is not an issue because it is not expected. Indeed, there is no reason to observe it in such conditions where the brightness temperature from the six neighbors to the synthesized field of view does not contribute to the complex visibilities. As a consequence, in the absence of any modeling error or noise, there is no significant difference observed between $\mathcal{T} = W \star T$ and the two retrieved maps $\tilde{T}_r$ and $T_r$, as expected. In both cases, the floor error as illustrated here by the Root Mean Squared Error (RMSE) [52] in the whole synthesized field of view is much less than 0.01 K, and it does not vary with the distance $h$ between the antenna array and the source in the boresight direction.

<figure><img src="image_23.png" /><figcaption>**Figure 20.** Brightness temperature maps $\tilde{T}_r$ (far-field, left) and $T_r$ (near-field, right) retrieved after inverting the complex visibilities $\tilde{V}$ and $V$ shown in Figure 19. In both cases, a Blackman windowing function $W$ has been used to dampen the reconstructed maps $\tilde{T}_r$ and $T_r$. To facilitate the comparison, the color scale is the same for these two images and for the brightness temperature $T$ shown in Figure 19: there is no significant difference between $\mathcal{T} = W \star T$ and the two retrieved maps shown here (in both cases, the RMSE is much less than 0.01 K).</figcaption></figure>
The second case considered here corresponds to a situation that has not been widely simulated nor reported so far: this is the case where the brightness temperature distribution of the observed scene extends beyond the field of view synthesized by the antenna array. In such a situation, field aliasing has to be observed, depending on the geometry of the antenna array and on the spacing between the elementary antennas. The brightness temperature distribution $T(\xi)$ of an extended source set to 100 K in the whole space in front of the antenna array of Figure 2 is shown in Figure 21, together with the corresponding complex visibilities $V$ and $\tilde{V}$.

<figure><img src="image_24.png" /><figcaption>**Figure 21.** Brightness temperature distribution $T(\xi)$ of an extended source that is set to 100 K in front of the antenna array of Figure 2. Shown in green are the corresponding complex visibilities $\tilde{V}$ when the antenna array is in the far-field region, while those in red are the complex visibilities $V$ when the phase center of the array is placed at a distance $h = 100$ m from the source in the boresight direction.</figcaption></figure>

Shown in Figure 22 are the brightness temperature maps $\tilde{T}_r$ and $T_r$ retrieved after inverting the complex visibilities $\tilde{V}$ and $V$ of Figure 21 according to (9) and (10). Here again, in every case, a Blackman windowing function $W$ has been used to damp the reconstructed maps $\tilde{T}_r$ and $T_r$. As the brightness temperature $T$ of the observed scene is not set to 0 K outside the field of view synthesized by the antenna array, field aliasing is expected from the six neighbors to that hexagon. This is exactly what can be observed in Figure 22, where only the central part of the hexagonal field of view is free from any aliasing artifact. However, even in the absence of any modeling error or noise, small differences can be observed between $T = W \star T$ and the two retrieved maps $\tilde{T}_r$ and $T_r$. More precisely, here, the floor error in that alias-free region of the synthesized field of view is about 0.14 K for the far-field case, whereas it is of the order of 0.53 K for the near-field one. In real situations, this floor error can be efficiently reduced down to a much lower level with the introduction of an artificial scene [53], as close as possible to the observed one, whose visibilities are removed prior to inversion and whose brightness temperature distribution is added back after inversion. However, referring back to Figure 4, such a difference at this raw level of reconstruction was already expected when comparing the singular values of $G$ and $\tilde{G}$ and especially those discarded prior to the pseudo-inversion, which are quite different. Finally, shown in Figure 23 are the variations of this floor error with the observing distance $h$ from the phase center of the antenna array. Although it decreases very rapidly to reach a plateau with a value very similar to that obtained in far-field conditions, it is still very difficult to find a value that makes a clear separation between near-field and far-field regions. With this figure of merit, one can again at most consider a transition zone from the Fresnel region to the Fraunhofer one between $2D^2/\lambda$ and $10D^2/\lambda$.
<figure><img src="image_25.png" /></figure>

**Figure 22.** Brightness temperature maps $\tilde{T}_r$ (far-field, **left**) and $T_r$ (near-field, **right**) retrieved after inverting the complex visibilities $\tilde{V}$ and $V$ shown in Figure 21. In both cases, a Blackman windowing function $W$ has been used to dampen the reconstructed maps $\tilde{T}_r$ and $T_r$. Shown in dashed brown are the horizons from the six neighbors to the hexagonal field of view in red, which are responsible for field-aliasing. To facilitate the comparison, the color scale is the same for these two images and for the brightness temperature $T$ shown in Figure 19. In the central part of the field of view free from field-aliasing, the distribution of the retrieved brightness temperature is 100.04 K ± 0.14 K (far-field, **left**) and 100.21 K ± 0.49 K (near-field, **right**); these numbers have to be compared to the expected 100 K.

<figure><img src="image_26.png" /></figure>

**Figure 23.** Root Mean Square Error (RMSE) in the alias-free field of view as a function of the distance $h$ between the source and the antenna array when observing an extended source with a constant temperature set to 100 K in the whole space in front of the instrument (red). Attention should focus on the decrease towards the floor error obtained in far-field conditions (0.14 K, green) as $h$ increases. For the Fraunhofer distance $2D^2/\lambda = 430$ m, this floor error is about 0.19 K, i.e., 35% larger than the far-field value. This excess remains below 5% as soon as $h$ exceeds $2 \text{ Km} \approx 10D^2/\lambda$.

5.5. *Radiometric Sensitivity*

The real and the imaginary parts of the complex visibilities $V$ and $\tilde{V}$ have been blurred by an additive Gaussian noise $\Delta V$ and $\Delta \tilde{V}$ with standard deviation $\sigma_{\Delta V} = \sigma_{\Delta \tilde{V}} = 0.1$ K. This is a typical value for SMOS, which should translate in the brightness temperature domain with a radiometric sensitivity on the order of 2.5 K [10]. This is what is observed in Figure 24 where the radiometric sensitivity in the boresight direction $\xi = (0,0)$ is equal to 2.55 K in far-field imaging conditions and to 2.57 K in near-field ones. Moreover, the radiometric sensitivity map does not depend on the distance between the antenna array
and the scene observed by the instrument. Referring back to Figure 4, this is an expected
behavior when comparing the singular values of G and $\tilde{G}$, and especially the role played
by those kept in the pseudo-inversion, which are almost identical.

<figure><img src="image_27.png" /><figcaption>**Figure 24.** An example of random Gaussian noise $\Delta V$, or $\Delta \tilde{V}$, added to the complex visibilities $V$ and $\tilde{V}$ here with $\langle \Delta V \rangle = \langle \Delta \tilde{V} \rangle = 0$ K and $\sigma_{\Delta V} = \sigma_{\Delta \tilde{V}} = 0.1$ K (**left**) and the corresponding spatial distributions of the standard deviations of 1000 error maps $\Delta \tilde{T}_r$ (**middle**) and $\Delta T_r$ (**right**). In both cases, a Blackman windowing function $W$ has been used to dampen the reconstructed maps $\Delta \tilde{T}_r$ and $\Delta T_r$. The radiometric sensitivity in the boresight direction $\zeta = (0,0)$ is equal to $\sigma_{\Delta \tilde{T}} = 2.55$ K and $\sigma_{\Delta T} = 2.57$ K.</figcaption></figure>

**6. Conclusions**

A comprehensive comparative study between near-field and far-field imaging radiom-
etry in synthetic aperture interferometry has been conducted. At the modeling level, the
integral relation between the brightness temperature distribution of an observed scene
and the corresponding complex visibilities has been established for near-field imaging
conditions when considering spherical waves as well as for far-field ones with plane waves.
Three differences between the two regions have been clearly identified and their potential
impacts on the synthesis imaging have been reported. The first dissemblance is at the level
of the difference between the distances from any pixel in the field of view to any pair of
elementary antennas that is estimated in the far-field approximation by the scalar product
between the corresponding baseline and the viewing direction. The second difference is
the significant presence of a multiplicative fraction that tends towards 1 as soon as the
far-field conditions are met. The third dissemblance is at the level of the antennas' voltage
pattern where local viewing directions vary from one antenna to another in near-field
imaging, whereas they are assumed to be equal in the far-field approximation. All these
differences have been taken into account throughout this study, and their impact has been
quantified. From an algebraic point of view, the properties of the corresponding linear
operators have been compared with an emphasis on their singular values. As expected,
the modeling matrices are both rank-deficient with the same effective rank. However,
they exhibit significantly different distributions of their singular values: those kept in the
inversion are very similar, whereas those discarded prior to the inversion are not. As a
consequence, the propagation of random noise does not vary from far-field to near-field
conditions, as it is governed by the former singular values. This is not the case for the
reconstruction floor error, which depends on the latter ones. At the imaging synthesis level,
numerical simulations have been conducted for the antenna array of SMOS in order to
compare near-field and far-field imaging performances when observing a point source
or when looking at an extended source. In the first situation, no significant difference
has been observed provided the visibilities are inverted with the appropriate resolving
matrix. Nevertheless, when inverting near-field visibilities with the far-field operator,
strong errors have been observed in the reconstructed point-spread functions, as expected. In the second case, when the extended source is limited to the field of view synthesized by the antenna array, no difference has been found between near-field and far-field imaging: the reconstruction floor error is close to zero in both imaging conditions. On the contrary, when the brightness temperature distribution of an observed scene extends beyond the synthesized field of view, field aliasing is observed and the reconstruction floor error is at a higher level for near-field imaging than it is for the far-field one. Finally, with regard to the propagation of random noise, no difference has been found between near-field and far-field inversions, as expected.

Although it is admitted that there is no clear boundary between the Fresnel and the Fraunhofer regions, the well-known distance $2D^2/\lambda$ is often used as a criterion for separating the near-field region from the far-field one. This criterion has been historically established for elementary antennas with diameter $D$ and operating at a wavelength $\lambda$ outside any imaging context. It was not therefore unreasonable to take advantage of this comparative study to check how this distance is consistent with a collection of elementary antennas working all together for synthesis imaging from interferometric data and when $D$ is now the longest baseline of the antenna array. This is why parametric studies with the distance between the phase center of the antenna array and an observed source in the boresight direction have been conducted with the idea of confronting different figures of merit to this historical criterion. In any case, no sharp separation has been observed between the near-field and the far-field regions: one can at most consider a transition zone between the Fresnel and the Fraunhofer regions. However, as already suggested in the literature, it turns out that $2D^2/\lambda$ is not sufficiently strict for aperture synthesis imaging applications since it is always found to be among the lowest values of this transition interval between the two regions. Referring back one last time to Figures 13, 17 and 23, a criterion around $10D^2/\lambda$ might be more appropriate to preserve the imaging performances. However, it would be very presumptuous to claim to change the definition of this criterion solely on the basis of a comparative study conducted with a single antenna array. The lesson learned from this parametric study should be to always invert the complex visibilities of a synthetic aperture imaging radiometer with the most appropriate resolving matrix and to consider the far-field approximation only after quantifying its impact on the imaging performances such as the angular resolution, the floor error and the radiometric sensitivity.

**Author Contributions:** Conceptualization, E.A.; methodology, E.A.; software, E.A.; validation, E.A., L.Y. and N.J.; formal analysis, E.A., L.Y. and N.J.; investigation, E.A., L.Y. and N.J.; resources, E.A. and L.Y.; data curation, E.A.; writing—original draft preparation, E.A.; writing—review and editing, E.A. and L.Y.; visualization, E.A.; supervision, E.A.; project administration, E.A. and L.Y. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by the french national research council (CNRS) within the frame of the FANTASIOR (Flexible ANTenna Array for Synthetic Imaging in Observational Radiometry) project led by the Centre d'Études Spatiales de la BIOsphère (CESBIO) and funded by the french space agency (CNES) under the R&T action R-S24/OT-0003-101-01.

**Data Availability Statement:** The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

**Acknowledgments:** The authors are very grateful to Yann Kerr (CESBIO) and to Cédric Tourain (CNES) for their constant support and encouragement during this work. Data processing support was provided by CESBIO and CNES. The authors are indebted to the associate editor and to the four anonymous reviewers for their helpful reviews which greatly improved the manuscript.

**Conflicts of Interest:** Author Dr. Nicolas Jeannin is employed by the company Airbus Defense & Space. All authors declare that this research was conducted in the absence of any business or financial relationship that could be interpreted as a potential conflict of interest.
## Appendix A. Mutual Coherence of an Incoherent Source

The expression (1) of the complex visibility for a pair of separated antennas observing the same brightness-temperature distribution is established in this annex with the approach and the reasoning found in radio astronomy [13] but with the symbols and the notations used in this note to facilitate understanding for remote sensing readers. Furthermore, the voltage pattern of the two elementary antennas is here taken into account, whereas it is not in chapter 15 of [13] nor in chapter 10 of [12], which are, however, the most cited references for applications of the Van Cittert-Zernike theorem in the radio domain and in that of optics, respectively.

The quasi-monochromatic electromagnetic fields at the elementary antennas $A_p$ and $A_q$ coming from a source element dS in the direction $(\theta, \phi)$ are given by:

$$E_p(\theta, \phi, t) = \mathcal{E}(\theta, \phi, t - \frac{r_p}{c}) \frac{e^{-2j\pi v_o \left(t - \frac{r_p}{c}\right)}}{r_p}, \qquad (A1)$$

and

$$E_q(\theta, \phi, t) = \mathcal{E}(\theta, \phi, t - \frac{r_q}{c}) \frac{e^{-2j\pi v_o \left(t - \frac{r_q}{c}\right)}}{r_q}, \qquad (A2)$$

where $\mathcal{E}(\theta, \phi, t)$ is the complex amplitude of the electromagnetic field at the level of the source dS found in the direction $(\theta, \phi)$: $|\mathcal{E}|$ is the strength and $\arg(\mathcal{E})$ is the phase of the radiation, $r_p \equiv \|\mathbf{r}_p\| = \|\mathbf{SA}_p\|$ and $r_q \equiv \|\mathbf{r}_q\| = \|\mathbf{SA}_q\|$ are the distances from dS to $A_p$ and $A_q$, respectively. The exponential term in Equations (A1) and (A2) represents the phase change when traveling from dS to $A_p$ and $A_q$, whereas $1/r_p$ and $1/r_q$ represents the attenuation of the amplitude during the same paths. In the words of R. Feynman in chapter 20 of [54], each expression (A1) and (A2) “... represents a general spherical wave travelling outward from the origin at the speed c.” There is a similar expression with $(t + r/c)$ instead of $(t - r/c)$, which “...also represents a spherical wave, but one which travels inward from large r toward the origin,” but it is here obviously without any interest.

The complex cross-correlation [55] with zero time offset of these two electromagnetic fields, when accounting for the attenuation by the voltage patterns $\mathcal{F}_p$ and $\mathcal{F}_q$ of the two elementary antennas, as illustrated in Figure 1, is

$$\begin{aligned} & \langle \mathcal{F}_p(\theta_p, \phi_p) E_p(\theta, \phi, t) \mathcal{F}_q^*(\theta_q, \phi_q) E_q^*(\theta, \phi, t) \rangle \\ &= \mathcal{F}_p(\theta_p, \phi_p) \mathcal{F}_q^*(\theta_q, \phi_q) \langle \mathcal{E}(\theta, \phi, t - \frac{r_p}{c}) \mathcal{E}^*(\theta, \phi, t - \frac{r_q}{c}) \rangle \frac{e^{\frac{2j\pi v_o \frac{r_p - r_q}{c}}{r_p r_q}}}{r_p r_q}, \\ &= \mathcal{F}_p(\theta_p, \phi_p) \mathcal{F}_q^*(\theta_q, \phi_q) \langle \mathcal{E}(\theta, \phi, t) \mathcal{E}^*(\theta, \phi, t - \frac{r_q - r_p}{c}) \rangle \frac{e^{\frac{2j\pi v_o \frac{r_p - r_q}{c}}{r_p r_q}}}{r_p r_q}, \end{aligned}$$

As soon as $|r_q - r_p|/c$ is small compared to the reciprocal of the receiver bandwidth, it can be neglected in the complex amplitude $\mathcal{E}^*$ so that the previous equation now reads

$$\begin{aligned} & \langle \mathcal{F}_p(\theta_p, \phi_p) E_p(\theta, \phi, t) \mathcal{F}_q^*(\theta_q, \phi_q) E_q^*(\theta, \phi, t) \rangle \\ &= \mathcal{F}_p(\theta_p, \phi_p) \mathcal{F}_q^*(\theta_q, \phi_q) \langle \mathcal{E}(\theta, \phi, t) \mathcal{E}^*(\theta, \phi, t) \rangle \frac{e^{\frac{2j\pi v_o \frac{r_p - r_q}{c}}{r_p r_q}}}{r_p r_q}, \end{aligned}$$
where the quantity $\langle \mathcal{E}(\theta, \phi, t) \mathcal{E}^*(\theta', \phi', t) \rangle = \langle |\mathcal{E}(\theta, \phi, t)|^2 \rangle$ is nothing but the time-averaged intensity of the source $\mathcal{I}(\theta, \phi)$:

$$ \begin{aligned} & \langle \mathcal{F}_p(\theta_p, \phi_p) E_p(\theta, \phi, t) \mathcal{F}_q^*(\theta_q, \phi_q) E_q^*(\theta, \phi, t) \rangle \\ & = \mathcal{F}_p(\theta_p, \phi_p) \mathcal{F}_q^*(\theta_q, \phi_q) \mathcal{I}(\theta, \phi) \frac{e^{2j\pi v_o \frac{r_p - r_q}{c}}}{r_p r_q}. \end{aligned} \quad (A3) $$

The source is assumed to be spatially incoherent, which means that cross-contributions like $\langle \mathcal{E}(\theta, \phi, t) \mathcal{E}^*(\theta', \phi', t) \rangle$, coming from different source elements $dS$ and $dS'$ are equal to zero. As a consequence, the complex visibility associated with the antenna pair $A_p$ and $A_q$ is obtained by integrating the previous expression (A3) over the source:

$$ V_{pq} = \iint_{\text{source}} \mathcal{F}_p(\theta_p, \phi_p) \mathcal{F}_q^*(\theta_q, \phi_q) \mathcal{I}(\theta, \phi) \frac{e^{2j\pi v_o \frac{r_p - r_q}{c}}}{r_p r_q} dS, \quad (A4) $$

which is nothing but Equation (1) since $2\pi v_o/c = 2\pi/\lambda_o = k_o$ and the intensity $\mathcal{I}$ of the source is related to its brightness temperature $T$ with Planck's law [56] according to $\mathcal{I} = 2v_o^2 kT/c^2$ in the Rayleigh-Jeans approximation at low frequencies [57,58].

To be complete and without any ambiguity, it is important to examine the assumptions and the limitations that are involved, or not, in this expression. The only point that could raise a question is the magnitude of $|r_q - r_p|/c$ compared to that of the reciprocal of the receiver's bandwidth $\Delta v$. For an antenna array operating in the 1400-1427 MHz protected band, $c/\Delta v$ is about 10 m, which means that this assumption is valid for a compact ground instrument. For a wider instrument, the receiver's bandwidth has to be decomposed into smaller ones [14] to push this limit.

## References

1. Balanis, C.A. Antenna Theory: Analysis and Design, 3rd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2005.

2. Selvan, K.T.; Janaswamy, R. Fraunhofer and Fresnel Distances: Unified Derivation for Aperture Antennas. IEEE Antennas Propag. Mag. 2017, 59, 12–15. [CrossRef]

3. Haupt, R.L. Antenna Arrays: A Computational Approach, 1st ed.; John Wiley & Sons: Hoboken, NJ, USA, 2010.

4. Laursen, B.; Skou, N. Synthetic Aperture Radiometry Evaluated by a Two-Channel Demonstration Model. IEEE Trans. Geosci. Remote Sens. 1998, 36, 822–832. [CrossRef]

5. van Cittert, P.H. Die Wahrscheinliche Schwingungsverteilung in Einer von Einer Lichtquelle Direkt Oder Mittels Einer Linse Beleuchteten Ebene. Physica 1934, 1, 201–210. [CrossRef]

6. Zernike, F. The Concept of Degree of Coherence and its Application to Optical Problems. Physica 1934, 5, 785–795. [CrossRef]

7. Anterrieu, É.; Lafuma, P.; Jeannin, N. An Algebraic Comparison of Synthetic Aperture Interferometry and Digital Beam Forming in Imaging Radiometry. Remote Sens. 2022, 14, 2285. [CrossRef]

8. Bará, J.; Camps, A.; Torres, F.; Corbella, I. Angular resolution of two-dimensional, hexagonally sampled interferometric radiometers. Radio Sci. 1998, 33, 1459–1473. [CrossRef]

9. Duran, I.; Lin, W.; Corbella, I.; Torres, F.; Duffo, N.; Martín-Neira, M. SMOS floor error impact and mitigation on ocean imaging. In Proceedings of the IEEE International Geoscience And Remote Sensing Symposium (IGARSS 2015), Milan, Italy, 26–31 July 2015.

10. Camps, A.; Corbella, I.; Bará, J.; Torres, F. Radiometric sensitivity computation in aperture synthesis interferometric radiometry. IEEE Trans. Geosci. Remote Sens. 1998, 36, 680–685. [CrossRef]

11. Ryle, M.; Vonberg, D.D. Solar Radiation at 175 Mc/s. Nature 1946, 158, 339–340. [CrossRef]

12. Born, M.; Wolf, E. Principles of Optics, 7th ed.; Cambridge University Press and Assessment: Cambridge, UK, 1999.

13. Thompson, A.R.; Moran, J.W.; Swenson, G.W. Interferometry and Synthesis in Radio Astronomy, 3rd ed.; Springer International Publishing: Cham, Switzerland, 2017.

14. Fischman, M.A.; England, A.W. A technique for Reducing Fringe Washing Effects in L-Band Aperture Synthesis Radiometry. In Proceedings of the IEEE International Geoscience And Remote Sensing Symposium (IGARSS 2000), Honolulu, HI, USA, 24–28 July 2000.

15. Zatman, M. How Narrow is Narrowband? IEE Proc.-Radar Sonar Navig. 1998, 145, 85–91. [CrossRef]
16. Tanner, A.B.; Lambrigsten, B.H.; Gaier, T.M.; Torres, F. Near Field Characterization of the GeoSTAR Demonstrator. In Proceedings of the IEEE International Geoscience And Remote Sensing Symposium (IGARSS 2006), Denver, CO, USA, 31 July–4 August 2006.

17. Peichl, M.; Suess, H.; Suess, M.; Kern, S. Microwave Imaging of the Brightness Temperature Distribution of Extended Areas in the Near and Far Field Using Two-Dimensional Aperture Synthesis with High Spatial Resolution. Radio Sci. **1998**, *33*, 781–801.
[CrossRef]

18. Liangbing, C.; Qingxia, L.; Zubiao, X.; Ke, C.; Yaoting, Z. Near-Field Imaging of One-Dimensional Aperture Synthesis Radiometers. In Proceedings of the IEEE International Symposium on Microwave, Antenna, Propagation and EMC Technologies for Wireless Communications (MAPE 2009), Beijing, China, 27–29 October 2009.

19. Zhang, C.; Wu, J.; Liu, H.; Yan, J. Imaging Algorithm for Synthetic Aperture Interferometric Radiometer in Near Field. Sci. China Technol. Sci. **2011**, *54*, 2224–2231. [CrossRef]

20. Díez-García, R.; Martín-Neira, M. Antenna spacing and pattern differences: Their impact in MIRAS reconstruction error. In Proceedings of the 14th Specialist Meeting on Microwave Radiometry and Remote Sensing of the Environment (MicroRad 2016), Espoo, Finland, 11–14 April 2016.

21. Fu, P.; Hu, F.; Hu, H.; Zheng, T. A Wavenumber Domain Imaging Algorithm for Synthetic Aperture Interferometric Radiometry in Near-Field. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS 2020), Waikoloa, HI, USA, 26 September–2 October 2020.

22. Fu, P.; Zhu, D.; Hu, F.; Xu, Y.; Xia, H. A Near-Field Imaging Algorithm Based on Angular Spectrum Theory for Synthetic Aperture Interferometric Radiometer. IEEE Trans. Microw. Theory Tech. **2022**, *70*, 3606–3616. [CrossRef]

23. Hu, H.; Zhu, D.; Hu, F. A Novel Imaging Method Using Fractional Fourier Transform for Near-Field Synthetic Aperture Radiometer Systems. IEEE Geosci. Remote Sens. Lett. **2022**, *19*, 1–5. [CrossRef]

24. Wu, Y.; Li, Y.; Song, G.; Dou, H.; Li, H. A Near-Field Imaging Method Based on the Near-Field Distance for an Aperture Synthesis Radiometer. Remote Sens. **2024**, *16*, 767. [CrossRef]

25. Feigenbaum, L. Brook Taylor and the method of increments. Arch. Hist. Exact Sci. **1985**, *34*, 1–140. [CrossRef]

26. Chen, J.; Li, Y.; Wang, J.; Li, Y.; Zhang, Y. Regularization Imaging Algorithm with Accurate G Matrix for Near-Field MMW Synthetic Aperture Imaging Radiometer. Prog. Electromagn. Res. B **2014**, *58*, 193–203. [CrossRef]

27. Schnoering, G.; Kawaguchi, T.; Höller, C. Mixed Near-and Far-Field 2D Imaging with Sparse Arrays. In Proceedings of the IEEE Radar Conference (RadarConf24), Denver, CO, USA, 6–10 May 2024.

28. Corbella, I.; Duffo, N.; Vall-llossera, M.; Camps, A.; Torres, A. The Visibility Function in Interferometric Aperture Synthesis Radiometry. IEEE Trans. Geosci. Remote Sens. **2004**, *42*, 1677–1682. [CrossRef]

29. Martín-Neira, M.; Suess, M.; Kainulainen, J.; Martin-Porqueras, F. The Flat Target Transformation. IEEE Trans. Geosci. Remote Sens. **2008**, *46*, 613–620. [CrossRef]

30. Anterrieu, É. A Resolving Matrix Approach for Synthetic Aperture Imaging Radiometers. IEEE Trans. Geosci. Remote Sens. **2004**, *42*, 1649–1656. [CrossRef]

31. Chen, J.; Zhu, X.; Zhang, S.; Chen, X. General G-matrix imaging method for near-field millimeter-wave SAIR with any arrays. In Proceedings of the IEEE International Wireless Symposium (IWS 2018), Chengdu, China, 6–10 May 2018.

32. Hansen, P.C. Rank-Deficient and Discrete Ill-Posed Problems, 1st ed.; Society for Industrial & Applied Mathematics: Philadelphia, PA, USA, 1998.

33. Goodberlet, M.A. Improved Image Reconstruction Techniques for Synthetic Aperture Radiometers. IEEE Trans. Geosci. Remote Sens. **2000**, *38*, 1362–1366. [CrossRef]

34. Golub G.H.; Van Loan, C.F. Matrix Computations, 4th ed.; Johns Hopkins University Press: Baltimore, MD, USA, 2013.

35. Hansen, P.C. The Truncated SVD as a Method for Regularization. BIT Numer. Math. **1987**, *27*, 534–553. [CrossRef]

36. Fay, T.H.; Kloppers, P.H. The Gibbs' Phenomenon. Int. J. Math. Educ. Sci. Technol. **2001**, *32*, 73–89. [CrossRef]

37. Anterrieu, É.; Waldteufel, P.; Lannes, A. Apodization Functions for 2D Hexagonally Sampled Synthetic Aperture Imaging Radiometers. IEEE Trans. Geosci. Remote Sens. **2002**, *40*, 2531–2542. [CrossRef]

38. McMullan, K.D.; Brown, M.A.; Martín-Neira, M.; Rits, W.; Ekholm, S.; Lemanczyk, J. SMOS: The Payload. IEEE Trans. Geosci. Remote Sens. **2008**, *46*, 594–605. [CrossRef]

39. Barré, H.; Duesmann, B.; Kerr, Y.H. SMOS: The Mission and the System. IEEE Trans. Geosci. Remote Sens. **2008**, *46*, 587–593. [CrossRef]

40. Pivnenko, S.; Nielsen, J.M.; Cappellin, C.; Lemanczyk, G.; Breinbjerg, O. High-Accuracy Calibration of the SMOS Radiometer Antenna Patterns at the DTU-ESA Spherical Near-Field Antenna Test Facility. In Proceedings of the IEEE International Geoscience And Remote Sensing Symposium (IGARSS 2007), Barcelona, Spain, 23–28 July 2007.

41. Bayle, F.; Wigneron, J.-P.; Kerr, Y.H.; Waldteufel, P.; Anterrieu, E.; Orlhac, J.-C.; Chanzy, A.; Marloie, O.; Bernardini, M.; Sobjaerg, S.; Calvet, J.-C.; Goutoule, J.-M.; Skou, N. Two-dimensional synthetic aperture images over a land surface scene. IEEE Trans. Geosci. Remote Sens. **2002**, *40*, 710–714. [CrossRef]

42. Talone, M.; Gourrion, J.; Sabia, R.; Gabarró, C.; Gonzalez, V.; Camps, A.; Corbella, I.; Monerris, A.; Font J. SMOS' brightness temperatures validation: First results after the commissioning phase. In Proceedings of the IEEE International Geoscience And Remote Sensing Symposium (IGARSS 2010), Honolulu, HI, USA, 25–30 July 2010.

43. Kerr, Y.H.; Waldteufel, P.; Richaume, P.; Wigneron, J.-P.; Ferrazzoli, P.; Mahmoodi, A.; Al Bitar, A.; Cabot, F.; Gruhier, C.; Enache Juglea, S.; et al. The SMOS Soil Moisture Retrieval Algorithm. IEEE Trans. Geosci. Remote Sens. **2012**, *50*, 1384–1403. [CrossRef]
44. Zine, S.; Boutin, J.; Font, J.; Reul, N.; Waldteufel, P.; Gabarró, C.; Tenerelli, J.; Petitcolin, F.; Vergely, J.-L.; Talone, M.; et al. Overview of the SMOS Sea Surface Salinity Prototype Processor. *IEEE Trans. Geosci. Remote Sens.* **2008**, *46*, 621–645.

45. Thompson, A.R.; Clark, B.G.; Wade, C.M.; Napier, P.J. The Very Large Array. *Astrophys. J. Suppl. Ser.* **1980**, *44*, 151–167. [CrossRef]

46. Padin, S; Shepherd, M.C.; Cartwright, J.K.; Keeney, R.G.; Mason, B.S.; Pearson, T.J.; Readhead, A.C.S.; Schaal, W.A.; Sievers, J.; Udomprasert, P.S.; et al. The Cosmic Background Imager. *Publ. Astron. Soc. Pac.* **2002**, *114*, 83–97.

47. Colliander, A.; Tauriainen, S.; Auer, T.I.; Kainulainen, J.; Uusitalo, J.; Toikka, M.; Hallikainen, M.T. MIRAS Reference Radiometer: A fully polarimetric noise injection radiometer. *IEEE Trans. Geosci. Remote Sens.* **2005**, *43*, 1135–1143. [CrossRef]

48. The MathWorks Inc.: Help Center. Available online: https://fr.mathworks.com/help/matlab/ref/eps.html (accessed on 1 August 2023).

49. Wang, Z.; Bovik, A.C. Mean squared error: Love it or leave it? A new look at signal fidelity measures. *IEEE Signal Process. Mag.* **2009**, *26*, 98–117. [CrossRef]

50. Wang, Z.; Bovik, A.C.; Sheikh, H.R.; Simoncelli, E.P. Image quality assessment: From error visibility to structural similarity. *IEEE Trans. Image Process.* **2004**, *13*, 600–612. [CrossRef]

51. Torres, F.; Abril, J.; Nova, E.; Broquetas, A.; Jofre, L. Far field condition in passive interferometry for security screening applications. In Proceedings of IEEE the European Conference on Antennas and Propagation (EuCAP 2010), Barcelona, Spain, 12–16 April 2010.

52. Huber P.J.; Ronchetti, E.M. *Robust Statistics*, 2nd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2009.

53. Khazâal, A.; Carfantan, H.; Anterrieu, É. On the reduction of the systematic error in imaging radiometry by aperture synthesis: A new approach for the SMOS space mission. *IEEE Geosci. Remote Sens. Lett.* **2009**, *6*, 47–51. [CrossRef]

54. Feynman, R.P.; Leighton, R.B.; Sands, M. *The Feynman Lectures on Physics, Vol. II: Mainly Electromagnetism and Matter*, New Millennium ed.; Basic Books: New-York, NY, USA, 2011.

55. Bracewell, R. *The Fourier Transform and Its Applications*, 3rd ed.; McGraw-Hill: New-York, NY, USA, 1999.

56. Planck, M. Ueber das Gesetz der Energieverteilung im Normalspectrum. *Ann. Phys.* **1901**, *309*, 553–563. [CrossRef]

57. Rayleigh, J. Remarks upon the Law of Complete Radiation. *Philos. Mag.* **1900**, *49*, 539–540. [CrossRef]

58. Jeans, J. On the Partition of Energy between Matter and Aether. *Philos. Mag.* **1905**, *10*, 91–98. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.