Unverified 2026

Canonical Push-Forward SDF Encoder

Implementation & benchmark of arXiv:2608.02306 — The Push-Forward Transform for Continuous and Robust Comparison of Dynamic Shapes

Usefulness6/10
Difficulty6/10
Novelty6/10

Source paper: The Push-Forward Transform for Continuous and Robust Comparison of Dynamic Shapes arXiv:2608.02306 ⓘ · analyzed Aug 31, 2026

AI-generated research hypothesis, automatically tested. Not peer-reviewed.

Idea description

Warp each observed shape and its interior scalar fields into a fixed reference domain using a learned diffeomorphism, then process the resulting canonical SDF with a CNN or 3D encoder. The representation should be insensitive to translation, rotation, reflection, scale, and re-parameterization when the correspondence module is constrained to absorb those transformations rather than shape changes.

Formulas

$$\partial S=\{x\in\Omega:L_S(x)=0\},\qquad S=\{x\in\Omega:L_S(x)>0\},\qquad S^c=\{x\in\Omega:L_S(x)<0\}.$$
$$\Psi:\overline{S_r}\to\overline{S},\qquad f_{S^*}(u)=f_S(\Psi(u)),\quad u\in S_r.$$
$$L_{S^*}(u)=L_S(\Psi(u)),\qquad \mathcal{L}_{\mathrm{PF}}=\frac{1}{|S_r|}\sum_{u\in S_r}\left|L_{S_1}(\Psi_1(u))-L_{S_2}(\Psi_2(u))\right|.$$
$$\mathcal{L}_{\mathrm{diff}}=\frac{1}{|S_r|}\sum_{u\in S_r}\left[\operatorname{softplus}\!\left(m-\det J_{\Psi}(u)\right)+\operatorname{softplus}\!\left(m-\det J_{\Psi^{-1}}(\Psi(u))\right)\right],\quad J_{\Psi}(u)=\frac{\partial\Psi(u)}{\partial u}.$$

Mathematical statement

Definition 1 models a shape boundary as the zero level set of a $C^1$ function $L_S:\Omega\to\mathbb{R}$ with nonzero gradient on the boundary: $\partial S=\{x:L_S(x)=0\}$, with interior $S=\{x:L_S(x)>0\}$ and exterior $S^c=\{x:L_S(x)<0\}$. Definition 2 introduces a diffeomorphism $\Psi:\overline{S_r}\to\overline{S}$ from a reference shape domain $S_r$ to the observed shape domain $S$, and transports a scalar field $f_S:S\to\mathbb{R}$ to reference coordinates. The extracted display is truncated, so the implementable scalar-field pullback is $f_{S^*}(u)=f_S(\Psi(u))$ for $u\in S_r$; for an SDF this becomes $L_{S^*}(u)=L_S(\Psi(u))$. Here $u$ is a reference-domain coordinate, $x=\Psi(u)$ is the corresponding observed-domain coordinate, and $J_\Psi(u)=\partial\Psi(u)/\partial u$ is the deformation Jacobian. The diffeomorphism assumption requires one-to-one, onto, differentiable mapping with nonsingular Jacobian, preventing folds and collapsed regions. The neural adaptation compares transported fields pointwise on $S_r$ and penalizes invalid Jacobians.

Implementation notes

1. Integration point: add a canonicalization front-end before a 2D CNN, 3D sparse CNN, or point/voxel encoder. The input is an SDF sampled on a fixed reference grid $u\in S_r$; optional intensity or molecular channels are warped with the same map. Predict a displacement field $v_\theta(u)$ and define $\Psi_\theta(u)=u+v_\theta(u)$, followed by differentiable grid sampling of the observed SDF at $\Psi_\theta(u)$. For rigid nuisance removal, predict translation, scale, and rotation first, and compose that transform with the nonrigid map. 2. Pseudocode: sample reference grid $u$; predict $v=\mathrm{DeformNet}(L_S)$; set $x=u+v$; compute $L^*=\mathrm{grid\_sample}(L_S,x)$; encode $z=\mathrm{Encoder}(L^*)$; optimize task loss plus $\lambda_{PF}\|L^*_1-L^*_2\|_1$ for augmented views of the same shape, Jacobian fold penalty $\mathcal L_{diff}$, and optionally inverse consistency $\|\Psi^{-1}(\Psi(u))-u\|^2$. 3. The mathematical quantities computed from the paper are the level-set zero crossing, transported field, Jacobian determinant, and reference-domain comparison. The deformation is estimated by the neural network; estimate Jacobians with autodiff or centered finite differences. Use a smooth SDF target and clamp distances to $[-d_{max},d_{max}]$. 4. First experiment: train a small 2D shape classifier on rotated, reflected, translated, and rescaled MPEG-7 or synthetic silhouettes. Compare a raw-image CNN, a standard spatial transformer, and the PF-T SDF encoder at equal parameter count. Measure accuracy under held-out transformations, embedding cosine consistency between augmented views, and deformation fold rate. A successful result is substantially smaller accuracy degradation under unseen transformations and lower within-class embedding variance without reducing clean accuracy. Then repeat on a small ShapeNet subset with a 3D sparse encoder.

Verification

This idea has not been verified yet.

Verification happens in two stages: Stage 1 — a mechanism check on a toy system confirms the claimed mathematical phenomenon reproduces; Stage 2 — a benchmark implements the idea on a real (small) neural network task and compares it against a tuned baseline over 8 paired seeds with a permutation test.

Artifacts

Artifacts unavailable.