Spatial Ambiguity Function (SAF)

• Spatial Ambiguity Function is a metric that characterizes an array’s ability to distinguish spatial positions and directions by comparing shifted aperture fields.
• It captures the inner-product response between displaced replicas of the aperture field, revealing couplings between linear and angular resolutions, especially in near-field regimes.
• SAF analysis informs practical trade-offs in array design by addressing aperture geometry, aliasing constraints, and Cramér–Rao bounds for localization accuracy.

The Spatial Ambiguity Function (SAF) quantifies the capability of an array (aperture) to distinguish different spatial configurations of a monochromatic signal source, subject to shifts in spatial position and direction of arrival. It is the spatial analog of the classical Doppler–delay ambiguity function, underlying the resolution limits and coupling between linear displacement and spatial frequency (or angle), including effects of aperture geometry, radiating backgrounds, and array sampling. The SAF is foundational in the analysis and design of sensor arrays for localization, synthetic aperture radar, and large-scale multi-antenna systems, particularly as next-generation arrays operate in near-field regimes where simple far-field approximations are inadequate (Yuryev et al., 2014, Monnoyer et al., 9 Dec 2025).

1. Mathematical Definition and Fundamental Properties

The classical definition of the SAF begins with a two-dimensional aperture with complex field-weight distribution $h(\mathbf r)$, $\mathbf r \in \mathbb R^2$, and received spatial signal $$H(\mathbf r) = \int_{\mathbb R^2} A(\boldsymbol\eta) e^{-j\,\boldsymbol\eta \cdot \mathbf r} d^2\eta$$, with $A(\boldsymbol\eta)$ the spatial frequency (directional) spectrum and $\lambda$ the wavelength. To analyze the effect of spatial (linear) and directional (spatial-frequency) perturbations, the function

$$\Xi(\Delta\mathbf x,\, \Delta\boldsymbol\eta) = \int_R h(\mathbf r + \tfrac{1}{2} \Delta\mathbf x)\, h^*(\mathbf r - \tfrac{1}{2} \Delta\mathbf x)\, e^{-j\,\Delta\boldsymbol\eta \cdot \mathbf r}\, d^2r,$$

is introduced, normalized such that $\Xi(\mathbf 0, \mathbf 0) = 1$. In one spatial dimension: $$\Xi(\Delta x,\; \Delta k) = \int_{-L/2}^{+L/2} h(x + \tfrac{\Delta x}{2})\, h^*(x - \tfrac{\Delta x}{2})\, e^{-j \Delta k x}\, dx.$$ This function characterizes the inner-product response between two shifted replicas of the aperture field, modulated by a relative phase gradient. The SAF quantifies the array’s ability to discriminate between a source at $(x, k)$ and one at $(x+\Delta x, k+\Delta k)$ in position and spatial frequency (Yuryev et al., 2014).

2. Far-Field and Near-Field Forms

In the far-field (Fraunhofer) regime, the SAF simplifies by identifying $\mathbf r \in \mathbb R^2$0, with $\mathbf r \in \mathbb R^2$1 the angular shift: $\mathbf r \in \mathbb R^2$2 This form admits classical "sinc" behavior for uniform illumination: $\mathbf r \in \mathbb R^2$3 yielding main-lobe widths $\mathbf r \in \mathbb R^2$4, $\mathbf r \in \mathbb R^2$5.

In the near-field (Fresnel) regime, phase curvature is significant: $\mathbf r \in \mathbb R^2$6 The resulting SAF peak is no longer at $\mathbf r \in \mathbb R^2$7 but is displaced along $\mathbf r \in \mathbb R^2$8, reflecting the coupling between displacement and spatial frequency. The lobe widths scale as $\mathbf r \in \mathbb R^2$9, $$H(\mathbf r) = \int_{\mathbb R^2} A(\boldsymbol\eta) e^{-j\,\boldsymbol\eta \cdot \mathbf r} d^2\eta$$0 (Yuryev et al., 2014, Monnoyer et al., 9 Dec 2025).

3. Role of Extended and Reradiated Backgrounds

If the spatial scene reradiates or scatters with random, spatially-extended amplitude pattern $$H(\mathbf r) = \int_{\mathbb R^2} A(\boldsymbol\eta) e^{-j\,\boldsymbol\eta \cdot \mathbf r} d^2\eta$$1, the measured SAF manifests as a product of (i) the point-target SAF and (ii) the background-induced field’s autocorrelation,

$$H(\mathbf r) = \int_{\mathbb R^2} A(\boldsymbol\eta) e^{-j\,\boldsymbol\eta \cdot \mathbf r} d^2\eta$$2

where $$H(\mathbf r) = \int_{\mathbb R^2} A(\boldsymbol\eta) e^{-j\,\boldsymbol\eta \cdot \mathbf r} d^2\eta$$3 or, in the spatial domain, $$H(\mathbf r) = \int_{\mathbb R^2} A(\boldsymbol\eta) e^{-j\,\boldsymbol\eta \cdot \mathbf r} d^2\eta$$4 for average background-induced field $$H(\mathbf r) = \int_{\mathbb R^2} A(\boldsymbol\eta) e^{-j\,\boldsymbol\eta \cdot \mathbf r} d^2\eta$$5. This causes the overall SAF volume to shrink relative to the ideal value, and modifies the resolvable widths: the array’s ability to distinguish linear or angular displacements depends jointly on the aperture and on the spatial statistics of the background (Yuryev et al., 2014).

4. Resolution, Cramér–Rao Bounds, and Estimate Coupling

Resolution in $$H(\mathbf r) = \int_{\mathbb R^2} A(\boldsymbol\eta) e^{-j\,\boldsymbol\eta \cdot \mathbf r} d^2\eta$$6 and $$H(\mathbf r) = \int_{\mathbb R^2} A(\boldsymbol\eta) e^{-j\,\boldsymbol\eta \cdot \mathbf r} d^2\eta$$7 (or angle) is analytically derived from the second-order derivatives of $$H(\mathbf r) = \int_{\mathbb R^2} A(\boldsymbol\eta) e^{-j\,\boldsymbol\eta \cdot \mathbf r} d^2\eta$$8 at the origin: $$H(\mathbf r) = \int_{\mathbb R^2} A(\boldsymbol\eta) e^{-j\,\boldsymbol\eta \cdot \mathbf r} d^2\eta$$9 Given received SNR $A(\boldsymbol\eta)$0, the Cramér–Rao lower bounds on unbiased estimation error variances are

$A(\boldsymbol\eta)$1

There is also a cross-covariance

$A(\boldsymbol\eta)$2

Closed-form results for variance and correlation coefficients can be expressed in terms of the root-mean-square (rms) widths $A(\boldsymbol\eta)$3 (aperture field), $A(\boldsymbol\eta)$4 (aperture spectrum), and $A(\boldsymbol\eta)$5 (background spectrum width): $A(\boldsymbol\eta)$6 and coupling coefficient $A(\boldsymbol\eta)$7 (Yuryev et al., 2014).

5. Discrete Array Sampling, Aliasing, and Grating Lobes

When the aperture is discretely sampled, the SAF must account for the induced aliasing structure. For a grid $A(\boldsymbol\eta)$8 with element spacing $A(\boldsymbol\eta)$9, the discrete SAF reads: $\lambda$0 where $\lambda$1 is the matched-filter kernel between putative source locations $\lambda$2 and $\lambda$3.

Taking the spatial Fourier transform yields

$\lambda$4

Thus, discrete sampling folds the spatial spectrum $\lambda$5 every $\lambda$6 and generates grating lobes wherever the support of $\lambda$7 exceeds $\lambda$8 (Monnoyer et al., 9 Dec 2025).

The soft bandlimit

$\lambda$9

quantifies the maximum local spatial frequency. The aliasing-free condition is $$\Xi(\Delta\mathbf x,\, \Delta\boldsymbol\eta) = \int_R h(\mathbf r + \tfrac{1}{2} \Delta\mathbf x)\, h^*(\mathbf r - \tfrac{1}{2} \Delta\mathbf x)\, e^{-j\,\Delta\boldsymbol\eta \cdot \mathbf r}\, d^2r,$$0, defining the aliasing-free region (AFR) in which ambiguous grating lobes do not appear. For canonical arrays:

• For a Uniform Linear Array (ULA), the AFR is an "eye-shaped" region described in closed form in the infinite-aperture limit. Its width and height are set by the intersection of the bandlimit $$\Xi(\Delta\mathbf x,\, \Delta\boldsymbol\eta) = \int_R h(\mathbf r + \tfrac{1}{2} \Delta\mathbf x)\, h^*(\mathbf r - \tfrac{1}{2} \Delta\mathbf x)\, e^{-j\,\Delta\boldsymbol\eta \cdot \mathbf r}\, d^2r,$$1 and the sampling interval $$\Xi(\Delta\mathbf x,\, \Delta\boldsymbol\eta) = \int_R h(\mathbf r + \tfrac{1}{2} \Delta\mathbf x)\, h^*(\mathbf r - \tfrac{1}{2} \Delta\mathbf x)\, e^{-j\,\Delta\boldsymbol\eta \cdot \mathbf r}\, d^2r,$$2.
• For a Uniform Circular Array (UCA), the AFR is a disk of radius $$\Xi(\Delta\mathbf x,\, \Delta\boldsymbol\eta) = \int_R h(\mathbf r + \tfrac{1}{2} \Delta\mathbf x)\, h^*(\mathbf r - \tfrac{1}{2} \Delta\mathbf x)\, e^{-j\,\Delta\boldsymbol\eta \cdot \mathbf r}\, d^2r,$$3, where $$\Xi(\Delta\mathbf x,\, \Delta\boldsymbol\eta) = \int_R h(\mathbf r + \tfrac{1}{2} \Delta\mathbf x)\, h^*(\mathbf r - \tfrac{1}{2} \Delta\mathbf x)\, e^{-j\,\Delta\boldsymbol\eta \cdot \mathbf r}\, d^2r,$$4 encodes the angular geometry (Monnoyer et al., 9 Dec 2025).

6. Trade-offs in Aperture Design, Thinning, and Aliasing

Increasing the array element spacing $$\Xi(\Delta\mathbf x,\, \Delta\boldsymbol\eta) = \int_R h(\mathbf r + \tfrac{1}{2} \Delta\mathbf x)\, h^*(\mathbf r - \tfrac{1}{2} \Delta\mathbf x)\, e^{-j\,\Delta\boldsymbol\eta \cdot \mathbf r}\, d^2r,$$5 narrows the AFR; grating lobes appear once the maximum local spatial frequency $$\Xi(\Delta\mathbf x,\, \Delta\boldsymbol\eta) = \int_R h(\mathbf r + \tfrac{1}{2} \Delta\mathbf x)\, h^*(\mathbf r - \tfrac{1}{2} \Delta\mathbf x)\, e^{-j\,\Delta\boldsymbol\eta \cdot \mathbf r}\, d^2r,$$6 exceeds $$\Xi(\Delta\mathbf x,\, \Delta\boldsymbol\eta) = \int_R h(\mathbf r + \tfrac{1}{2} \Delta\mathbf x)\, h^*(\mathbf r - \tfrac{1}{2} \Delta\mathbf x)\, e^{-j\,\Delta\boldsymbol\eta \cdot \mathbf r}\, d^2r,$$7. Enlarging aperture length or radius (increasing $$\Xi(\Delta\mathbf x,\, \Delta\boldsymbol\eta) = \int_R h(\mathbf r + \tfrac{1}{2} \Delta\mathbf x)\, h^*(\mathbf r - \tfrac{1}{2} \Delta\mathbf x)\, e^{-j\,\Delta\boldsymbol\eta \cdot \mathbf r}\, d^2r,$$8 or $$\Xi(\Delta\mathbf x,\, \Delta\boldsymbol\eta) = \int_R h(\mathbf r + \tfrac{1}{2} \Delta\mathbf x)\, h^*(\mathbf r - \tfrac{1}{2} \Delta\mathbf x)\, e^{-j\,\Delta\boldsymbol\eta \cdot \mathbf r}\, d^2r,$$9) raises the maximum frequency and further shrinks the AFR, even as it improves point-target resolution. The design of sparse (thinned) or extremely large-scale arrays requires careful balance between the desired resolution and the tolerable region of ambiguity-free operation (Monnoyer et al., 9 Dec 2025).

Schematic guidelines:

• To maximize resolution in $\Xi(\mathbf 0, \mathbf 0) = 1$0, use large $\Xi(\mathbf 0, \mathbf 0) = 1$1 or a scene with large background spectral support $\Xi(\mathbf 0, \mathbf 0) = 1$2.
• To maximize angular resolution, use large $\Xi(\mathbf 0, \mathbf 0) = 1$3 (increasing $\Xi(\mathbf 0, \mathbf 0) = 1$4) or ensure the scene provides high spatial position diversity ($\Xi(\mathbf 0, \mathbf 0) = 1$5 small).
• For synthetic-aperture applications, the time-growing $\Xi(\mathbf 0, \mathbf 0) = 1$6 improves angular resolution only if the correlation function $\Xi(\mathbf 0, \mathbf 0) = 1$7 remains stationary.
• In the near field, nontrivial coupling between linear and angular estimates requires either multi-frequency or multi-baseline measurements to separate effects.

7. Illustrative Examples and Comparison to Classical Far-Field AF

Canonical cases clarify essential SAF behaviors:

• In the Fraunhofer (far-field) limit, the spatial-frequency shift is constant; the SAF reduces to a one-dimensional sinc response in angle, with straight and periodic grating lobes appearing for $\Xi(\mathbf 0, \mathbf 0) = 1$8.
• In the Fresnel (near-field) regime, SAF mainlobe width is broadened and the peak shifts along the line $\Xi(\mathbf 0, \mathbf 0) = 1$9, reflecting linear–angular estimate coupling and chirp-like spectral properties.
• Thinned or sparse arrays manifest curved aliasing fronts in the SAF, as predicted quantitatively by the soft bandlimit $$\Xi(\Delta x,\; \Delta k) = \int_{-L/2}^{+L/2} h(x + \tfrac{\Delta x}{2})\, h^*(x - \tfrac{\Delta x}{2})\, e^{-j \Delta k x}\, dx.$$0 and allegiance to the AFR (Monnoyer et al., 9 Dec 2025).

A summary table encapsulates key differences:

Regime	SAF Form	Grating Lobe Structure
Far-field	Sinc in angle (linear spatial freq.)	Straight, periodic; $$\Xi(\Delta x,\; \Delta k) = \int_{-L/2}^{+L/2} h(x + \tfrac{\Delta x}{2})\, h^*(x - \tfrac{\Delta x}{2})\, e^{-j \Delta k x}\, dx.$$1
Near-field	Chirped/shifted-sinc, mainlobe at $$\Xi(\Delta x,\; \Delta k) = \int_{-L/2}^{+L/2} h(x + \tfrac{\Delta x}{2})\, h^*(x - \tfrac{\Delta x}{2})\, e^{-j \Delta k x}\, dx.$$2	Curved, array- and geometry-dependent
Thinned array (NF)	Discrete-spectrum, aliasing controlled by $$\Xi(\Delta x,\; \Delta k) = \int_{-L/2}^{+L/2} h(x + \tfrac{\Delta x}{2})\, h^*(x - \tfrac{\Delta x}{2})\, e^{-j \Delta k x}\, dx.$$3	Multiple, non-uniform, eye/disk AFR

The transition from far- to near-field involves a shift from constant to locally-varying spatial frequency, with corresponding changes in resolution, ambiguity, and grating lobe morphology. By leveraging closed-form expressions for $$\Xi(\Delta x,\; \Delta k) = \int_{-L/2}^{+L/2} h(x + \tfrac{\Delta x}{2})\, h^*(x - \tfrac{\Delta x}{2})\, e^{-j \Delta k x}\, dx.$$4 and the AFRs, one can predict and control these phenomena in array system design (Yuryev et al., 2014, Monnoyer et al., 9 Dec 2025).