DSP in Phased Array 3D Imaging

Chris Bore · Watch Now · DSP Online Conference 2020 · Duration: 29:03

Phase is often the poor relation to magnitude in Fourier Transform analysis but can be the crucial measurement in some applications. Such is the case with phased array methods for imaging: in radar, security and medical imaging applications. The key to such applications is physics modelling and measurement of phase to reconstruct wave propagation and hence reconstruct 3D or 4D images of reflecting, refracting or scattering objects. This talk offers an overview of phased array imaging from the viewpoint of DSP, modelled as a 3D Fourier or Laplace Transform problem.

What this presentation is about and why it matters

This talk explains how phased-array measurements of scattered waves can be turned into three‑dimensional images using digital signal processing. Instead of treating magnitude as the dominant observable, the talk emphasizes phase and the spatial structure of waves (the so‑called k‑space or spatial frequency). Using radio‑wave breast imaging as a motivating example, the speaker shows how measurements made with an antenna array and a vector network analyzer (VNA) form a frequency‑domain sampling of the scene that can be interpreted and inverted mathematically as a 3D Fourier (or, in more general settings, Laplace) transform.

Why this matters for engineers: many practical imaging systems (radar, sonar, medical microwave imaging, Wi‑Fi sensing, lidar, ultrasound) measure complex-valued responses that include phase. Treating those measurements correctly — often directly in the frequency domain — leads to more physically faithful reconstructions, avoids common DSP pitfalls (for example, naïve zero‑padding tricks that confuse time and space), and yields better designs for acquisition and processing.

Who will benefit the most from this presentation

  • DSP engineers who design and implement imaging and inverse‑problem algorithms.
  • RF and antenna engineers working with S‑parameters, VNAs and phased arrays.
  • Students of radar, sonar, or medical imaging who want intuition linking wave physics and Fourier/DSP methods.
  • System designers deciding what to measure (time‑domain pulses vs. frequency sweeps) and how to process limited, noisy data.
  • Researchers interested in practical considerations: sampling in k‑space, multipath, and refractive effects.

What you need to know

The presentation is accessible if you have a basic background in the following areas. Below are the key concepts and a few compact mathematical pointers to help you follow the talk more easily.

  • Complex phasors and waves: A time‑harmonic travelling wave is conveniently written as a phasor. You will see expressions like $E\,e^{j(\omega t - kD)}$, where $\omega$ is angular frequency, $k$ is the spatial wave number and $D$ is distance along the propagation path. Understanding complex exponentials and Euler's identity will help.
  • Wave number (k): The wave number measures phase change per unit distance: $k = 2\pi/\lambda$. Think of $k$ as the spatial analog of angular frequency; it indexes oscillation in space rather than time.
  • S‑parameters and VNAs: A VNA measures complex ratios between ports. An S_{M,N}(\omega) measurement at one frequency is a complex phasor that summarizes amplitude and phase of the combined paths from transmitter N to receiver M.
  • Frequency‑domain vs time‑domain: The VNA gives direct frequency‑domain samples. A pulse in time corresponds to a wideband spectrum; conversely, sparse frequency samples produce poor native time resolution. Practitioners often transform to time via an IDFT, but the talk argues for direct frequency‑domain (k‑space) processing for imaging.
  • Imaging as an inverse Fourier problem: Under common simplifying assumptions (single scattering, isotropic scatterers, straight rays, homogeneous medium) the measured scattering matrix can be written as a spatial Fourier transform of the scattering coefficient field $s(P)$. Reconstruction then resembles applying an inverse 3D Fourier transform (or a Laplace inversion for more general path models): schematically $S_{M,N}(\omega)\approx\int s(P)\,e^{-j k D(P)}\,dV$ and inversion aims to recover $s(P)$.
  • Practical DSP issues: sampling in angle/frequency (k‑space coverage), zero‑padding and interpolation, aperture geometry (hemisphere vs planar array), and the effect of refractive index variations and multipath on phase are all important.

Glossary

  • Phase — the angular part of a complex phasor; carries spatial position information for waves.
  • Phasor — complex representation of a sinusoid, $A\,e^{j\phi}$, used to track amplitude and phase.
  • Wave number (k) — spatial frequency, $k=2\pi/\lambda$, giving radians of phase per meter.
  • Scattering coefficient s(P) — a local property at point P that determines how strongly that point scatters incident energy.
  • S‑parameters — measured complex ratios (received/transmitted) between antenna ports as a function of frequency.
  • Vector Network Analyzer (VNA) — instrument that sweeps frequency and measures complex transfer between ports (amplitude and phase).
  • k‑space — the domain of spatial frequencies (wave numbers); sampling here controls image resolution and coverage.
  • Inverse problem — recovering the unknown scattering field from measured S‑parameters (an ill‑posed or underdetermined problem unless constrained).
  • Zero‑padding — adding zeros in the frequency domain before an IDFT; it interpolates the time‑domain result but does not create new physical information.
  • Multipath — when waves take multiple routes to the receiver; leads to sums of phasors with different path lengths and phases.

A few nice words about the presentation

The speaker blends physics, RF practice and DSP intuition in a clear, practitioner‑oriented way. He emphasizes physical meaning over formal algebra where helpful, yet is willing to use compact mathematics to reveal structure (for example, the connection between measured S‑parameters and a 3D Fourier/Laplace transform). If you want a concrete bridge between what a VNA measures and how to turn those complex spectra into a spatial image — including practical cautions about sampling, zero‑padding and assumptions you must check — this talk is a concise, honest and visually minded guide.

Formatting help
italicssurround text with
*asterisks*
boldsurround text with
**two asterisks**
hyperlink
[hyperlink](https://example.com)
or just a bare URL
code
surround text with
`backticks`
strikethroughsurround text with
~~two tilde characters~~
quote
prefix with
>

Unknown
Score: 0 | 5 years ago | no reply

Good and informative
myinsite

Ilias
Score: 0 | 5 years ago | no reply

Very nice perspective! Excellent!

Danilo
Score: 0 | 5 years ago | 1 reply

Nice one.

ChrisBoreSpeaker
Score: 0 | 5 years ago | no reply

Glad you enjoyed it

Brewster
Score: 1 | 5 years ago | 2 replies

Hi Chris, You need come out of retirement :-) Though you obviously solved the resolution problem, would a frequency chirp instead of fixed frequencies have any benefits or issues?

ChrisBoreSpeaker
Score: 0 | 5 years ago | no reply

Also, the VNA approach preserves all the data, despite the high frequencies, and so we can go back and do different reconstructions on processing on the raw data set, which you couldn't practically do with time domain because you normally reduce the data

ChrisBoreSpeaker
Score: 0 | 5 years ago | no reply

Hello Brewster, nice to see you!
The chirp might allow faster data gathering but it loses the direct frequency domain data gathering: which does not suffer from Gibbs phenomenon and all that, because we actually do measure each frequency in CW.