Sampling theorems for inverse problems on Riemannian manifolds
Pith reviewed 2026-05-18 22:59 UTC · model grok-4.3
The pith
Sampling theorems give explicit reconstruction error bounds for inverse problems on Riemannian manifolds from noisy pointwise samples.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We consider inverse problems of reconstructing f from y = Ff + noise on a Riemannian manifold M without boundary, using only pointwise samples y_j = (Ff)(x_j) + η_j where {x_j} is a Marcinkiewicz-Zygmund family. Sampling theorems are derived that supply explicit bounds on the reconstruction error in terms of n, the smoothness of f and properties of F. The case of F being a convolution on a compact two-point homogeneous space is studied in detail, yielding as a corollary a sampling theorem for convolutions on the two-dimensional sphere, illustrated with four examples of relevant measurements.
What carries the argument
Marcinkiewicz-Zygmund families of sampling points that satisfy the inequalities needed to control reconstruction error for the operator F.
If this is right
- The error bound decreases explicitly with larger n or greater smoothness of f when the points form a Marcinkiewicz-Zygmund family.
- Convolutions on compact two-point homogeneous spaces admit detailed error estimates from the sampling theorems.
- A sampling theorem holds for convolutions on the two-dimensional sphere as a direct corollary.
- Four concrete examples of terrestrial and celestial measurements obtain explicit error guarantees from the derived bounds.
Where Pith is reading between the lines
- Similar explicit bounds could be obtained on other Riemannian manifolds once suitable Marcinkiewicz-Zygmund point families are identified.
- The dependence on n and smoothness could be used to decide how many measurement locations are needed for a target accuracy in sphere-based data applications.
- The proof techniques for the convolution case may extend to other classes of operators F on the same spaces.
Load-bearing premise
The chosen points must form a Marcinkiewicz-Zygmund family on the manifold so that the required sampling inequalities hold.
What would settle it
A calculation on the two-dimensional sphere showing that reconstruction error for a smooth test function under a known convolution operator does not improve with n at the rate given by the bounds would disprove the theorems.
Figures
read the original abstract
We consider inverse problems consisting of the reconstruction of an unknown signal $f$ from noisy measurements $y=Ff+\text{noise}$, where $Ff$ is a function on a Riemannian manifold without boundary $\mathcal M$. We consider the case when only pointwise samples are available, namely $y_j = (Ff)(x_j)+\eta_j$, where $\{x_j\}_{j=1}^n\subseteq\mathcal M$ is a Marcinkiewicz-Zygmund family. We derive sampling theorems providing explicit bounds on the reconstruction error depending on $n$, the smoothness of $f$ and the properties of $F$. We study in detail the case when $F$ is a convolution on a compact two-point homogeneous space. As a corollary, we state a sampling theorem for convolutions on the two-dimensional sphere, and discuss four relevant examples related to terrestrial and celestial measurements.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops sampling theorems for inverse problems on Riemannian manifolds without boundary. Given noisy pointwise samples y_j = (Ff)(x_j) + η_j where {x_j} forms a Marcinkiewicz-Zygmund family, it derives explicit bounds on the reconstruction error in terms of n, the smoothness of f, and properties of the operator F. Detailed analysis is provided for convolution operators F on compact two-point homogeneous spaces, yielding a corollary for convolutions on the 2-sphere together with four examples drawn from terrestrial and celestial measurements.
Significance. If the claimed explicitness of the error bounds holds without hidden manifold-dependent constants, the results would offer a useful theoretical tool for stable reconstruction from discrete samples in manifold signal processing. The specialization to two-point homogeneous spaces and the sphere corollary provide concrete, applicable statements that could inform measurement design in geophysics and astronomy. The manuscript's strength lies in linking abstract sampling inequalities to operator-specific bounds, but this value is conditional on the quantification of the underlying Marcinkiewicz-Zygmund constants.
major comments (2)
- [Abstract / MZ-family paragraph] Abstract and the paragraph introducing the Marcinkiewicz-Zygmund family: the reconstruction-error bounds are stated to depend explicitly on n, smoothness of f, and properties of F, yet they rest on sampling inequalities whose lower and upper constants are not shown to be expressible in terms of n and intrinsic geometric invariants of a general M (e.g., injectivity radius or curvature bounds). If these constants remain unspecified or absorbed, the explicitness claim does not hold for arbitrary Riemannian manifolds even though it may hold for the symmetric convolution cases.
- [General sampling theorem] Section deriving the general sampling theorem (presumably §3 or Theorem 3.1): the error estimate is obtained by combining the sampling inequality with an approximation or regularization term; the final expression must be checked to confirm that no manifold-specific constants from the MZ family enter the bound in a non-explicit way, otherwise the dependence claimed in the abstract is incomplete.
minor comments (2)
- [Introduction] Notation for the manifold M and the operator F should be introduced consistently in the first section rather than appearing first in the abstract.
- [Sphere corollary / examples] The four examples on the sphere would benefit from a short table summarizing the corresponding operators F, smoothness assumptions, and resulting error rates.
Simulated Author's Rebuttal
We thank the referee for the careful reading and insightful comments on our manuscript. The points raised concern the precise meaning of 'explicit' bounds in the general Riemannian setting versus the specialized convolution cases. We address each major comment below and indicate the revisions we will make to clarify the dependence on Marcinkiewicz-Zygmund constants.
read point-by-point responses
-
Referee: [Abstract / MZ-family paragraph] Abstract and the paragraph introducing the Marcinkiewicz-Zygmund family: the reconstruction-error bounds are stated to depend explicitly on n, smoothness of f, and properties of F, yet they rest on sampling inequalities whose lower and upper constants are not shown to be expressible in terms of n and intrinsic geometric invariants of a general M (e.g., injectivity radius or curvature bounds). If these constants remain unspecified or absorbed, the explicitness claim does not hold for arbitrary Riemannian manifolds even though it may hold for the symmetric convolution cases.
Authors: We agree that the abstract's phrasing could be read as claiming fully explicit constants solely in terms of n and basic geometric invariants for arbitrary manifolds. In the general theorem the error bounds are expressed explicitly in terms of the lower and upper MZ constants of the given sampling family together with n, the smoothness parameter, and operator properties of F. These MZ constants are intrinsic to the chosen point set on M and are not further reduced to curvature bounds in the general case, as that would require additional assumptions on the distribution of points. In the subsequent analysis for compact two-point homogeneous spaces we do derive explicit upper and lower bounds on the MZ constants using the representation theory and volume growth of those spaces, yielding completely explicit expressions for the sphere corollary. We will revise the abstract and the introductory paragraph on MZ families to state that the general bounds depend explicitly on the MZ constants of the sampling set, while the convolution results on two-point homogeneous spaces provide fully explicit constants in terms of n and the geometry. revision: partial
-
Referee: [General sampling theorem] Section deriving the general sampling theorem (presumably §3 or Theorem 3.1): the error estimate is obtained by combining the sampling inequality with an approximation or regularization term; the final expression must be checked to confirm that no manifold-specific constants from the MZ family enter the bound in a non-explicit way, otherwise the dependence claimed in the abstract is incomplete.
Authors: We have re-examined the proof of the general sampling theorem. The reconstruction error is bounded by a term involving the MZ lower constant times the noise level plus a term involving the MZ upper constant times the approximation error of the regularized solution. No additional manifold-dependent constants are introduced beyond those already present in the MZ inequality itself. The final expression is therefore explicit once the MZ constants are fixed by the sampling family. To address the referee's concern we will add a short remark after the statement of the general theorem clarifying that the MZ constants are regarded as known quantities characterizing the sampling set and that their explicit evaluation is carried out in the later sections for the symmetric-space cases. revision: yes
Circularity Check
No circularity; bounds derived from assumed MZ sampling inequalities
full rationale
The paper assumes {x_j} forms a Marcinkiewicz-Zygmund family on M and derives sampling theorems with error bounds depending on n, smoothness of f, and properties of F. No load-bearing step reduces by the paper's own equations or self-citation to its inputs by construction. The MZ property supplies the sampling inequalities as an external assumption; constants are absorbed into the stated bounds without redefinition or fitting. Special cases (convolutions on two-point homogeneous spaces and the sphere) are handled separately with the same structure. The derivation is self-contained against the given assumptions and does not exhibit self-definitional, fitted-prediction, or uniqueness-imported patterns.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The collection {x_j} is a Marcinkiewicz-Zygmund family on the manifold M.
- domain assumption F is a bounded linear operator (specifically a convolution when studied in detail) on the manifold.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AbsoluteFloorClosure.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We derive sampling theorems providing explicit bounds on the reconstruction error depending on n, the smoothness of f and the properties of F. ... {x_j} forms a Marcinkiewicz-Zygmund family on the manifold M
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
compact two-point homogeneous space ... convolution operator F ... Fourier multiplier ... Sobolev spaces H^σ(S^d)
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Abhishake, T. Helin, and N. Mücke.Statistical inverse learning problems with random observations, pages 201–244. De Gruyter, Berlin, Boston, 2025
work page 2025
-
[2]
K. Abraham and R. Nickl. On statistical Calderón problems.Math. Stat. Learn., 2(2):165–216, 2019
work page 2019
- [3]
-
[4]
Nonlinearity, 36(1):734–808, 2023
G.S.Alberti, A.Arroyo, andM.Santacesaria.Inverseproblemsonlow-dimensionalmanifolds. Nonlinearity, 36(1):734–808, 2023
work page 2023
-
[5]
G. S. Alberti, A. Felisi, M. Santacesaria, and S. I. Trapasso. Compressed sensing for inverse problems and the sample complexity of the sparse radon transform.Journal of the European Mathematical Society, 2025
work page 2025
- [6]
-
[7]
G. S. Alberti and M. Santacesaria. Infinite dimensional compressed sensing from anisotropic measurements and applications to inverse problems in PDE.Appl. Comput. Harmon. Anal., 50:105–146, 2021
work page 2021
-
[8]
G. S. Alberti and M. Santacesaria. Infinite-dimensional inverse problems with finite measurements.Arch. Ration. Mech. Anal., 243(1):1–31, 2022
work page 2022
-
[9]
J. J. Benedetto and P. J. S. G. Ferreira, editors.Modern sampling theory. Applied and Numerical Harmonic Analysis. Birkhäuser Boston, Inc., Boston, MA, 2001. Mathematics and applications
work page 2001
-
[10]
M. Bertero, P. Boccacci, and C. De Mol.Introduction to inverse problems in imaging. CRC press, 2021
work page 2021
- [11]
-
[12]
G. Blanchard and N. Mücke. Optimal rates for regularization of statistical inverse learning problems.Found. Comput. Math., 18(4):971–1013, 2018
work page 2018
-
[13]
L. Brandolini, B. Gariboldi, and G. Gigante. Irregularities of Distribution on Two-Point Homogeneous Spaces. In The Mathematical Heritage of Guido Weiss, Appl. Numer. Harmon. Anal., pages 101–125. Birkhäuser/Springer, Cham, 2025
work page 2025
-
[14]
L. Brandolini, B. Gariboldi, G. Gigante, and A. Monguzzi. Single radius spherical cap discrepancy on compact two-point homogeneous spaces.Math. Z., 308(3):Paper No. 52, 16, 2024
work page 2024
-
[15]
T. A. Bubba, M. Burger, T. Helin, and L. Ratti. Convex regularization in statistical inverse learning problems. Inverse Probl. Imaging, 17(6):1193–1225, 2023
work page 2023
-
[16]
T. A. Bubba and L. Ratti. Shearlet-based regularization in statistical inverse learning with an application to x-ray tomography.Inverse Problems, 38(5):Paper No. 054001, 43, 2022
work page 2022
- [17]
-
[18]
A. Ebner and M. Haltmeier. Convergence rates for the joint solution of inverse problems with compressed sensing data.Inverse Problems, 39(1):Paper No. 015011, 16, 2023
work page 2023
- [19]
-
[20]
H. W. Engl, M. Hanke, and A. Neubauer.Regularization of inverse problems, volume 375 ofMathematics and its Applications. Kluwer Academic Publishers Group, Dordrecht, 1996
work page 1996
-
[21]
F. Filbir and H. N. Mhaskar. A quadrature formula for diffusion polynomials corresponding to a generalized heat kernel.J. Fourier Anal. Appl., 16(5):629–657, 2010
work page 2010
-
[22]
F. Filbir and H. N. Mhaskar. Marcinkiewicz-Zygmund measures on manifolds.J. Complexity, 27(6):568– 596, 2011
work page 2011
-
[23]
G. B. Folland.A course in abstract harmonic analysis. Textbooks in Mathematics. CRC Press, Boca Raton, FL, second edition, 2016
work page 2016
-
[24]
K.-J. Förster. Inequalities for ultraspherical polynomials and application to quadrature. InProceedings of the Seventh Spanish Symposium on Orthogonal Polynomials and Applications (VII SPOA) (Granada, 1991), volume 49, pages 59–70, 1993
work page 1991
- [25]
-
[26]
B. Gariboldi and G. Gigante. Optimal asymptotic bounds for designs on manifolds.Anal. PDE, 14(6):1701– 1724, 2021
work page 2021
-
[27]
G. Gigante and P. Leopardi. Diameter bounded equal measure partitions of Ahlfors regular metric measure spaces. Discrete Comput. Geom., 57(2):419–430, 2017
work page 2017
-
[28]
M. Giordano and R. Nickl. Consistency of Bayesian inference with Gaussian process priors in an elliptic inverse problem.Inverse Problems, 36(8):085001, 35, 2020
work page 2020
- [29]
-
[30]
D. M. Healy, Jr., H. Hendriks, and P. T. Kim. Spherical deconvolution.J. Multivariate Anal., 67(1):1–22, 1998
work page 1998
-
[31]
Helgason.Groups and geometric analysis, volume 83 ofMathematical Surveys and Monographs
S. Helgason.Groups and geometric analysis, volume 83 ofMathematical Surveys and Monographs. Amer- ican Mathematical Society, Providence, RI, 2000. Integral geometry, invariant differential operators, and spherical functions, Corrected reprint of the 1984 original
work page 2000
-
[32]
T. Helin. Least squares approximations in linear statistical inverse learning problems.SIAM J. Numer. Anal., 62(4):2025–2047, 2024
work page 2025
-
[33]
E. Herrholz and G. Teschke. Compressive sensing principles and iterative sparse recovery for inverse and ill-posed problems.Inverse Problems, 26(12):125012, 24, 2010
work page 2010
-
[34]
R. Hielscher and M. Quellmalz. Optimal mollifiers for spherical deconvolution. Inverse Problems, 31(8):085001, 28, 2015. SAMPLING THEOREMS FOR INVERSE PROBLEMS ON RIEMANNIAN MANIFOLDS 29
work page 2015
-
[35]
Hörmander.The analysis of linear partial differential operators
L. Hörmander.The analysis of linear partial differential operators. III. Classics in Mathematics. Springer, Berlin, 2007. Pseudo-differential operators, Reprint of the 1994 edition
work page 2007
-
[36]
Isakov.Inverseproblems for partial differentialequations, volume 127 ofApplied Mathematical Sciences
V. Isakov.Inverseproblems for partial differentialequations, volume 127 ofApplied Mathematical Sciences. Springer, Cham, third edition, 2017
work page 2017
-
[37]
Statistical inverse problems: discretization, model reductionand inversecrimes
J.Kaipio and E.Somersalo. Statistical inverse problems: discretization, model reductionand inversecrimes. J. Comput. Appl. Math., 198(2):493–504, 2007
work page 2007
-
[38]
B. Kaltenbacher, A. Neubauer, and O. Scherzer.Iterative regularization methods for nonlinear ill-posed problems, volume 6 ofRadonSeriesonComputationalandAppliedMathematics. Walter de Gruyter GmbH & Co. KG, Berlin, 2008
work page 2008
-
[39]
H. Kekkonen, M. Lassas, and S. Siltanen. Analysis of regularized inversion of data corrupted by white Gaussian noise.Inverse Problems, 30(4):045009, 18, 2014
work page 2014
-
[40]
G. Kerkyacharian, T. M. Pham Ngoc, and D. Picard. Localized spherical deconvolution.Ann. Statist., 39(2):1042–1068, 2011
work page 2011
-
[41]
P. T. Kim and J.-Y. Koo. Optimal spherical deconvolution.J. Multivariate Anal., 80(1):21–42, 2002
work page 2002
-
[42]
H. J. Landau. Necessary density conditions for sampling and interpolation of certain entire functions.Acta Mathematica, 117:37–52, 1967
work page 1967
- [43]
-
[44]
D. Lazzaro, S. Morigi, and L. Ratti. Oracle-net for nonlinear compressed sensing in electrical impedance tomography reconstruction problems.J. Sci. Comput., 101(2):Paper No. 49, 29, 2024
work page 2024
-
[45]
J. Li, Y. Ling, J. Geng, and H. Wang. Weighted leastℓp approximation on compact Riemannian manifolds. J. Fourier Anal. Appl., 30(5):Paper No. 57, 29, 2024
work page 2024
- [46]
-
[47]
M. Maggioni and H. N. Mhaskar. Diffusion polynomial frames on metric measure spaces.Appl. Comput. Harmon. Anal., 24(3):329–353, 2008
work page 2008
-
[48]
J. D. McEwen and Y. Wiaux. A novel sampling theorem on the sphere.IEEE Trans. Signal Process., 59(12):5876–5887, 2011
work page 2011
-
[49]
H. N. Mhaskar, F. J. Narcowich, and J. D. Ward. Spherical Marcinkiewicz-Zygmund inequalities and positive quadrature.Math. Comp., 70(235):1113–1130, 2001
work page 2001
-
[50]
F. c. Monard and P. Stefanov. Sampling the X-ray transform on simple surfaces.SIAM J. Math. Anal., 55(3):1707–1736, 2023
work page 2023
-
[51]
Natterer.The mathematics of computerized tomography, volume 32 ofClassics in Applied Mathematics
F. Natterer.The mathematics of computerized tomography, volume 32 ofClassics in Applied Mathematics. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2001. Reprint of the 1986 original
work page 2001
-
[52]
Nickl.Bayesiannon-linear statistical inverseproblems
R. Nickl.Bayesiannon-linear statistical inverseproblems. Zurich Lectures in Advanced Mathematics. EMS Press, Berlin, [2023]©2023
work page 2023
-
[53]
R. Nickl and G. P. Paternain. On some information-theoretic aspects of non-linear statistical inverse prob- lems. InICM—International Congress of Mathematicians. Vol. 7. Sections 15–20, pages 5516–5538. EMS Press, Berlin, [2023]©2023
work page 2023
-
[54]
H. Nyquist. Certain topics in telegraph transmission theory.Transactions of the American Institute of Electrical Engineers, 47(2):617–644, 1928
work page 1928
-
[55]
R. Penrose. On best approximation solutions of linear matrix equations.Proc. Cambridge Philos. Soc., 52:17–19, 1956
work page 1956
-
[56]
E. A. Rakhmanov, E. B. Saff, and Y. M. Zhou. Minimal discrete energy on the sphere.Math. Res. Lett., 1(6):647–662, 1994
work page 1994
-
[57]
R. Ramlau and O. Scherzer, editors.The Radon transform—the first 100 years and beyond, volume 22 of Radon Series on Computational and Applied Mathematics. Walter de Gruyter GmbH & Co. KG, Berlin, [2019]©2019. 30 G. S. ALBERTI, E. DE VITO, B. GARIBOLDI, AND G. GIGANTE
work page 2019
-
[58]
A. Rastogi, G. Blanchard, and P. Mathé. Convergence analysis of Tikhonov regularization for non-linear statistical inverse problems.Electron. J. Stat., 14(2):2798–2841, 2020
work page 2020
-
[59]
O. Scherzer, editor. Handbook of mathematical methods in imaging. Vol. 1, 2, 3. Springer, New York, second edition, 2015
work page 2015
-
[60]
C. E. Shannon. Communication in the presence of noise.Proceedings of the IRE, 37(1):10–21, 1949
work page 1949
-
[61]
M. Simeoni. Functional penalised basis pursuit on spheres.Appl. Comput. Harmon. Anal., 53:1–53, 2021
work page 2021
-
[62]
M. M. J.-A. Simeoni.FunctionalInverseProblems on Spheres: Theory,Algorithms and Applications. PhD thesis, EPFL, Lausanne, 2020
work page 2020
-
[63]
M. M. Skriganov. Point distributions in two-point homogeneous spaces.Mathematika, 65(3):557–587, 2019
work page 2019
- [64]
-
[65]
P. Stefanov and S. Tindel. Sampling linear inverse problems with noise.Asymptot. Anal., 132(3-4):331–382, 2023
work page 2023
-
[66]
E. M. Stein and G. Weiss.Introduction to Fourier analysis on Euclidean spaces, volume No. 32 ofPrinceton Mathematical Series. Princeton University Press, Princeton, NJ, 1971
work page 1971
-
[67]
Szegő.Orthogonal polynomials, volume Vol
G. Szegő.Orthogonal polynomials, volume Vol. XXIII of American Mathematical Society Colloquium Publications. American Mathematical Society, Providence, RI, fourth edition, 1975
work page 1975
-
[68]
A. C. M. van Rooij and F. H. Ruymgaart. Regularized deconvolution on the circle and the sphere. In Nonparametric functional estimation and related topics (Spetses, 1990), volume 335 ofNATO Adv. Sci. Inst. Ser. C: Math. Phys. Sci., pages 679–690. Kluwer Acad. Publ., Dordrecht, 1991
work page 1990
- [69]
-
[70]
H.-C. Wang. Two-point homogeneous spaces.Ann. of Math. (2), 55:177–191, 1952
work page 1952
-
[71]
E. T. Whittaker. On the functions which are represented by the expansions of the interpolation theory. Proceedings of the Royal Society of Edinburgh, 35:181–194, 1915
work page 1915
-
[72]
J. A. Wolf.Spaces of constant curvature. University of California, Department of Mathematics, Berkeley, CA, second edition, 1972. 1MaLGa Center, Department of Mathematics, University of Genoa, Via Dodecaneso 35, 16146, Genov a, Italy Email address:giovanni.alberti@unige.it Email address:ernesto.devito@unige.it 2Dipartimento di Ingegneria Gestionale, dell’...
work page 1972
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.