---
title: Kinetic Super-Resolution in Dynamic Imaging
url: https://www.emergentmind.com/topics/kinetic-super-resolution
type: topic
---

# Kinetic Super-Resolution in Dynamic Imaging

Kinetic super-resolution designates a family of inverse problems in which resolution enhancement is coupled to motion, temporal evolution, or governing dynamics rather than treated as purely spatial interpolation. In the literature summarized here, the term covers recovery of moving point sources from low-frequency multitime measurements, motion-aware reconstruction from subpixel translations or motion blur, and coarse-to-fine reconstruction of turbulent or PDE states using physics-informed or generative models [1803.03251, 2505.15961, 1611.07767, 2502.10876, 2210.17319, 2410.20884, 2311.09740, 2507.03567, 2508.13770]. A related but distinct usage appears in interferometric phase estimation, where “super-resolution” refers to narrower phase fringes rather than spatiotemporal state recovery; that distinction is explicit in the phase-measurement literature [1705.05609].

## 1. Scope and defining characteristics

Across these works, kinetic super-resolution is characterized by the direct use of kinematic or dynamical structure inside the reconstruction map. Motion may appear as a latent variable to be estimated jointly with position, as a known measurement operator induced by sensor trajectory, as temporal consistency constraints between high-resolution frames, or as a PDE-constrained evolution law on the reconstructed state. This suggests that the adjective “kinetic” refers less to a single algorithm than to a class of formulations in which time dependence is part of the inverse problem itself.

A concise typology appears below.

| Setting | Core formulation | Representative papers |
|---|---|---|
| Moving sparse sources | Recover \((x,v,w)\) jointly in phase space from multitime low-frequency data | [1803.03251] |
| Motion-aware image/video SR | Use subpixel translations, motion blur, or inter-frame optical flow as reconstruction constraints | [2505.15961], [1611.07767], [2502.10876] |
| Dynamics-aware PDE/turbulence SR | Infer fine fields from coarse observations using PDE residuals, solver unrolling, or generative conditional models | [2210.17319], [2410.20884], [2311.09740], [2507.03567], [2508.13770] |

The same literature also emphasizes what kinetic super-resolution is not. It is not merely geometric upsampling, because the target is often a physically admissible high-resolution state rather than a visually sharper field. It is also not identical to conventional super-resolution defined from artificially downsampled ground truth. In CFD-oriented work, this distinction is central: coarse solver outputs are not equivalent to downsampled fine solutions, because downsampling preserves more of the underlying physics than an authentic coarse-grid simulation [2311.09740].

## 2. Phase-space lifting and dynamic spike super-resolution

A mathematically explicit formulation of kinetic super-resolution appears in the moving-spike setting. The unknown signal is a time-varying measure
$$
\mu_t = \sum_{i=1}^N w_i\,\delta_{x_i+v_i t}, \qquad t\in[-\delta,\delta],
$$
with discrete observation times
$$
t_k = kT,\qquad k=-K,\dots,K,
$$
and low-frequency measurements
$$
y_k = F\mu_{t_k}.
$$
In the low-frequency Fourier case,
$$
(F\mu_t)_\ell = \int_{[0,1]^d} e^{-2\pi i x\cdot \ell}\,d\mu_t(x), \qquad |\ell|\le f_c.
$$
The key idea is to lift the problem from physical space to phase space \((x,v)\), replacing framewise reconstruction plus tracking by a single sparse inverse problem on
$$
\mathcal{Q} = \{(x,v)\in\mathbb{R}^{2d}:\ x+kTv\in[0,1]^d,\ \forall k\in[-K,K]\},
$$
with sparse measure
$$
w=\sum_{i=1}^N w_i\,\delta_{(x_i,v_i)}.
$$
The reconstruction program becomes
$$
\min_{\nu\in\mathcal{M}(\mathcal{Q})} \|\nu\|_{\mathrm{TV}}
\quad\text{subject to}\quad G\nu = y. \tag{6}
$$
This replaces the standard two-step pipeline—static localization at each time, then tracking—by simultaneous recovery of positions, velocities, and weights [1803.03251].

The theory is organized around dynamical dual certificates of the form
$$
q(x,v)=\sum_{k=-K}^K\sum_{|\ell|\le f_c} c_{k,\ell}\,e^{2\pi i \ell\cdot(x+kTv)}. \tag{9}
$$
Exact recovery follows when \(q(x_i,v_i)=n_i\) and \(|q(x,v)|<1\) away from the support. A notable obstruction is the existence of ghost particles, defined as alternative trajectories that pass through the same observed locations at all sampled times. Theorem 5 states that if each frame is sufficiently separated so that a static dual certificate exists and the configuration has no ghost particles, then the true dynamic measure is the unique solution of the dynamical TV problem. Proposition 7 shows that, for particles drawn independently from absolutely continuous distributions, ghost particles are absent with probability 1. The noisy discrete setting is also treated: in 1D, under \(f_c\ge 128\) and a discrete analogue of the no-ghost condition, Theorem 10 yields
$$
\|\hat w - w\|_1 \le C\,(\mathrm{SRF}_x)^2\,\delta,
$$
with
$$
\mathrm{SRF}_x = \frac{1}{\Delta x\,f_c}.
$$

The numerical experiments use 1D low-frequency Fourier measurements with \(f_c=20\) and \(K=2\) (five frames). They report that dynamical reconstruction succeeds more often than static reconstruction, especially when particles are close and static recovery is difficult. The same framework is then applied to ultrafast ultrasound localization microscopy, where microbubbles act as moving point reflectors; the proposed workflow selects time intervals with approximately constant total \(\ell^2\)-norm and applies the dynamical reconstruction to each interval, yielding super-resolved vessel imaging together with blood-flow velocity estimation [1803.03251].

## 3. Structured motion as a measurement resource

A second branch of kinetic super-resolution treats motion itself as an informative measurement operator. In structured-motion imaging, the low-resolution data are acquired under known sensor motion, either as a grid of subpixel translations or as continuous motion during exposure. For a moving sensor trajectory \((x(t),y(t))\), the recorded low-resolution image satisfies
$$
I[i,j] = \int_0^1 \int_0^\Delta \int_0^\Delta g(x(t)+i\Delta+x,\; y(t)+j\Delta+y)\,dx\,dy\,dt.
$$
The moving-sensor image can be represented as a convolution with an occupancy map \(q\) induced by the motion. With super-resolution factor \(f\), interlacing all \(f^2\) subpixel shifts yields
$$
H = J \otimes B,
$$
so the problem reduces to deconvolution with a box filter. For moving-sensor data the model becomes
$$
I = (J \otimes Q \otimes B)\downarrow_f.
$$
The central claim is that motion blur is not necessarily a nuisance: with high-precision motion information, sparse image priors, and convex optimization, it can aid reconstruction. Numerical experiments show that pseudo-random motion can reconstruct a high-resolution target from a single low-resolution image, and the reported real-data system works up to super-resolution factor \(f=8\) [2505.15961].

In video super-resolution, motion can be built directly into the high-resolution reconstruction rather than used only for pairwise frame alignment. A variational model jointly reconstructs a batch of high-resolution frames \(u=(u^1,\ldots,u^n)\) from low-resolution frames \(f=(f^1,\ldots,f^n)\):
$$
\min_{u} \sum_{i=1}^{n} \|D(b *u^i)-f^i \|_1 + \alpha \inf_{u=w+z} \big(R_{\text{temp}}(w) + R_{\text{spat}}(z)\big).
$$
The temporal and spatial regularizers are coupled through the motion-compensated discrepancy
$$
W(u^i,u^{i+1})(x) = \frac{u^{i}(x) - u^{i+1}(x + \boldsymbol{v}^i(x))}{h},
$$
and the scale \(h\) is chosen automatically. The essential consequence is linear, rather than quadratic, growth in the number of motion-estimation problems, because only neighboring-frame flows are required. On the reported \(4\times\) benchmark, the proposed method attains 29.19 dB average PSNR, compared with 29.13 dB for VDSR, and it reaches 23.97 dB on the synthetic tube example [1611.07767].

Classical multi-frame formulations remain important in this lineage. A standard observation model writes each low-resolution frame as
$$
Y_k = D_k M_k B_k X + n_k,
$$
or equivalently \(Y_k = W_k X + n_k\), where \(M_k\) is the motion matrix, \(B_k\) the blur operator, and \(D_k\) the decimation matrix. One thesis-level treatment combines Horn–Schunck optical flow with TV-regularized reconstruction,
$$
L(X)=\sum_{k=1}^N \|H_k X - Y_k\|^2 + \lambda\, TV(X), \qquad H_k = D_k M_k B_k,
$$
and solves the resulting problem by majorization-minimization with conjugate-gradient inner steps. The method yields visually improved high-resolution outputs from multiple \(64\times64\) low-resolution frames, but the implementation is explicitly described as not real-time, with reported processing around 1.5 hours in MATLAB experiments [2502.10876].

## 4. Physics-informed reconstruction from sparse or coarse observations

In dynamical systems and turbulence, kinetic super-resolution is often defined by the requirement that reconstructed fine scales remain consistent with the governing PDE. One formulation learns a map
$$
f_\theta:\; u(\Omega_L,t)\mapsto u(\Omega_H,t),
$$
from sparse observations on a coarse grid \(\Omega_L\) to a high-resolution state on \(\Omega_H\), while minimizing a combined observation and residual loss:
$$
\mathcal{L}_\theta = \alpha \mathcal{L}_{\mathcal{O}} + \mathcal{L}_{\mathcal{R}}.
$$
The observation term enforces agreement with the known low-resolution samples, whereas the residual term penalizes violation of
$$
\partial_t \mathbf{u} - \mathcal{N}(\mathbf{u};\lambda)=0.
$$
For Kolmogorov flow, the governing system is the 2D incompressible Navier–Stokes equation
$$
\nabla \cdot \mathbf{u}=0,
$$
$$
\partial_t \mathbf{u} + \mathbf{u}\cdot\nabla \mathbf{u} = -\nabla p + \nu \Delta \mathbf{u} + g(\mathbf{x}),
$$
with
$$
\nu = \frac{1}{34}, \qquad g(\mathbf{x}) = [\sin(4x_2),0]^\top.
$$
The implementation uses a VDSR architecture with bi-cubic upsampling, residual learning, periodic padding, a differentiable pseudospectral discretization, forward Euler time integration with \(\Delta t = 0.005\), and residual evaluation in the Fourier domain. The reported setup uses \(\Omega_H\in\mathbb{R}^{150\times150}\), \(\Omega_L\in\mathbb{R}^{10\times10}\), 2048 training windows, 256 validation windows, \(\tau=2\), Adam, learning rate \(3\times10^{-4}\), and \(\alpha=10^5\). Average relative \(\ell^2\) errors are 0.0872 for the physics-informed CNN, 0.2091 for bi-linear interpolation, and 0.1717 for bi-cubic interpolation [2210.17319].

A distinct but related strategy puts dynamics directly into the training loss through a differentiable solver. For forced two-dimensional turbulence on the torus, the governing equation is
$$
\partial_t \omega + \mathbf u \cdot \boldsymbol{\nabla} \omega = \frac{1}{Re}\Delta \omega - \alpha \omega - n \cos(n y),
$$
with velocity recovered from
$$
\Delta \psi = -\omega, \qquad u = \partial_y \psi,\quad v = -\partial_x \psi.
$$
Instead of supervising the network by a high-resolution target at \(t=0\), the method trains an initial condition whose forward evolution remains consistent with observations. The coarse-only loss is
$$
\mathscr L_{TC} = \frac{1}{N_S N_T}\sum_{j=1}^{N_S}\sum_{k=1}^{N_t}\| \mathcal C \circ \boldsymbol\varphi_{t_k}(\mathbf u_j) - \mathcal C \circ \boldsymbol\varphi_{t_k}\circ \mathcal N_{\boldsymbol \Theta} \circ \mathcal C ( \mathbf u_j) \|^2.
$$
This requires a fully differentiable solver in the loop; the implementation uses JAX-CFD with a pseudospectral solver, high-resolution simulations on \(128\times128\) for \(Re=100\) and \(512\times512\) for \(Re=1000\), forcing mode \(n=4\), and \(\alpha=0.1\) at \(Re=1000\). The model has reconstruction errors similar to standard supervised super-resolution despite using no high-resolution reference data. At \(Re=1000\) and coarse factor \(M=32\), the coarse-only error at \(t=0\) is only about \(1.5\times\) the high-resolution-trained model’s error, and after forward evolution the errors become comparable to, and sometimes smaller than, those from high-resolution-trained models. The same work reports that the learned model outperforms variational data assimilation for initial state estimation, with assimilation errors exceeding 40% at the harsher \(M=32\) setting [2410.20884].

These formulations treat super-resolution as recovery of missing physics rather than filling in between pixels. In both cases, temporal evolution supplies the admissibility criterion for the fine-scale reconstruction.

## 5. Coarse-grid to fine-grid PDE prediction

Another major formulation redefines super-resolution specifically for PDE-based computation. The mapping is written as
$$
f: C \rightarrow F,
$$
where \(C\) is a coarse-grid CFD solution and \(F\) is a fine-grid CFD solution. The defining claim is that coarse inputs should come from actual coarse-grid simulations rather than from artificial downsampling of fine data. A physics-infused UNet is used to learn the nonlinear relationship between coarse mesh data
$$
C(x_1,x_2,\dots,x_n)\in \mathbb{R}^m
$$
and fine mesh data
$$
F(x_1,x_2,\dots,x_n)\in \mathbb{R}^d,\qquad m \ll d.
$$
The total loss is
$$
\mathcal{L}_{total} = \mathcal{L}_{data} + \mathcal{L}_{physics},
$$
where the physics term is an MSE over convective and diffusive or conductive derivative terms computed by 2nd-order finite differences [2311.09740].

The architecture consists of four contracting convolution blocks, a bottleneck, four expansion blocks with skip connections, a \(1\times1\) output projection, and a bilinear upsampler. The formulation is demonstrated on 2D Burgers’ equation, methane combustion, and an industrial heat exchanger. For Burgers’ equation,
$$
\partial_t \mathbf{u}+\mathbf{u}\cdot \nabla \mathbf{u}=\nu\left(\nabla^2 \mathbf{u}\right),
$$
the reported grid mapping is \(51\times51 \rightarrow 401\times401\). For methane combustion the mapping is \(50\times20 \rightarrow 500\times100\), and for the heat exchanger it is \(30\times30 \rightarrow 480\times480\). The Burgers case uses
$$
\mathcal{L}_{physics}=\alpha_{conv}\mathcal{L}_{conv}+\alpha_{diff}\mathcal{L}_{diff},
\qquad
\alpha_{conv}=10^{-4},\qquad \alpha_{diff}=5\times 10^{-10},
$$
while the methane and heat-exchanger cases use analogous weighted convective and diffusive or conductive penalties [2311.09740].

The reported numerical results are explicit. For Burgers’ \(U_x\), bilinear interpolation has RMSE 0.4927, bicubic 0.5133, UNet 0.0283, and PIUNet 0.0207; for \(U_y\), PIUNet reaches RMSE 0.0225 and \(R^2=0.9965\). The paper also reports mean \(L_2\) error of PIUNet on the fine mesh as 0.0097, compared with coarse-mesh \(L_2\) error 0.0446, corresponding to roughly 78% reduction. In methane combustion, the adiabatic flame temperature RMSE decreases from 78.317 for bilinear and 77.374 for bicubic to 20.954 for PIUNet, with about 73% improvement in \(T_{adia}\) RMSE versus baseline interpolation; methane mass fraction RMSE improves by about 20.47%. In the heat-exchanger case, PIUNet reduces RMSE approximately by 36.8% for \(T_1\), 30.7% for \(T_2,T_3\), 35.04% for \(T_{m_1}\), and 20.06% for \(T_{m_2},T_{m_3}\). The corresponding coarse-solver-plus-network workflow yields about 24× speedup for Burgers, about 93× for methane combustion, and about 235× for the heat exchanger [2311.09740].

This line of work frames kinetic super-resolution as a surrogate fine-mesh solver. The output is not a photographic refinement but an approximate fine-grid PDE state intended to preserve convective transport, diffusion or conduction, and boundary-condition structure.

## 6. Turbulence, subgrid physics, and generative reconstruction

In turbulence, kinetic super-resolution is often assessed by whether it reconstructs spectra, vorticity statistics, dissipation, and localized high-frequency structure. One particle-aware formulation uses a conditional GAN for two-way coupled particle-laden turbulent flows. The generator reconstructs the high-resolution velocity field as
$$
\hat{\bm{\xi}} = \{\hat u,\hat v,\hat w\} = G(\bm{\xi}_{LR}, k_{sgs}, \rho_{p,eff}),
$$
where
$$
k_{sgs} = 0.5\left(\overline{u_i u_i} - \overline{u_i}\,\overline{u_i}\right)
$$
is subgrid kinetic energy and
$$
\rho_{p,eff} = \frac{N_p m_p}{V_{cell}}
$$
is effective particle mass density. The discriminator is conditioned on the low-resolution input and stationary-wavelet-transform detail coefficients,
$$
D = D\left(\bm{\xi}_{HR}, \bm{\xi}_{LR}, SWT(\bm{\xi}_{HR}), SWT(\bm{\xi}_{LR})\right),
$$
with
$$
SWT(\bm{\xi}) = \{LH_2(\bm{\xi}), HL_2(\bm{\xi}), HH_2(\bm{\xi})\}.
$$
The generator uses a deep RRDB network with 16 RRDB blocks, each with 3 RDBs; the discriminator is a U-Net with spectral normalization. Training and testing use DNS from forced homogeneous isotropic turbulence and decaying turbulence, with \(\mathrm{St}_{\eta} = 0.6, 1, 6\), \(ML = 0.49, 0.75\), \(128\times128\) slices, 16,000 training samples per case, 800 test samples per case, total training set 320,000 samples, and total test set 16,000 samples. The ablation that masks particle conditioning at inference shows that large-scale structures remain somewhat reasonable but high-wavenumber energy decays too quickly and vorticity structures become less accurate. Including particle input improves \(NRMSE_{sp}(E_{1D,sgs})\) from \(6.8\times10^{-2}\) to \(4.6\times10^{-2}\), stress-tensor error from \(1.3\times10^{-1}\) to \(1.1\times10^{-1}\), and velocity error from \(4.8\times10^{-2}\) to \(4.4\times10^{-2}\). The same model reproduces both positive subgrid dissipation and backscatter, in contrast to the standard Smagorinsky closure
$$
\varepsilon_{sgs,SMAG} = 2 (C_s \Delta_s)^2 |\overline{S}|^3,
$$
which is always positive [2507.03567].

A generative alternative uses stochastic interpolants for two-dimensional turbulence. The target field \(u\) is sampled on a \(128\times128\) grid, while the coarse input is obtained by low-pass filtering with \(k_{\text{cutoff}}=8\) and downsampling to \(16\times16\). The conditional generative task is to learn
$$
\rho(x_1\mid x_0),
$$
with interpolant
$$
I_\tau = \alpha_\tau x_0 + \beta_\tau x_1 + \sigma_\tau W_\tau,
$$
and SDE
$$
dX_\tau = b_\theta(X_\tau,x_0,\tau)\,d\tau + \sigma_\tau\,dW_\tau, \qquad X_{\tau=0}=x_0.
$$
For the turbulence application,
$$
\alpha_\tau = 1-\tau,\qquad \beta_\tau = \tau^2,\qquad \sigma_\tau = 0.1(1-\tau).
$$
Inference uses a Heun SDE integrator with 100 pseudo-timesteps. The main methodological extension is a patch-wise strategy: the \(128\times128\) field is split into 16 patches of size \(32\times32\), reconstructed in two stages by a free-generator and a cond-generator arranged in a checkerboard-like procedure. This reduces patch-boundary artifacts, especially in the vorticity field. The reported diagnostics include the kinetic energy spectrum, vorticity,
$$
\omega=\nabla\times \mathbf{u},
$$
local dissipation
$$
\epsilon = 2\mu s_{ij}s_{ij},\qquad s_{ij}=\left(\frac{\partial u_i}{\partial x_j}+\frac{\partial u_j}{\partial x_i}\right),
$$
and PDFs of spatially averaged kinetic energy, vorticity skewness, and spatially averaged dissipation. For dissipation, the base field has KL divergence 0.5480, while SI reduces it to 0.0106 or 0.0040 depending on the variant; Wasserstein-1 drops from 2.5809 to 0.3013 or 0.1094. The paper states that stochastic interpolants outperform flow-matching and diffusion baselines across a range of metrics, and that the patch-wise approach often matches or exceeds full-field reconstruction [2508.13770].

Taken together, these results place subgrid physics at the center of kinetic super-resolution. The relevant target is not only pointwise accuracy but the recovery of unresolved energy, dissipation, intermittency, and particle-modulated fine-scale structure.

## 7. Conceptual distinctions, misconceptions, and recurring limitations

Several recurring distinctions structure the field. First, dynamic or kinetic super-resolution is repeatedly contrasted with static reconstruction followed by tracking. In the moving-spike literature, the standard pipeline fails when particles are too close in any single frame, ignores temporal information during reconstruction, and requires a separate tracking stage; phase-space lifting is proposed specifically to avoid those drawbacks [1803.03251].

Second, motion is not uniformly treated as degradation. Structured-motion imaging argues that motion blur can be helpful for super-resolution and that pseudo-random motion can encode usable spatial information in a single blurred exposure [2505.15961]. Video super-resolution similarly uses motion to couple unknown high-resolution frames directly, rather than merely transporting low-resolution evidence frame by frame [1611.07767]. At the same time, classical motion-based methods remain limited by motion-estimation accuracy, brightness constancy assumptions, occlusions, and computational cost [2502.10876, 1611.07767].

Third, several PDE-oriented works reject the assumption that super-resolution should be trained on downsampled fine solutions. The coarse CFD literature explicitly states that downsampling fine data is not the correct analogue for PDE super-resolution, because true coarse-grid simulations embody different physics [2311.09740]. A related misconception is that turbulent super-resolution is equivalent to image sharpening. Physics-informed and dynamics-in-the-loss methods instead define success through PDE consistency, spectral fidelity, or future evolution under a solver [2210.17319, 2410.20884].

Finally, the term “super-resolution” itself is field-dependent. In deterministic phase measurements, super-resolution means interference fringes narrower than the usual half-wavelength periodicity, while super-sensitivity means phase uncertainty below the shot-noise limit,
$$
\Delta \phi_{\rm SNL} \sim \frac{1}{\sqrt{N}}.
$$
That work uses a coherent state, a squeezed vacuum state, homodyne detection, and a dichotomic windowing strategy
$$
\Pi_0 = \int_{-a}^{a} dp\, \lvert p \rangle\langle p\rvert,\qquad \Pi_1 = I - \Pi_0,
$$
to achieve both effects simultaneously, reporting 430 photons, a 22-fold improvement in phase resolution, and a 1.7-fold improvement in sensitivity. The paper’s distinction between fringe width and estimation precision clarifies that “super-resolution” is not a universal synonym for better inference; its technical meaning depends on the measurement problem being posed [1705.05609].

In the kinetic literature proper, a plausible implication is that future progress will continue to depend on embedding the correct dynamics into the inverse problem. The most successful formulations here do not simply increase pixel density. They reconstruct trajectories, exploit structured motion, or infer fine-scale states constrained by PDEs, spectra, and multiscale statistics.

Source: https://www.emergentmind.com/topics/kinetic-super-resolution