---
title: Fourier Magnetic Imaging
url: https://www.emergentmind.com/topics/fourier-magnetic-imaging
type: topic
---

# Fourier Magnetic Imaging

Fourier magnetic imaging denotes a class of measurement and reconstruction strategies in which magnetic, spin, or current information is encoded, sampled, or inverted in Fourier space rather than obtained solely by direct real-space rastering. In the nitrogen-vacancy (NV) implementation introduced by Arai et al., pulsed magnetic field gradients phase-encode spatial information on NV electronic spins in wavenumber or k-space followed by a fast Fourier transform to yield real-space images with nanoscale resolution, wide field-of-view (FOV), and compressed sensing speed-up [1409.2749]. Closely related Fourier-domain constructions appear in spread-spectrum MRI, nanoscale nuclear-spin MRI, magnetic current imaging from stray-field maps, and lensless magnetic imaging schemes based on holography, ghost imaging, and Fourier-plane optical processing [1110.5870], [1302.2977], [2507.11853], [1801.10046], [2212.10183], [2212.02975].

## 1. Fourier-domain encoding and image formation

At its most basic, Fourier magnetic imaging follows the standard magnetic resonance relation between an object and its reciprocal-space signal. In Fourier MRI, the measured signal at k-space location $(k_x,k_y)$ is
$$
S(k_x,k_y)=\iint_{\mathbb R^2}\rho(x,y)\,e^{-2\pi i\,(x\,k_x+y\,k_y)}\,dx\,dy
=\bigl(\mathcal F\{\rho\}\bigr)(k_x,k_y),
$$
with $\rho(x,y)$ the transverse magnetization distribution in the slice [2506.18638]. In nanoscale pulsed MRI, the same structure appears as
$$
S(k)=\int \rho(r)\,e^{-i\,k\cdot r}\,dr,\qquad
k(t)=\gamma\int_0^t G(t')\,dt',
$$
where the encoding wave-vector is accumulated by the applied time-varying gradient $G(t)$ [1302.2977].

For NV-center Fourier magnetic imaging, the phase acquired by a spin at position $\mathbf r$ is
$$
\phi(\mathbf r)=\gamma\int_0^T \mathbf G(t)\cdot\mathbf r\,dt \approx 2\pi\,\mathbf k\cdot\mathbf r,
$$
and the measured fluorescence at each sampled point in k-space is
$$
s(\mathbf k)\propto \cos\bigl(2\pi\,\mathbf k\cdot\mathbf r+\theta\bigr),
$$
with $\theta$ carrying any additional phase, for example from a local AC field [1409.2749]. This formal analogy with MRI is exact at the level of reciprocal-space encoding.

The standard reciprocal-space relations between sampling density, FOV, and spatial resolution also carry over. For NV Fourier magnetic imaging,
$$
\mathrm{FOV}\approx \frac{2\pi}{\Delta k}, \qquad \Delta r\approx \frac{\pi}{k_{\max}},
$$
so the k-space step $\Delta k$ sets the FOV and the maximum extent $k_{\max}$ sets the nominal spatial resolution [1409.2749]. In Fourier MRI formulated distributionally, rectangular sampling on a grid with interval $\Delta k$ generates periodic replications of the image with period $\mathrm{FOV}=1/\Delta k$ in each direction; to avoid aliasing one chooses $\Delta k\le 1/\mathrm{FOV}$ [2506.18638]. Truncating k-space to $|k|\le k_{\max}$ corresponds to multiplying the signal by a rectangular window, hence to a sinc point-spread function in image space with width $\Delta x\approx 1/(2k_{\max})$; tapered or Gaussian windowing reduces ringing at the expense of resolution [2506.18638].

These relations clarify that Fourier magnetic imaging is not defined by a single sensor platform. What unifies the methods is the use of reciprocal-space sampling, reciprocal-space inference, or Fourier-plane optical processing as the primary route from measurement to magnetic structure.

## 2. NV-center implementations

The original NV implementation uses optically detected magnetic resonance with pulsed field gradients generated by on-chip microcoils. A 3–5 $\mu$s green laser pulse at $\lambda=532$ nm initializes the NV into $|0\rangle$, a dynamical-decoupling microwave sequence such as spin-echo or $2n$-pulse CPMG/XY protects coherence during a total free-precession duration $\tau$, pulsed magnetic field gradients are applied during each free-precession interval, and a final $\pi/2$ pulse maps the accumulated phase into fluorescence read out in the 640–800 nm band [1409.2749]. The gradient coils are anti-Helmholtz pairs separated by $\sim100$ $\mu$m; a 1 A pulse generates a $\sim0.7$ G/$\mu$m gradient at the center, uniform within 1% over $15\,\mu\mathrm m\times15\,\mu\mathrm m$ [1409.2749].

Uniform k-space grids were demonstrated in one and two dimensions. In a one-dimensional example, $k\in[0,\,0.144\,\mathrm{nm}^{-1}]$ was sampled in 58 steps of $\Delta k\approx2.5\times10^{-3}\,\mathrm{nm}^{-1}$, giving FOV $\approx2.5$ $\mu$m and $\Delta r\approx11$ nm from the standard Fourier-imaging relations [1409.2749]. In two dimensions, grids such as $21\times21$ or $41\times41$ points were acquired. After tapered-cosine windowing with 10% taper, a discrete Fourier transform yields a complex real-space image whose magnitude reveals NV positions and whose phase encodes local AC field information [1409.2749].

Representative performance metrics establish the scale of the method. Arai et al. reported 3.5 nm spatial resolution for 1D single-NV imaging with SNR $\approx13$, 30 nm pixel resolution for 2D multiple-NV imaging resolving a 121 nm separation, field-of-view up to $15\,\mu\mathrm m\times15\,\mu\mathrm m$ in hybrid real + k imaging with dynamic range $\mathrm{FOV}/\Delta r\sim500$, magnetic sensitivity of $\sim1$ $\mu$T Hz$^{-1/2}$ per NV, and gradient sensitivity of $\sim14$ nT/nm Hz$^{-1/2}$ [1409.2749]. The same work implemented compressed sensing by randomly sampling only $M$ out of $N=2048$ k-points, using $\ell_1$ minimization with a partial DFT matrix. For $M=128$, the acceleration factor was $R=16$, reconstructed NV separations agreed within error bars with full-sample values, and phase-difference field estimates of $(5.0\pm0.5)\times10^{3}$ nT were obtained versus $(6.3\pm0.8)\times10^{3}$ nT for full sampling [1409.2749].

Later work pushed the approach to the sub-nanometer regime. A compact ambient platform with thermal drift compensation generated a pulsed magnetic field gradient of up to 13.5 G/$\mu$m and achieved localization of a single NV center with a spatial resolution of $0.28\pm0.10$ nm and a magnetic field measurement deviation of 9 nT [2603.22718]. In that implementation, a spin-echo sequence with longest gradient-pulse evolution $2\tau=500$ $\mu$s and maximum current $I_{\max}=10$ mA produced $G_{\max}=(3.26\pm0.06)$ G/$\mu$m at the NV location, corresponding to $k_{\max}\approx2.2834$ nm$^{-1}$ and a theoretical $\Delta x\approx1.38$ nm, while the observed full-width at half-maximum was $0.28\pm0.10$ nm [2603.22718]. The same platform reported $T_2=1.2$ ms and magnetic sensitivity $\eta_B\approx0.7$ $\mu$T/$\sqrt{\mathrm{Hz}}$ [2603.22718].

## 3. MRI, spread spectrum, and optimal Fourier sampling

Fourier magnetic imaging inherits much of its reconstruction theory from MRI. In standard under-sampled MRI, with $x\in\mathbb C^N$ the vectorized image and $\alpha=\Psi^*x$ an $S$-sparse expansion in an orthonormal sparsity basis $\Psi$, the measurement model is
$$
y = F_\Omega x + \text{noise} = F_\Omega \Psi \alpha + \text{noise},
$$
where $F_\Omega$ is the row-restricted Fourier operator [1110.5870]. In the spread-spectrum variant, a diagonal unit-modulus modulation $M$ or $C$ is applied before Fourier sampling:
$$
y = F_\Omega M x = F_\Omega M\Psi\alpha =: A_\Omega^c \alpha.
$$
Digital random pre-modulation uses Rademacher or Steinhaus sequences; the analog version uses a quadratic-phase chirp $c(\tau)=\exp(i\pi w|\tau|^2)$ [1110.5870].

The central effect of pre-modulation is coherence reduction. Without modulation, the mutual coherence between Fourier sensing rows and sparsity atoms can be high, forcing $m\sim N$. With random modulation, Puy et al. show that, with probability at least $1-\varepsilon$,
$$
\mu(F\,C,\Psi)\le \beta(F,\Psi)\sqrt{2\log(2N^2/\varepsilon)},
$$
where the modulus-coherence is
$$
\beta(F,\Psi)=\max_{i,j}\sqrt{\sum_{k=1}^N |F_{k,i}^*\,\Psi_{k,j}|^2}.
$$
If $F$ is the unit-magnitude Fourier basis, then for any $\Psi$ one has $\beta(F,\Psi)=1/\sqrt N$, so the modulated coherence is essentially $O(\sqrt{\log N}/\sqrt N)$ [1110.5870].

This leads to the universality claim. Standard compressed sensing gives recovery of every $S$-sparse $\alpha$ by $\ell_1$ minimization provided
$$
m \ge C\,N\,\mu^2\,S\,\log^4(N).
$$
After random pre-modulation this becomes
$$
m \ge \text{const}\cdot S\,\log^5(N),
$$
and for a universal sensing basis such as Fourier or Hadamard the number of measurements scales as
$$
m \ge C_\rho\,S\,\log^5(N),
$$
independent of the choice of sparsity basis $\Psi$ [1110.5870]. Numerical phase-transition experiments at $N=1024$ showed that, with random pre-modulation, Dirac, Haar, and Fourier sparsity cases collapse onto the same Donoho–Tanner optimal curve [1110.5870]. In numerical MRI examples using $m\approx0.4N$ random phase-encoded lines, chirp modulation with $\bar w\approx0.1$ improved reconstructions from PSNR $\approx20$ dB and SSIM $\approx0.65$ without chirp to PSNR $\approx30$ dB and SSIM $\approx0.90$ with chirp, while per-iteration complexity remained $O(N_w\log N_w)$ [1110.5870].

Two later developments address reconstruction beyond standard $\ell_1$ recovery. First, uncertainty quantification for Fourier MRI can be sharpened by reweighting without-replacement sampling so that the virtual Gram matrix satisfies $\mathbb E[(\tilde A)^*\tilde A/n]=I$. In a Shepp–Logan phantom experiment with $x^0\in\mathbb C^{32768}$, $m=0.6N$, and complex Gaussian noise with $\mathrm{SNR}\approx4.5\%$, reweighted debiasing at $\lambda=15\lambda_0$ reduced $\|\hat x^u-x^0\|_\infty$ from $\approx0.0616$ to $\approx0.0258$, reduced $\|R\|_\infty$ from $\approx0.0303$ to $\approx0.0117$, and achieved confidence-interval coverage of $97.85\%$ overall and $97.77\%$ on support [2407.13575]. Second, random anisotropic sampling combined with dualizable shearlet frames yields asymptotic optimality for cartoon-like functions, with
$$
\|f-R(f,\Delta_J)\|_2 \lesssim |\Delta_J|^{-1 + C\rho},
$$
and numerical experiments on $512\times512$ images reported 5–7 dB PSNR gain over wavelet-based schemes [1510.05029].

## 4. Fourier inversion of magnetic field maps into current density

A distinct branch of Fourier magnetic imaging reconstructs current density from measured magnetic fields. For a two-dimensional current density $\mathbf J(x',y')$ in the plane $z=0$, the out-of-plane field at height $z$ is given by the Biot–Savart integral
$$
B_z(x,y,z)=\frac{\mu_0}{4\pi}\iint
\frac{J_x(x',y')(y-y')-J_y(x',y')(x-x')}
{\bigl[(x-x')^2+(y-y')^2+z^2\bigr]^{3/2}}
\,dx'\,dy'
$$
in the SPIM formulation [2507.11853]. Its two-dimensional Fourier transform satisfies
$$
\widetilde B_z(k_x,k_y;z)
=
\frac{\mu_0}{2}\,e^{-kz}\,
\bigl[k_y\widetilde J_x(k_x,k_y)-k_x\widetilde J_y(k_x,k_y)\bigr],
\qquad
k=\sqrt{k_x^2+k_y^2},
$$
together with the divergence-free constraint
$$
-ik_x\widetilde J_x-ik_y\widetilde J_y=0
$$
[2507.11853]. The current-to-field operator is therefore diagonal in Fourier space, and the real-space convolution is replaced by multiplication by an inversion kernel.

The SPIM framework adds spatial pre-processing before the FFT-based inversion. Stage I rotates lock-in detector channels $I$ and $Q$ by a global angle $\varphi$ chosen to maximize the total variation of the new in-phase image $I'$. In the reported 3D-spiral SQUID scan, $\varphi\approx5.8^\circ$, which sharpens $I'$ by 0.3% and suppresses $Q'$ by 25% [2507.11853]. Stage II applies affine alignment with small rotation $\theta$ or skew $\alpha$; values of $\theta\approx0.3^\circ$ and $\alpha\approx0.3^\circ$ eliminated visible misalignment [2507.11853]. Stage III converts $I'(x,y)$ into magnetic field using the SQUID calibration constant $3.632$ $\mu$T/V, computes $\widetilde B_z$ via FFT, forms $\widetilde J_x,\widetilde J_y$ with a hard cut-off $k\le k_w$, and uses inverse FFT to recover $J_x(x,y)$ and $J_y(x,y)$. The best results were found for $k_w\approx 3/z$ with $z\approx120$ $\mu$m [2507.11853].

The classical Fourier route is effective but intrinsically ill-conditioned because the kernel decays exponentially at large $\|k\|$. In the wavelet-based analysis of two-dimensional magnetic current imaging, the kernel is
$$
g(k;z)=\frac{\mu_0 d}{2}\exp\bigl(-2\pi z\|k\|\bigr),
$$
so naïve deconvolution amplifies noise strongly at high spatial frequencies [2408.16550]. Fourier methods therefore use low-pass filters such as cosine taper, Gaussian cut-off, or sharp cut-off, but this forces a trade-off between noise and resolution. The reported consequences are explicit: at $z=5\,\mu\mathrm m$ and $\sigma=1.25\,\mu\mathrm T$, Fourier-filtered error can exceed 30%, whereas an $L1$-curl, divergence-free wavelet method keeps error below 15%; more generally, across all tested noise levels and standoffs, the new method reduces relative $L^2$ error by roughly a factor of two compared to the Fourier approach [2408.16550]. The Fourier formulation remains the baseline against which such regularized inversions are measured.

## 5. Lensless and Fourier-plane magnetic imaging

Related methods recover magnetic structure from Fourier-plane intensity, holographic interference, or higher-order correlations rather than from gradient phase encoding. In thermal-neutron Fourier-transform ghost imaging, a spatially incoherent polarized neutron beam is split into a sample arm with a bucket detector and a reference arm with a position-sensitive detector. The key observable is the fourth-order correlation function
$$
G^{(2,2)}(\xi_r,\xi_t)=\langle \psi^\dagger(\xi_r)\psi^\dagger(\xi_t)\psi(\xi_t)\psi(\xi_r)\rangle,
$$
and, using fermionic Wick factorization and Pauli anti-symmetry, the covariance of intensity fluctuations becomes negative:
$$
\langle\Delta I_r(\xi_r)\Delta I_t(\xi_t)\rangle
=
-\,v_i v_f
\left|
\int I_0\,h_r^*(\xi_r,\eta)\,h_t(\xi_t,\eta)\,d\eta
\right|^2.
$$
With the choice $d_r=d_1+d_2$, this yields a coincidence signal directly proportional to $|\mathcal F[\mathcal P S^p](q)|^2$, the modulus-squared of the lateral Fourier transform of the projected scattering potential, including both atomic and magnetic terms [1801.10046]. The method is lensless, and the stated resolution limit is the neutron de Broglie wavelength [1801.10046].

Fourier-transform holography (FTH) provides a different lensless route. The measured hologram is
$$
H(q)=|O(q)+R(q)|^2=|O|^2+|R|^2+O(q)R^*(q)+O^*(q)R(q),
$$
so a known reference isolates a phase-locked cross term and overcomes the phase problem in one single step of calculation [2212.10183]. In magnetic tomography, the phase contrast obeys
$$
\Delta\Phi(x,y)=\Phi_+(x,y)-\Phi_-(x,y)
=
-\,r_e\,\lambda\,\rho\,f_m'\int (\hat m\cdot\hat k)\,dz,
$$
making the reconstructed phase directly proportional to the projection of magnetization along the beam [2212.10183]. This approach produced a 3D full-vectorial image of a 800 nm-thick extended Fe/Gd multilayer in a 5 $\mu$m-diameter circular field of view with a resolution of approximately 80 nm [2212.10183].

Fourier-plane optical processing can also be used to disambiguate magnetic spectra. In NV-ensemble vector magnetometry, Fourier-plane amplitude masks and a linear polarizer generate four measurements $\boldsymbol\chi = W\mathbf c$, where $W$ is a $4\times4$ circulant weight matrix; inversion gives the isolated ODMR lineshapes for each NV orientation even when the eight resonances overlap [1705.09241]. This enables vector magnetic imaging at arbitrarily low fields and extends the field-dynamic range by a factor $\gtrsim3$ without increasing $B_{\rm bias}$ [1705.09241]. The measured crosstalk is $\sim1\%$, while the relative SNR is approximately 0.47 for NA 1.49/oil and approximately 0.25 for NA 0.75/air [1705.09241].

A camera-based Fourier-space MOKE setup uses a high-NA objective to map a broad angular spectrum of incident and reflected wave-vectors onto detector pixels and fits the resulting intensity maps to first-order Kerr expressions. With left and right circular input polarization and no analyzing optics, the method retrieves the three magnetization components together with the optical and magneto-optical constants [2212.02975]. In simulations with realistic camera noise, a single shot gives magnetization orientation error $\approx15^\circ$, ten averages give $\approx2.5^\circ$, and 250 averages give $\lesssim0.5^\circ$; errors in optical and Voigt constants fall correspondingly from the 20–30% range to $\lesssim0.1$–2% depending on parameter and averaging [2212.02975].

## 6. Performance envelope, limitations, and evolving inverse models

Across the surveyed literature, Fourier magnetic imaging spans several experimentally distinct regimes. The following representative values are explicitly reported.

| Implementation | Fourier quantity | Representative performance |
|---|---|---|
| NV Fourier magnetic imaging | Pulsed-gradient k-space of NV electronic spins | 3.5 nm in 1D; 30 nm in 2D; FOV up to $15\,\mu\mathrm m\times15\,\mu\mathrm m$ [1409.2749] |
| Sub-nanometer NV localization | Spin-echo k-space encoding with pulsed gradients | $0.28\pm0.10$ nm spatial resolution; 9 nT magnetic field deviation [2603.22718] |
| Nanoscale Fourier-transform MRI | Pulsed nuclear-spin Fourier encoding | Two-dimensional projection with approximately 10-nm resolution [1302.2977] |
| 3D magnetic FTH tomography | Reciprocal-space hologram and inverse FT | 800 nm-thick Fe/Gd, 5 $\mu$m FOV, approximately 80 nm resolution [2212.10183] |

Several limitations recur. In gradient-based Fourier imaging, achievable $k_{\max}$ is bounded by gradient strength, coherence time, and thermal or timing stability; in the NV case, stronger gradients and longer $T_2$ are identified as the route below 1 nm, while in nanoscale nuclear-spin MRI spin relaxation during long encoding reduces high-$k$ contrast [1409.2749], [1302.2977]. In Fourier inversion of magnetic field maps, the forward operator is a low-pass filter, so any formal inverse must control noise amplification with cut-offs or regularization, producing the familiar blur-versus-noise trade-off [2408.16550], [2507.11853]. In Fourier-plane optical decomposition, mask throughput and half-pupil blocking reduce sensitivity and introduce anisotropic PSF elongation of about $1.4\times$ along the blocked axis [1705.09241]. In lensless reciprocal-space methods, phase retrieval is a central issue unless a reference architecture such as FTH is used [2212.10183].

Recent work indicates a shift from purely kinematic Fourier inversion toward physics-informed Fourier models. In a Fourier-space approach to magnetization reconstruction from NV stray-field measurements, the forward model uses FFT-based stray-field calculations and Fourier-space upward continuation inside a variational functional
$$
J[m,h]=L_{\mathrm{data}}[m,h]+\lambda\,E_{\mathrm{total}}[m],
$$
where $E_{\mathrm{total}}$ includes exchange, demagnetizing, anisotropy, and Zeeman energies [2602.17180]. On synthetic data with Néel-type walls at $h_{\mathrm ref}=80$ nm and 3% Gaussian noise, the recovered height converges within 2 nm of 80 nm over a wide $\lambda$-plateau; for $\alpha\le5\%$ noise the height error remains below 5 nm [2602.17180]. On $Fe_{3-x}GaTe_2$ data with FOV $800\times800$ nm and $100\times100$ pixels, the optimal reconstruction gives $h^*\approx75$–80 nm, reproduces the measured stray-field map to within $<5\%$ RMS error, typically converges in 500–1 000 iterations, and yields final stray-field RMS residual $\sim0.02$ MA/m [2602.17180]. This suggests a convergence between explicit Fourier encoders, FFT-based forward operators, and micromagnetically constrained inverse problems.

The resulting picture is that Fourier magnetic imaging is not a single apparatus but a reciprocal-space paradigm. In one branch, pulsed gradients encode position directly into spin phase and the image is obtained by inverse Fourier transformation. In another, the magnetic observable is reconstructed from Fourier-domain field kernels, holograms, or Fourier-plane optical signatures. In both branches, the decisive quantities are the accessible bandwidth in reciprocal space, the conditioning of the inverse map, and the extent to which physical priors can be incorporated without erasing high-spatial-frequency magnetic information.

Source: https://www.emergentmind.com/topics/fourier-magnetic-imaging