REVIEW 3 major objections 5 minor 34 references
Landmark shape spaces with induced metrics
T0 review · 3 major / 5 minor · reviewed 2026-07-31 · grok-4.5
Pith's one-line read A screened elasticity operator unifies collision-free landmark metrics with rigid-motion shape invariance.
desk verdict Real operator-level unification of Kendall quotients and LDDMM landmarks; completeness is solid for free rigid-quotient flow, thinner once scale is constrained. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The screened elasticity operator L_s = −(Id−σ²Δ)^{s−1} ∘ div ∘ sym∇ (and its Green kernel k_s): it sits between ordinary Sobolev operators and the classical elasticity/Killing operator, keeps null-space exactly se(d), and yields a conditionally positive definite cometric on the annihilator of rigid-motion orbits that descends to the quotient.
What would settle it
Match configurations that differ only by local rigid motions inside disjoint regions: the new kernel must transport those interiors rigidly while standard identity-containing kernels do not; if landmarks still collide for s ≥ 1 + d/2, or local rigid motion is still penalized on the quotient, the unification fails.
Extended reading notes
Core claim
The screened elasticity operator has null-space exactly the rigid motions, is equivariant under rigid motions, and for sufficient order induces a stratified Riemannian metric on landmark configurations modulo rigid motions (and on a fixed-scale submanifold) that prevents collisions, is independent of landmark count, preserves local rigid motions, and removes global ones—unifying the two classical landmark geometries.
Load-bearing premise
That smoothing the classical elasticity operator while keeping only rigid motions in its null space is enough for complete, usable geodesics once the order is high enough, even though scale is not built into the operator and must be forced by a separate constraint.
Editorial extensions
If this is right
- Landmark matching can use one metric that both removes global rigid motion and preserves local rigidity, without a separate preregistration step.
- Metrics on landmark shapes become compatible across different numbers of landmarks because they descend from an ambient diffeomorphism metric.
- With enough screening, geodesic completeness and non-collision hold on the rigid-motion quotient for the stated dimension and order range.
- Scale-normalized shape comparison becomes a constrained Hamiltonian system on the unit-radius submanifold of the quotient.
- The same operator idea extends, in principle, from landmarks to surfaces and images so rigid and non-rigid registration sit in one metric framework.
Reading between the lines
- Flow-based stochastic shape models that previously ignored rigid invariance could adopt this operator so evolutionary or growth trajectories quotient global pose automatically.
- Because the kernel does not vanish at infinity in 2D, numerical matching on large domains may need careful truncation or boundary handling not fully settled by the present experiments.
- The stratified orbit-type structure implies matching algorithms should monitor rank drops of landmark spans; crossing singular strata could silently change the effective geometry.
- If the same null-space design can be ported to discrete mesh or image operators, unified rigid-plus-deformable registration might replace two-stage pipelines in imaging practice.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper constructs landmark shape spaces that combine Kendall-type invariances (quotient by SE(d) and fixed scale) with LDDMM-type regularity by inducing metrics from a screened elasticity operator L_s = −(Id−σ²Δ)^{s−1} ∘ div ∘ sym∇. This operator is SE(d)-equivariant with null-space exactly se(d); its Green’s kernel k_s is conditionally positive definite w.r.t. se(d) and, for s ≥ 1 + d/2, sufficiently regular that free Hamiltonian trajectories with initial momentum in the annihilator exist for all time without landmark collision. The induced stratified cometric on Land_n/SE(d) preserves local rigid motions while removing global ones; scale is imposed by restricting to the holonomic slice R̄ = 1 and integrating a constrained Hamiltonian system (RATTLE). Matching and geodesic evolution are demonstrated numerically with k_2.
Significance. If the construction is as claimed, it fills a genuine gap between Kendall shape spaces (invariances, but Euclidean and n-dependent, no collision control) and LDDMM landmark metrics (regular, ambient, collision-free, but typically not rigid-motion free and penalizing local rigid motion). The screened elasticity operator is a clean, well-motivated choice sitting between Sobolev operators (too small null-space) and pure elasticity/Beppo-Levi (too singular or too large null-space). The derivations of null-space, equivariance, CPD structure, Noether conservation of rigid-momentum maps, and horizontal lifts are standard geometric-mechanics arguments carried out carefully. Explicit kernels for s = 2, a constrained integrator, and illustrative matching experiments (including local-rigid preservation vs Gaussian/Matérn) make the proposal usable. The main applied payoff—unified rigid-plus-deformable registration and n-independent metrics with Kendall invariances—is real if the fixed-scale geometry is placed on the same footing as the free case.
major comments (3)
- [§4.5, §5, App. A.11, Abstract] Abstract and §1 claim the new spaces achieve both collision prevention and fixed scale as defining properties of the unification. Proposition §4.5 proves all-time existence and no-collision only for the free Hamiltonian on Land_n with α(0) ∈ se(d)(q)^ann and s ≥ 1 + d/2 (via C¹ regularity of k_s). Scale is not L_s-invariant (§5); it is imposed by the holonomic constraint R̄ = 1 whose momentum map does not Poisson-commute with H (nonzero bracket computed in App. A.11). Geodesics on the actual Kendall-type object Land_n^1 are therefore only a constrained system with Lagrange multiplier. Global existence, constraint regularity, and non-collision for that constrained flow are not proved—only locally set up and illustrated with RATTLE (§6–7). Either extend the completeness argument to the constrained system, or qualify the abstract/intro claims so that collision-free completeness is stated on
- [§3.2, §4.5, App. A.10] Section 4.5 states that free Hamiltonian trajectories project to geodesics on Land_n “as long as they do not cross more singular strata,” and the orbit-space stratification (§3.2) is otherwise left open for geodesic behavior. For the unification with Kendall (where singular strata and orbit types are part of the classical picture), the paper should either (i) cite/adapt known convexity or non-crossing results for proper isometric actions under the extended positive-definite cometric of Corollary A.10 and verify the hypotheses, or (ii) explicitly restrict the main theorems and numerical claims to the regular stratum S_k (open dense) and state that crossing into lower-rank configurations is not controlled. As written, the stratified Riemannian structure is defined, but geodesic completeness of the quotient including strata is not established.
- [§6–7, Figures 2–8] The numerical section supports local-rigid preservation and constraint stability but does not stress-test the load-bearing geometric claims. Experiments use only k_2, modest n, and qualitative grid/trajectory plots; there is no systematic check of long-time constraint drift vs step size, no comparison of geodesic distance or matching residual against a rigid-pre-registered Sobolev/LDDMM baseline on the same pairs, and no probe of near-collision or near-singular (low-rank) configurations where the completeness gap would matter. Adding at least one quantitative baseline comparison and a near-degeneracy experiment would substantially strengthen the applied half of the contribution without changing the theory.
minor comments (5)
- [§4.2, App. A.5] Equation (4)–(5) and App. A.5 give k_2 explicitly; for s > 2 the text defers to root perturbation / repeated-root partial fractions. A short statement of the leading regularity (C^k with k in terms of s,d) already appears, but a pointer to which formula is used in code for general s would help reproducibility.
- [Figure 2] Figure 2 caption compares k_2 to Gaussian and Matérn-3/2 with “rather small” σ; reporting the actual σ (and μ, λ) used in each panel would make the local-rigid claim easier to reproduce.
- [§3–5] Notation switches between Land_n, Land_n bar, Land_n^1, and g_{L_s} / g_1; a small notation table early in §3–4 would reduce friction.
- [throughout] Typos/style: “Mat´ ern” spacing, “Lam´ e”, “att=0”, “a to fixed-scale sphere” (§5), and occasional missing spaces before citations. Standard copy-edit pass.
- [§1.1, §2.2] Related work on affine/rigid handling in LDDMM (Glaunès 2005, Younes 2006/2010) is cited; a one-sentence contrast with asymptotically affine kernels vs exact se(d) null-space would clarify novelty for readers from that line.
Circularity Check
No circularity: honest operator construction whose claimed properties follow from null-space, equivariance, and kernel regularity by standard geometric mechanics, not by fitting or self-definitional prediction.
full rationale
The paper defines the screened elasticity operator L_s = −(Id−σ²Δ)^{s−1} ∘ div ∘ sym∇ so that its null-space is exactly se(d), inherits SE(d)-equivariance from metric operations, and obtains a stratified cometric on Land_n/SE(d) by the usual right-invariant descent / CPD construction. Local rigid-motion preservation, removal of global rigid motions, ambient n-independence, and C^1 non-collision for free Hamiltonian flow (s ≥ 1 + d/2) are consequences of that definition and of standard landmark-kernel arguments, not renamings of fitted constants or uniqueness theorems imported from the authors’ prior work. Self-citations (Michor diffeomorphism groups, LDDMM lineage, completeness for positive-definite kernels) supply background Lie-group and shape-space facts; the central operator and quotient metric are derived in-place. Scale normalization is imposed separately as a holonomic constraint and is not claimed to follow from free L_s-invariance—gaps there are completeness/correctness issues, not circular reduction of a prediction to its inputs. No fitted-input-as-prediction loop, no self-definitional X⇔Y, no load-bearing uniqueness-from-authors step.
Assumptions & free parameters
free parameters (4)
- screening length σ
- Sobolev order s =
s=2 used in experiments
- Lamé parameters μ, λ
- integration step size / step count and optimizer settings
assumptions (5)
- domain assumption Right-invariant weak Riemannian metrics on Diff_A(R^d) induce metrics on landmark orbits via Riemannian submersion when the inertia operator yields a usable kernel/cometric.
- standard math Classical sym∇ / linear elasticity has null-space exactly se(d) and is SE(d)-equivariant via Euclidean metric operations.
- domain assumption Only SE(d) (not Aff(d) or full similarity) should be quotiented inside the operator; scale is fixed by restricting to R̄=1.
- domain assumption For C¹ kernels with s ≥ 1 + d/2, Hamiltonian trajectories with α(0) in se(d)(q)^{ann} exist for all time and landmarks do not collide.
- ad hoc to paper Stratified orbit space Land_n/SE(d) carries the induced cometric on annihilators as a stratified Riemannian metric suitable for geodesic shooting away from singular strata.
invented entities (2)
-
screened elasticity operator L_s and kernels k_s
-
quotient landmark shape space Land_n with metric g_{L_s} and scale slice Land_n_1
Cite this review
Pith. "Pith review of Landmark shape spaces with induced metrics." pith.science (2026). https://pith.science/paper/ZTKO7BU3
@misc{pith2026260728064,
author = {Pith},
title = {Pith review of: Landmark shape spaces with induced metrics},
year = {2026},
howpublished = {\url{https://pith.science/paper/ZTKO7BU3}},
note = {Machine review of arXiv:2607.28064}
}
read the original abstract
We present a unification of Kendall's landmark shape spaces, where rigid motions are factored out and scale fixed on landmark configurations equipped with Euclidean geometry, with landmark configuration spaces carrying Riemannian metrics descending from right-invariant Sobolev metrics on the diffeomorphism group. The resulting new landmark shape spaces achieve the defining properties of both approaches: The regularity of the descending metric prevents landmarks from colliding, the metric is defined in the ambient space independent of the number of landmarks, local rigid transformations are preserved, global rigid motions are removed, and scale fixed. To achieve this, we define a particular Sobolev-type operator, the screened elasticity operator, whose null-space consists exactly of the rigid motions, we show how this operator descends to achieve the desired geometry, and we present approaches to solving matching problems and computing geodesics numerically. The resulting construction allows the use of landmark configuration spaces with sufficiently regular metrics in applications while retaining the shape invariances that are a hallmark of Kendall's shape spaces.
Figures
Figures from the paper (6 more)
Reference graph
Works this paper leans on
-
[1]
Higher-order fractional Laplacians : An overview
Nicola Abatangelo. Higher-order fractional Laplacians : An overview. Bruno Pini Mathematical Analysis Seminar , 12(1):53--80, 2021
2021
-
[2]
Adams and Lars I
David R. Adams and Lars I. Hedberg. Function Spaces and Potential Theory . Springer Science & Business Media, November 1999
1999
-
[3]
Dmitri Alekseevsky, Andreas Kriegl, Mark Losik, and Peter W. Michor. The R iemannian geometry of orbit spaces---the metric, geodesics, and integrable systems. Publ. Math. Debrecen , 62(3-4):247--276, 2003. Dedicated to Professor Lajos Tam \'a ssy on the occasion of his 80th birthday
2003
-
[4]
velocity
Hans C Andersen. Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations. Journal of Computational Physics , 52(1):24--34, October 1983
1983
-
[5]
Martin Bauer, Martins Bruveris, and Peter W. Michor. Overview of the Geometries of Shape Spaces and Diffeomorphism Groups . J. Math. Imaging Vis. , 50(1-2):60--97, September 2014
2014
-
[6]
Bookstein
Fred L. Bookstein. Morphometric Tools for Landmark Data : Geometry and Biology . Cambridge University Press, June 1997
1997
-
[7]
Fast Large Deformation Matching with the Energy Distance Kernel , May 2025
Siwan Boufadene, Fran c ois-Xavier Vialard, and Jean Feydy. Fast Large Deformation Matching with the Energy Distance Kernel , May 2025
2025
-
[8]
Droske and M
M. Droske and M. Rumpf. A variational approach to nonrigid morphological image registration. SIAM Journal on Applied Mathematics , 64(2):668--687, 2004
2004
Show all 34 references
-
[9]
Shape Metrics Based on Elastic Deformations
Matthias Fuchs, Bert J \"u ttler, Otmar Scherzer, and Huaiping Yang. Shape Metrics Based on Elastic Deformations . Journal of Mathematical Imaging and Vision , 35(1):86--102, 2009
2009
-
[10]
Transport Par Diff\'eomorphismes de Points, de Mesures et de Courants Pour La Comparaison de Formes et l'anatomie Num\'erique
Joan Glaun \`e s. Transport Par Diff\'eomorphismes de Points, de Mesures et de Courants Pour La Comparaison de Formes et l'anatomie Num\'erique . PhD thesis, Universit\'e Paris 13, France, 2005
2005
-
[11]
General Pattern Theory : A Mathematical Study of Regular Structures
Ulf Grenander. General Pattern Theory : A Mathematical Study of Regular Structures . Oxford University Press, USA, February 1994
1994
-
[12]
Preston, and Stefan Sommer
Karen Habermann, Stephen C. Preston, and Stefan Sommer. Characterization of geodesic completeness for landmark space. Discrete and Continuous Dynamical Systems , 50(0):224--240, 2025
2025
-
[13]
Landmark matching via large deformation diffeomorphisms
Sarang Joshi and Michael Miller. Landmark matching via large deformation diffeomorphisms. Image Processing, IEEE Transactions on , 9(8):1357--1370, 2000
2000
-
[14]
D. G. Kendall. The Diffusion of Shape . Advances in Applied Probability , 9(3):428--430, 1977
1977
-
[15]
David G. Kendall. Shape Manifolds , Procrustean Metrics , and Complex Projective Spaces . Bull. London Math. Soc. , 16(2):81--121, March 1984
1984
-
[16]
Huiling Le and David G. Kendall. The Riemannian Structure of Euclidean Shape Spaces : A Novel Environment for Statistics . The Annals of Statistics , 21(3):1225--1271, September 1993
1993
-
[17]
Strongly Elliptic Systems and Boundary Integral Equations
William Charles Hector McLean. Strongly Elliptic Systems and Boundary Integral Equations . Cambridge University Press, January 2000
2000
-
[18]
Matrix-valued Kernels for Shape Deformation Analysis
Mario Micheli and Joan Alexis Glaun \`e s. Matrix-valued Kernels for Shape Deformation Analysis . Geometry, Imaging and Computing , 1(1):57--139, 2014
2014
-
[19]
Marsden and Thomas J
Jerrold E. Marsden and Thomas J. R. Hughes. Mathematical Foundations of Elasticity . Dover Publications, New York, 1994
1994
-
[20]
Peter W. Michor. Topics in Differential Geometry . American Mathematical Soc., 2008
2008
-
[21]
Peter W. Michor. Manifolds of Mappings for Continuum Mechanics . In Reuven Segev and Marcelo Epstein, editors, Geometric Continuum Mechanics , pages 3--75. Springer International, Cham, 2020
2020
-
[22]
Michor and David Mumford
Peter W. Michor and David Mumford. A zoo of diffeomorphism groups on R ^n . Annals of Global Analysis and Geometry , 44(4):529--540, December 2013
2013
-
[23]
Michor, and David Mumford
Mario Micheli, Peter W. Michor, and David Mumford. Sectional Curvature in Terms of the Cometric , with Applications to the Riemannian Manifolds of Landmarks . SIAM Journal on Imaging Sciences , 5(1):394--433, January 2012
2012
-
[24]
Marsden and Tudor S
Jerrold E. Marsden and Tudor S. Ratiu. Introduction to Mechanics and Symmetry , volume 17 of Texts in Applied Mathematics . Springer New York, New York, NY, 1999
1999
-
[25]
Mathematical Theory of Elastic and Elasto-Plastic Bodies
Jindřich Nečas and Ivan Hlaváček, editors. Mathematical Theory of Elastic and Elasto-Plastic Bodies . February 2017
2017
-
[26]
Some basic solutions in strain gradient elasticity theory of an arbitrary order
D Rogula. Some basic solutions in strain gradient elasticity theory of an arbitrary order. Archiwum Mechaniki Stosowanej , 25:43--68, 1973
1973
-
[27]
Carl Edward Rasmussen and Christopher K. I. Williams. Gaussian Processes for Machine Learning . MIT Press, Cambridge, Mass., 2008
2008
-
[28]
A nonlinear elastic shape averaging approach
Martin Rumpf and Benedikt Wirth. A nonlinear elastic shape averaging approach. SIAM Journal on Imaging Sciences , 2(3):800--833, 2009
2009
-
[29]
Stochastics of shapes and Kunita flows
Stefan Sommer, Gefan Yang, and Elizabeth Louise Baker. Stochastics of shapes and Kunita flows. Harvard Data Science Review , 2026. arXiv:2512.11676
2026 arXiv
-
[30]
Diffeomorphisms Groups and Pattern Matching in Image Analysis
Alain Trouv \'e . Diffeomorphisms Groups and Pattern Matching in Image Analysis . International Journal of Computer Vision , 28(3):213--221, July 1998
1998
-
[31]
Metamorphoses Through Lie Group Action
Alain Trouv \'e and Laurent Younes. Metamorphoses Through Lie Group Action . Foundations of Computational Mathematics , 5(2):173--198, April 2005
2005
-
[32]
Scattered Data Approximation
Holger Wendland. Scattered Data Approximation . Cambridge Monographs on Applied and Computational Mathematics . Cambridge University Press, Cambridge, 2004
2004
-
[33]
Combining geodesic interpolating splines and affine transformations
Laurent Younes. Combining geodesic interpolating splines and affine transformations. IEEE transactions on image processing: a publication of the IEEE Signal Processing Society , 15(5):1111--1119, May 2006
2006
-
[34]
Shapes and Diffeomorphisms
Laurent Younes. Shapes and Diffeomorphisms . Springer, 2010
2010
Reviewed July 31, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.