Near-Field Diffraction Tomography

Near-Field Diffraction Tomography

The Fourier diffraction theorem as presented in Section s02 assumes far-field measurements -- observation points many wavelengths from the object, where wavefronts are approximately planar. In many RF imaging scenarios -- medical microwave, through-wall radar, and especially extremely large MIMO (XL-MIMO) arrays -- the aperture is comparable to the target range, placing the system squarely in the near field.

This section examines near-field phase corrections, Fresnel-zone analysis, and the connection to near-field beamfocusing in communications (Telecom Ch 33). The near field offers both a challenge (invalid plane-wave assumption) and an opportunity (access to evanescent wave information for sub-wavelength resolution).

Definition:

Near-Field and Far-Field Regions

For an aperture (or object) of diameter DD operating at wavelength Ξ»\lambda:

  • Far field (Fraunhofer region): r>dff=2D2/Ξ»r > d_{\mathrm{ff}} = 2D^2/\lambda. Wavefronts are approximately planar; the FDT applies directly.

  • Fresnel (radiating near-field) region: 0.62D3/Ξ»<r<2D2/Ξ»0.62\sqrt{D^3/\lambda} < r < 2D^2/\lambda. Wavefronts are spherical; phase curvature is significant.

  • Reactive near field: r<0.62D3/Ξ»r < 0.62\sqrt{D^3/\lambda}. Evanescent fields carry spatial frequencies beyond 2k02k_0, enabling sub-wavelength resolution.

Example (XL-MIMO): A 2 m aperture at 28 GHz (Ξ»=1.07\lambda = 1.07 cm): dff=2(2)2/0.0107=747d_{\mathrm{ff}} = 2(2)^2/0.0107 = 747 m. Any target within several hundred meters is in the Fresnel zone.

,

Definition:

Near-Field Forward Model

In the near field, the incident wave from a point source at st\mathbf{s}_{t} is spherical:

Einc(r)=ejk0∣rβˆ’st∣4Ο€βˆ£rβˆ’st∣E_{\text{inc}}(\mathbf{r}) = \frac{e^{jk_0|\mathbf{r} - \mathbf{s}_{t}|}}{4\pi|\mathbf{r} - \mathbf{s}_{t}|}

(not a plane wave ejk0k^iβ‹…re^{jk_0\hat{\mathbf{k}}_i\cdot\mathbf{r}}). The scattered field at receiver rr\mathbf{r}_{r} becomes:

Es(rr)=k02∫G(rr,rβ€²) χ(rβ€²) G(rβ€²,st) drβ€².E_s(\mathbf{r}_{r}) = k_0^2 \int G(\mathbf{r}_{r}, \mathbf{r}')\,\chi(\mathbf{r}')\, G(\mathbf{r}', \mathbf{s}_{t})\,d\mathbf{r}'.

This is a double convolution with the Green's function, not a single Fourier relationship. The Fourier diffraction theorem does not directly apply; instead, we must either (1) correct the near-field data before applying the FDT, or (2) use a near-field-aware reconstruction algorithm.

Theorem: Near-Field to Far-Field Transformation

For a planar measurement surface at distance dd from the object, the near-field data can be transformed to equivalent far-field data via:

E~ff(kβŠ₯)=E~nf(kβŠ₯) eβˆ’jkzd,\tilde{E}_{\text{ff}}(\mathbf{k}_\perp) = \tilde{E}_{\text{nf}}(\mathbf{k}_\perp)\,e^{-jk_z d},

where kz=k02βˆ’βˆ£kβŠ₯∣2k_z = \sqrt{k_0^2 - |\mathbf{k}_\perp|^2} and E~\tilde{E} denotes the spatial Fourier transform over the measurement aperture.

After this transformation, the standard FDT-based reconstruction can be applied.

Each plane-wave component of the field propagates with a phase shift ejkzde^{jk_z d} over distance dd. The NF-FF transformation removes this propagation phase, back-propagating the measured field to the object plane.

Proof
Plane-wave decomposition

The near-field data Enf(ρ,d)E_{\text{nf}}(\boldsymbol{\rho}, d) is decomposed into plane-wave components:

Enf(ρ,d)=∫E~nf(kβŠ₯) ejkβŠ₯⋅ρ dkβŠ₯.E_{\text{nf}}(\boldsymbol{\rho}, d) = \int \tilde{E}_{\text{nf}}(\mathbf{k}_\perp)\, e^{j\mathbf{k}_\perp \cdot \boldsymbol{\rho}}\, d\mathbf{k}_\perp.

Remove propagation phase

Each component has accumulated phase ejkzde^{jk_z d} during propagation from the object plane to the measurement plane. Multiplying by eβˆ’jkzde^{-jk_z d} removes this phase:

E~ff(kβŠ₯)=E~nf(kβŠ₯) eβˆ’jkzd.\tilde{E}_{\text{ff}}(\mathbf{k}_\perp) = \tilde{E}_{\text{nf}}(\mathbf{k}_\perp)\,e^{-jk_z d}.

The result is the field spectrum at the object plane, which the far-field FDT can then invert. β– \quad\blacksquare

Limitation for evanescent components

For ∣kβŠ₯∣>k0|\mathbf{k}_\perp| > k_0, kzk_z is imaginary and eβˆ’jkzd=e+Ξ³de^{-jk_z d} = e^{+\gamma d} with Ξ³=∣kβŠ₯∣2βˆ’k02\gamma = \sqrt{|\mathbf{k}_\perp|^2 - k_0^2}. This exponential amplification makes the transformation ill-conditioned for evanescent data -- regularization or truncation is needed.

Fresnel Phase Correction for Imaging

In the Fresnel zone, the dominant near-field effect is phase curvature: the spherical wavefront from a source at distance r0r_0 introduces a quadratic phase error across the aperture:

Δϕ(ρ)β‰ˆk0∣ρ∣22r0\Delta\phi(\boldsymbol{\rho}) \approx \frac{k_0 |\boldsymbol{\rho}|^2}{2 r_0}

where ρ\boldsymbol{\rho} is the transverse displacement from the aperture center. For an aperture of diameter DD:

Δϕmax⁑=k0D28r0=Ο€D24Ξ»r0.\Delta\phi_{\max} = \frac{k_0 D^2}{8 r_0} = \frac{\pi D^2}{4 \lambda r_0}.

When Δϕmax⁑β‰ͺΟ€/4\Delta\phi_{\max} \ll \pi/4 (the Fraunhofer condition), the phase curvature is negligible. When Δϕmaxβ‘βˆΌΟ€\Delta\phi_{\max} \sim \pi (Fresnel zone), the curvature must be compensated in the reconstruction.

Compensation: Apply the conjugate quadratic phase eβˆ’jk0∣ρ∣2/(2r0)e^{-jk_0|\boldsymbol{\rho}|^2/(2r_0)} to each measurement before standard FDT processing. This is equivalent to beamfocusing at range r0r_0 rather than beamsteering to a far-field direction.

Near-Field vs. Far-Field PSF Comparison

Compares the point spread function (PSF) for near-field and far-field imaging models at the selected distance and frequency.

Left: PSF using the far-field (plane-wave) model. Right: PSF using the near-field (spherical-wave) model.

As the measurement distance decreases below dffd_{\mathrm{ff}}, the far-field PSF degrades (broadens, develops asymmetric sidelobes) while the near-field PSF remains sharp. The dashed circle shows the βˆ’3-3 dB mainlobe width.

Parameters
1
10
0.3

Example: Near-Field Imaging with XL-MIMO

An XL-MIMO array with D=2D = 2 m aperture operates at 28 GHz (Ξ»=1.07\lambda = 1.07 cm). The target is at range r=50r = 50 m.

(a) Is the target in the near field or far field? (b) What is the Fresnel phase error across the aperture? (c) What happens if we ignore the near-field correction?

Solution
Far-field distance

dff=2D2/Ξ»=2(2)2/0.0107=747d_{\mathrm{ff}} = 2D^2/\lambda = 2(2)^2/0.0107 = 747 m. The target at 50 m is deep in the Fresnel zone (r/dff=0.067r/d_{\mathrm{ff}} = 0.067).

Fresnel phase error

Δϕmax⁑=Ο€D2/(4Ξ»r)=Ο€(2)2/(4Γ—0.0107Γ—50)=5.87\Delta\phi_{\max} = \pi D^2/(4\lambda r) = \pi (2)^2/(4 \times 0.0107 \times 50) = 5.87 rad β‰ˆ1.87Ο€\approx 1.87\pi. This is much larger than Ο€/4\pi/4 -- the near-field correction is essential.

Impact of ignoring correction

Without near-field phase compensation, the reconstructed image suffers from:

  • Position errors: Targets appear shifted by Ξ”xβ‰ˆD2/(8r)=0.5\Delta x \approx D^2/(8r) = 0.5 cm at the edge of the field of view.
  • Defocusing: The PSF mainlobe broadens by a factor βˆΌΞ”Ο•max⁑/Ο€β‰ˆ1.9\sim \Delta\phi_{\max}/\pi \approx 1.9, roughly halving the achievable resolution.
  • Sidelobe distortion: The PSF develops asymmetric sidelobes that create ghost targets.

Definition:

Sub-Wavelength Resolution from Evanescent Waves

In the reactive near field (r≲λr \lesssim \lambda), evanescent waves carry spatial frequency information beyond 2k02k_0. The evanescent field decays as:

Eevan∝eβˆ’βˆ£kβŠ₯∣2βˆ’k02 d.E_{\text{evan}} \propto e^{-\sqrt{|\mathbf{k}_\perp|^2 - k_0^2}\,d}.

A spatial frequency ∣kβŠ₯∣|\mathbf{k}_\perp| is detectable if the evanescent signal exceeds the noise floor:

∣kβŠ₯∣<k02+(ln⁑(SNR)d)2,|\mathbf{k}_\perp| < \sqrt{k_0^2 + \left(\frac{\ln(\text{SNR})}{d}\right)^2},

giving enhanced resolution:

Ξ΄nf=Ο€Kmax⁑nf<Ξ»4.\delta_{\text{nf}} = \frac{\pi}{K_{\max}^{\text{nf}}} < \frac{\lambda}{4}.

At d=Ξ»/10d = \lambda/10 with SNR=40\text{SNR} = 40 dB: Kmax⁑nfβ‰ˆ12k0K_{\max}^{\text{nf}} \approx 12k_0, giving Ξ΄β‰ˆΞ»/24\delta \approx \lambda/24.

Connection to Near-Field Beamfocusing (Telecom Ch 33)

In communications, XL-MIMO systems exploit near-field propagation for beamfocusing: directing energy to a specific point in space (range + angle), not just a direction. The beamfocusing weights are exactly the conjugate of the near-field steering vector:

vfocus=[…,β€…β€Šeβˆ’jk0∣rnβˆ’p0∣,β€…β€Šβ€¦]Tβˆ₯[…]βˆ₯\mathbf{v}_{\text{focus}} = \frac{[\ldots,\; e^{-jk_0|\mathbf{r}_{n} - \mathbf{p}_{0}|},\; \ldots]^T}{\|[\ldots]\|}

In imaging, the same operation appears as near-field matched filtering -- focusing the array response to each voxel p\mathbf{p} in the scene. The near-field imaging problem and the near-field beamfocusing problem are mathematically dual: imaging asks "what is at each point?" while beamfocusing asks "how do I concentrate energy at one point?" Both require the same spherical-wave phase model.

Common Mistake: Near-Field Calibration Is Critical

Mistake:

Near-field imaging amplifies modeling errors because the spherical-wave phase varies rapidly with position:

  • Position errors: A 1 mm error in sensor position at 10 GHz introduces phase error Δϕ=2Ο€(0.001)/(0.03)β‰ˆ0.21\Delta\phi = 2\pi(0.001)/(0.03) \approx 0.21 rad.

  • Mutual coupling: Closely spaced near-field sensors couple electromagnetically, distorting the measured fields.

  • Background subtraction: The incident field must be precisely subtracted to isolate EsE_s; near-field incident fields are harder to model than plane waves.

Correction:

Calibrate by imaging known reference objects (point scatterers or dielectric cylinders) and adjusting the forward model to match the measured data. Autofocus algorithms can also estimate and correct position errors from the data itself.

Common Mistake: Evanescent Wave Recovery Is Ill-Conditioned

Mistake:

The NF-FF transformation amplifies evanescent components by eΞ³de^{\gamma d} where Ξ³=∣kβŠ₯∣2βˆ’k02\gamma = \sqrt{|\mathbf{k}_\perp|^2 - k_0^2}. At ∣kβŠ₯∣=2k0|\mathbf{k}_\perp| = 2k_0 and d=Ξ»d = \lambda: Ξ³d=3β‹…2Ο€β‰ˆ10.9\gamma d = \sqrt{3} \cdot 2\pi \approx 10.9, giving amplification e10.9β‰ˆ54,000e^{10.9} \approx 54{,}000.

Any measurement noise in the evanescent band is amplified by the same factor, overwhelming the signal.

Correction:

Apply a spectral cutoff or Tikhonov regularization in the NF-FF transformation to suppress evanescent components beyond the SNR-limited bandwidth. Accept the resolution limit imposed by the measurement distance and SNR.

πŸ”§Engineering Note

Near-Field Imaging with Commercial XL-MIMO

5G/6G XL-MIMO panels with 256--1024 elements at 28 GHz or mmWave frequencies have apertures of 0.5--2 m, creating far-field distances of 50--750 m. This means that all indoor targets and most outdoor targets of interest are in the Fresnel zone.

Implications for joint communication-sensing: The same XL-MIMO array used for beamfocusing in communications can perform near-field imaging of the environment. The near-field phase model is needed for both applications.

Practical challenge: Computing the near-field steering vector for each voxel requires O(Nantβ‹…Q)O(N_{\text{ant}} \cdot Q) operations, where QQ is the number of voxels. For a 1024-element array imaging a 100Γ—100Γ—100100 \times 100 \times 100 scene: 10910^9 operations per frequency -- feasible on GPU but not in real time on CPU.

Near-Field vs. Far-Field Imaging

PropertyFar-FieldNear-Field (Fresnel)Near-Field (Reactive)
Distance criterionr>2D2/Ξ»r > 2D^2/\lambda0.62D3/Ξ»<r<2D2/Ξ»0.62\sqrt{D^3/\lambda} < r < 2D^2/\lambdar<0.62D3/Ξ»r < 0.62\sqrt{D^3/\lambda}
Wavefront modelPlanarSphericalSpherical + evanescent
FDT applicable?Yes, directlyAfter NF-FF correctionNo (need full near-field model)
Resolution limitΞ»/4\lambda/4Ξ»/4\lambda/4 (same propagating content)<Ξ»/4< \lambda/4 (evanescent gain)
CalibrationStandardCriticalExtremely critical
Typical RF scenarioRadar at long rangeIndoor MIMO, medicalNanoscale probing

Quick Check

A 50 cm aperture array operates at 10 GHz (Ξ»=3\lambda = 3 cm). A target is at 5 m range. Is the target in the near field or far field?

Far field, because 5Β m≫λ5\text{ m} \gg \lambda

Near field, because 5Β m<2D2/Ξ»=16.7Β m5\text{ m} < 2D^2/\lambda = 16.7\text{ m}

Far field, because the Fresnel number is less than 1

Cannot be determined without knowing the target size

Correction:
Near field, because 5Β m<2D2/Ξ»=16.7Β m5\text{ m} < 2D^2/\lambda = 16.7\text{ m}

dff=2(0.5)2/0.03=16.7d_{\text{ff}} = 2(0.5)^2/0.03 = 16.7 m. At 5 m, the target is in the Fresnel zone.

Quick Check

Evanescent waves in near-field imaging enable:

Longer measurement range

Sub-wavelength resolution beyond the diffraction limit

Higher SNR than propagating waves

Wider bandwidth without additional frequencies

Correction:
Sub-wavelength resolution beyond the diffraction limit

Evanescent waves carry spatial frequencies ∣kβŠ₯∣>k0|\mathbf{k}_\perp| > k_0, providing resolution finer than Ξ»/4\lambda/4.

Historical Note: Near-Field Microwave Imaging for Medicine

1990s--present

Near-field microwave imaging for medical diagnostics began in the 1990s with Paul Meaney and Keith Paulsen at Dartmouth, who built the first clinical prototype for breast cancer detection. Their system used a 16-element circular array at 300--900 MHz, imaging at ranges of 5--10 cm -- deep in the near field. The key insight was that near-field phase modeling, while computationally expensive, was essential for accurate reconstruction.

Today, near-field microwave imaging is an active research area for breast cancer screening (as a non-ionizing alternative to mammography), brain stroke detection, and bone quality assessment. The same near-field algorithms are now being adapted for XL-MIMO sensing in 6G communications.

Fresnel Zone

The region of space between the reactive near field and the Fraunhofer far field, where wavefronts are spherical but evanescent fields are negligible. In the Fresnel zone, the dominant near-field effect is quadratic phase curvature across the aperture.

Related: Fraunhofer Region (Far Field)

Fraunhofer Region (Far Field)

The region beyond r=2D2/Ξ»r = 2D^2/\lambda where wavefronts from an aperture of diameter DD are approximately planar. In the far field, the Fourier diffraction theorem applies directly without near-field corrections.

Related: Fresnel Zone

Key Takeaway

Most practical RF imaging systems operate in the near field where spherical-wave modeling is essential. The NF-FF transformation converts near-field data to far-field equivalent for FDT processing, but is ill-conditioned for evanescent components. XL-MIMO arrays have far-field distances of hundreds of meters, making near-field imaging the default regime. Near-field imaging and near-field beamfocusing (communications) are mathematically dual problems. Evanescent waves enable sub-wavelength resolution but only at very close range (d≲λd \lesssim \lambda).

Multi-Frequency Diffraction TomographyChapter Summary