Range-Doppler Processing

The Range-Doppler Map as a 2D Image

A pulsed-Doppler radar transmits NpN_p pulses with pulse repetition interval (PRI) TPRIT_{\text{PRI}}. The received echoes form a 2D data matrix β€” fast time (within each pulse, encoding range) and slow time (across pulses, encoding velocity). Processing this matrix with a 2D FFT produces the range-Doppler map, which is the simplest "image" a radar can form.

For RF imaging, the range-Doppler map is the matched-filter image AHy\mathbf{A}^{H}\mathbf{y} restricted to the range and Doppler dimensions.

Definition:

Fast-Time / Slow-Time Data Matrix

The radar collects NpN_p pulses, each containing NsN_s range samples. The data matrix is:

X∈CNsΓ—Np,\mathbf{X} \in \mathbb{C}^{N_s \times N_p},

where:

  • Column nn contains the nn-th pulse after demodulation.
  • Row mm corresponds to fast-time sample mm (range bin mm).
  • Element Xm,nX_{m,n} contains the signal from range bin mm, pulse nn.

Fast-time processing (across columns) resolves range via matched filtering. Slow-time processing (across rows) resolves Doppler via coherent integration (DFT).

Definition:

Range-Doppler Processing via 2D FFT

The range-Doppler map is formed by:

  1. Fast-time matched filtering: For each pulse nn, compute yn(Ο„)=rn(Ο„)βˆ—sβˆ—(βˆ’Ο„)y_n(\tau) = r_n(\tau) * s^*(-\tau) (or equivalently multiply in frequency and IFFT). This resolves targets in range.

  2. Slow-time DFT: For each range bin mm, compute the DFT across the NpN_p pulses: Zm[k]=βˆ‘n=0Npβˆ’1yn[m] eβˆ’j2Ο€kn/Np,Z_m[k] = \sum_{n=0}^{N_p - 1} y_n[m] \, e^{-j2\pi kn/N_p}, resolving targets in Doppler.

The result ∣Zm[k]∣2|Z_m[k]|^2 is the range-Doppler map with:

  • Range resolution: Ξ”R=c/(2W)\Delta R = c/(2W).
  • Doppler resolution: Ξ”fd=1/(NpTPRI)=1/TCPI\Delta f_d = 1/(N_p T_{\text{PRI}}) = 1/T_{\text{CPI}}.
  • Velocity resolution: Ξ”v=Ξ»/(2TCPI)\Delta v = \lambda/(2T_{\text{CPI}}).

Theorem: Coherent Integration Gain

Coherent integration of NpN_p pulses via the slow-time DFT increases the output SNR\text{SNR} by a factor of NpN_p relative to the single-pulse SNR\text{SNR}:

SNRcoh=Npβ‹…SNR1.\text{SNR}_{\text{coh}} = N_p \cdot \text{SNR}_{1}.

The integrated SNR\text{SNR} is:

SNRcoh=NpPtG2Ξ»2Οƒ(4Ο€)3R4kBTsWL.\text{SNR}_{\text{coh}} = \frac{N_p P_t G^2 \lambda^{2} \sigma}{(4\pi)^3 R^4 k_B T_s W L}.

Each pulse adds signal coherently (amplitude adds) while noise adds incoherently (power adds). Signal power grows as Np2N_p^2 while noise power grows as NpN_p, giving SNR\text{SNR} gain NpN_p.

Definition:

PRF Design and Ambiguity Trade-Offs

The pulse repetition frequency (PRF) fPRF=1/TPRIf_{\text{PRF}} = 1/T_{\text{PRI}} determines both the unambiguous range and unambiguous velocity:

  • Unambiguous range: Rua=c/(2fPRF)R_{\text{ua}} = c/(2 f_{\text{PRF}}).
  • Unambiguous velocity: vua=Ξ»fPRF/4v_{\text{ua}} = \lambda f_{\text{PRF}} / 4.

These are inversely related: Ruaβ‹…vua=cΞ»/8R_{\text{ua}} \cdot v_{\text{ua}} = c \lambda / 8.

This leads to three PRF regimes:

  • Low PRF: Unambiguous in range, ambiguous in velocity.
  • High PRF: Unambiguous in velocity, ambiguous in range.
  • Medium PRF: Ambiguous in both, but with manageable unfolding.

Example: Constructing a Range-Doppler Map

A pulsed radar operates at f0=10f_0 = 10 GHz (Ξ»=3\lambda = 3 cm), W=10W = 10 MHz, PRF =1= 1 kHz, Np=64N_p = 64 pulses. Two targets exist: Target A at RA=15R_A = 15 km, vA=30v_A = 30 m/s; Target B at RB=15.01R_B = 15.01 km, vB=βˆ’50v_B = -50 m/s. (a) Can the radar resolve these targets in range? (b) What are their Doppler frequencies? (c) Can the radar resolve them in Doppler?

Range-Doppler Map Visualization

Generate a range-Doppler map for configurable targets. Observe how bandwidth determines range resolution, CPI determines velocity resolution, and PRF determines ambiguity boundaries.

Parameters
20
64
1
15

Definition:

Range Migration

During the coherent processing interval (CPI) TCPI=NpTPRIT_{\text{CPI}} = N_p T_{\text{PRI}}, a target at radial velocity vv migrates through range bins:

Ξ”Rmig=vβ‹…TCPI.\Delta R_{\text{mig}} = v \cdot T_{\text{CPI}}.

When Ξ”Rmig>Ξ”R\Delta R_{\text{mig}} > \Delta R (i.e., migration exceeds one range bin), the target's energy smears across range bins, degrading coherent integration gain.

Range migration correction compensates this drift before Doppler processing. Common techniques include:

  • Keystone transform: resampling the data matrix to align range cells across pulses.
  • Range walk compensation: linear range shift per pulse based on estimated velocity.

Common Mistake: Forgetting Range-Doppler Coupling in LFM

Mistake:

Treating the range-Doppler map as having independent range and Doppler axes when using an LFM waveform.

Correction:

For LFM, a Doppler shift fdf_d causes an apparent range shift Ξ”R=fdc/(2ΞΌ)=fdcT/(2W)\Delta R = f_d c / (2\mu) = f_d cT / (2W). A target at 100 m/s with Ξ»=3\lambda = 3 cm and W/T=1012W/T = 10^{12} Hz/s appears shifted by Ξ”R=fd/(2ΞΌ)β‹…c\Delta R = f_d/(2\mu) \cdot c. For slow targets this is negligible; for fast targets or long CPIs, it must be corrected by the keystone transform or by adjusting the range-Doppler map axes.

Coherent Processing Interval (CPI)

The time duration TCPI=NpTPRIT_{\text{CPI}} = N_p T_{\text{PRI}} over which pulses are coherently integrated. Determines velocity resolution Ξ”v=Ξ»/(2TCPI)\Delta v = \lambda/(2T_{\text{CPI}}) and total integration gain NpN_p.

Quick Check

A radar at λ=3\lambda = 3 cm needs to detect targets at up to R=150R = 150 km with unambiguous velocity. What is the minimum PRF if the maximum target velocity is vmax⁑=300v_{\max} = 300 m/s?

20 kHz

40 kHz

1 kHz

10 kHz

Key Takeaway

The 2D FFT (fast-time matched filter + slow-time DFT) produces the range-Doppler map with resolution Ξ”R=c/(2W)\Delta R = c/(2W) in range and Ξ”v=Ξ»/(2TCPI)\Delta v = \lambda/(2T_{\text{CPI}}) in velocity. Coherent integration provides NpN_p gain in SNR\text{SNR}. Range migration must be corrected when vβ‹…TCPI>Ξ”Rv \cdot T_{\text{CPI}} > \Delta R.