arxiv: 2605.06558 · v2 · submitted 2026-05-07 · ⚛️ physics.med-ph

Recognition: 2 theorem links

· Lean Theorem

John Equation Constraints for the 3D X-ray Transform under a Cylindrical-Spherical Mixed Parameterization: Theoretical Derivation, Experimental Validation, and Application Analysis

Shaojie Tang , Zhiwei Qiao , Xuanqin Mou

Authors on Pith no claims yet

Pith reviewed 2026-05-11 01:46 UTC · model grok-4.3

classification ⚛️ physics.med-ph

keywords John equationX-ray transform3D CTconstraint equationsmixed parameterizationdata consistencytomographyprojection data

0 comments

The pith

The John equation under cylindrical-spherical mixed parameterization yields a complete system of constraint equations for the 3D X-ray transform.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper derives the explicit form of the John equation for the three-dimensional X-ray transform by adopting a mixed parameterization: the source point in cylindrical coordinates and the ray direction in spherical coordinates. Through transformations of the partial differential operators, application of the -1 homogeneity property, and algebraic simplification, it obtains a full set of constraint equations that the projection data must obey. Under the alignment condition where the ray direction is perpendicular to the source radial vector in the horizontal plane and the ray has no tilt, these constraints reduce to simpler differential relations carrying direct physical interpretations. A sympathetic reader would care because the results supply concrete mathematical tools for checking data consistency, calibrating geometry, and handling incomplete data in three-dimensional computed tomography systems.

Core claim

Under the cylindrical-spherical mixed parameterization with source a = (s cos θ, s sin θ, z_0) and direction d = ρ (-cos β sin α, cos β cos α, sin β), the John equation transforms, via partial differential operator changes and the -1 homogeneity condition, into a system of constraint equations on the X-ray transform data. In the special case α = θ and β = 0, the equations simplify to differential relations with clear physical meanings for the projection data.

What carries the argument

The cylindrical-spherical mixed parameterization of the source point and ray direction, combined with the transformed partial differential operators and the -1 homogeneity condition applied to the John equation.

If this is right

Provides explicit differential constraints for verifying consistency of measured projection data in 3D CT systems.
Supplies mathematical relations usable for calibrating geometric parameters of the imaging setup.
Offers tools to support reconstruction algorithms when projection data are incomplete.
Creates a direct link between the abstract John equation and concrete source-detector geometries used in practice.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The simplified relations under alignment conditions could be implemented as fast numerical filters for real-time quality control during CT scans.
The same transformation approach might be applied to other common parameterizations to obtain analogous physical interpretations.
These constraints could be incorporated into iterative reconstruction methods to penalize inconsistencies and suppress artifacts.

Load-bearing premise

The chosen cylindrical-spherical mixed parameterization must capture the full geometry of the 3D X-ray transform without loss of information or introduction of singularities, and the operator transformations and homogeneity steps must preserve the original John equation exactly.

What would settle it

Direct numerical evaluation of the X-ray transform for a known simple object such as a uniform sphere, followed by substitution into the derived constraint equations under the alignment condition to check whether the relations hold to machine precision.

read the original abstract

The John equation serves as the mathematical foundation of the X-ray transform, describing the intrinsic compatibility conditions that projection data must satisfy. In this paper, within three-dimensional (3D) Euclidean space, an innovative mixed parameterization scheme is adopted: the source point is represented using cylindrical coordinates a=(s cos{\theta},s sin{\theta},z_0), and the ray direction is represented using spherical coordinates d=\{rho}(-cos\{beta}sin{\alpha},cos\{beta}cos{\alpha},sin\{beta}). The specific form of the John equation under this geometric parameterization is systematically derived. Through detailed partial differential operator transformations, application of -1 homogeneity, and algebraic simplification, a complete system of constraint equations is obtained. In particular, under the special configurations where the ray direction is perpendicular to the radial direction of the source point in the horizontal plane (i.e., the so-called alignment condition:{\alpha} = {\theta}) and the ray has no tilt (\{beta} = 0), the constraint equations simplify to differential relations with clear physical meanings. This paper not only establishes a bridge between abstract mathematical theory and concrete imaging geometry, but also provides rigorous mathematical tools for data consistency verification, geometric parameter calibration, and incomplete-data reconstruction in 3D Computed Tomography (CT) systems. The research results are of great significance for advancing the mathematical theory and practical applications of CT imaging.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This paper spells out the John equation in a mixed cylindrical-spherical parameterization for 3D X-ray data and reduces it to simpler relations under the alignment condition, giving usable constraints for CT consistency checks.

read the letter

The core contribution is a coordinate-specific form of the John equation constraints derived through operator transformations, homogeneity reduction, and algebraic simplification. Under the alignment case where α equals θ and β is zero, the equations collapse to differential relations that tie directly to horizontal-plane geometry and zero-tilt rays. That step is the part most likely to see use in practice for verifying projection data or calibrating source positions in cone-beam CT setups. The derivation follows standard PDE techniques and the stress-test found no internal contradictions or missing terms, so the central math holds up on its own terms. The parameterization choice is reasonable for covering typical 3D imaging geometries, though any singularities at the cylinder axis would need explicit handling in the full write-up. Experimental validation is claimed, but its strength depends on the metrics and test cases shown later in the manuscript; if those are only illustrative, the practical payoff stays modest. Overall the work is incremental—an explicit transcription rather than a new theorem—but the resulting constraint system is concrete enough to be checked against real scanner data. Readers who already work with consistency conditions or incomplete-data reconstruction in medical or industrial CT will get the most out of it. The paper is coherent and engages the literature properly, so it clears the bar for a serious referee. I would send it to peer review with the expectation that reviewers will ask for the full expanded equations and any validation numbers.

Referee Report

2 major / 2 minor

Summary. The manuscript derives the specific form of the John equation for the 3D X-ray transform under a mixed cylindrical-spherical parameterization (source in cylindrical coordinates, direction in spherical coordinates). It applies partial differential operator transformations, -1 homogeneity reduction, and algebraic simplification to obtain a complete system of constraint equations, with further reduction to physically interpretable differential relations under the alignment conditions α=θ and β=0. The work claims applications to data consistency verification, geometric calibration, and incomplete-data reconstruction in 3D CT.

Significance. If the coordinate transformations preserve the John equation exactly and the experimental claims hold, the parameterization could supply practical, geometry-specific tools for consistency checks and calibration in CT systems. The approach is a direct but non-trivial specialization of a known equation, potentially useful for bridging abstract theory to concrete scanner geometries.

major comments (2)

[Abstract and theoretical derivation] Abstract and theoretical derivation sections: the manuscript asserts that the final constraint equations are obtained via operator transformations and homogeneity without loss of equivalence to the original John equation, yet provides no explicit intermediate equations, transformed operators, or final constraint expressions. This absence prevents verification that the cylindrical-spherical parameterization introduces no singularities or information loss, which is load-bearing for the central claim.
[Experimental validation and application analysis] Experimental validation and application analysis sections: the paper states that experimental validation and application analysis are performed, but supplies no data, figures, quantitative metrics (e.g., consistency error, reconstruction quality), or comparison against standard parameterizations. This undermines the asserted practical utility for CT systems.

minor comments (2)

[Abstract] Abstract contains repeated LaTeX artifacts (e.g., {alpha}, {beta}, {theta}) that impair readability; these should be rendered as proper math symbols.
[Title and abstract] The title and abstract refer to 'John Equation' with inconsistent capitalization; standardize to 'John's equation' or 'the John equation' throughout.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments on our manuscript. We address each major comment below and will incorporate the necessary revisions to improve clarity and substantiation of our claims.

read point-by-point responses

Referee: [Abstract and theoretical derivation] Abstract and theoretical derivation sections: the manuscript asserts that the final constraint equations are obtained via operator transformations and homogeneity without loss of equivalence to the original John equation, yet provides no explicit intermediate equations, transformed operators, or final constraint expressions. This absence prevents verification that the cylindrical-spherical parameterization introduces no singularities or information loss, which is load-bearing for the central claim.

Authors: We agree that the current manuscript describes the derivation process at a high level but does not display the explicit intermediate equations, the transformed partial differential operators, or the final constraint expressions. This limits independent verification of equivalence to the original John equation and the absence of singularities or information loss under the mixed parameterization. In the revised version, we will expand the theoretical derivation section to include all intermediate steps, the explicit forms of the transformed operators, the application of the -1 homogeneity reduction, and the resulting complete system of constraint equations. revision: yes
Referee: [Experimental validation and application analysis] Experimental validation and application analysis sections: the paper states that experimental validation and application analysis are performed, but supplies no data, figures, quantitative metrics (e.g., consistency error, reconstruction quality), or comparison against standard parameterizations. This undermines the asserted practical utility for CT systems.

Authors: We acknowledge that the manuscript claims experimental validation and application analysis but does not provide supporting data, figures, quantitative metrics, or comparisons. This weakens the demonstration of practical utility for data consistency verification, calibration, and reconstruction. In the revision, we will add a dedicated experimental section with simulation or phantom data results, including figures, metrics such as consistency errors and reconstruction quality (e.g., RMSE, PSNR), and direct comparisons against standard parameterizations to substantiate the claims. revision: yes

Circularity Check

0 steps flagged

No significant circularity: standard coordinate transformation of known equation

full rationale

The derivation consists of applying partial differential operator transformations, -1 homogeneity, and algebraic simplification to the pre-existing John equation under a cylindrical-spherical mixed parameterization. These are explicit change-of-variables steps on an externally known PDE; no parameters are fitted to data, no predictions are generated from subsets of the same data, and no self-citation chain is invoked as load-bearing justification. The alignment-condition simplifications are direct substitutions, not self-definitional. The work is therefore self-contained against external mathematical benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The derivation starts from the known John equation in 3D Euclidean space and applies coordinate changes and homogeneity; no free parameters, new physical entities, or ad-hoc assumptions beyond standard PDE calculus are indicated in the abstract.

axioms (2)

domain assumption The John equation holds as the intrinsic compatibility condition for the 3D X-ray transform.
This is the starting mathematical object that is transformed.
domain assumption The chosen cylindrical-spherical parameterization is a valid global or locally valid chart for source points and ray directions.
Invoked when the source is written as (s cos θ, s sin θ, z0) and direction as spherical angles.

pith-pipeline@v0.9.0 · 5570 in / 1288 out tokens · 38535 ms · 2026-05-11T01:46:37.147378+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Through detailed partial differential operator transformations, application of -1 homogeneity, and algebraic simplification, a complete system of constraint equations is obtained.
IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

the source point is represented using cylindrical coordinates a=(s cos θ, s sin θ, z₀), and the ray direction is represented using spherical coordinates d=ρ(-cos β sin α, cos β cos α, sin β)

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

26 extracted references · 26 canonical work pages

[1]

On the determination of functions by their integral values along certain manifolds

Radon J. On the determination of functions by their integral values along certain manifolds. Ber. Saechsibbe Akad. Wiss., 1917, 69: 262-278. Translated by Analogic Corp., 1976

work page 1917
[2]

The mathematics of computerized tomography

Natterer, F. The mathematics of computerized tomography. Society for Industrial and Applied Mathematics. 2001

work page 2001
[3]

The ultrahyperbolic differential equation with four independent variables

John, F. The ultrahyperbolic differential equation with four independent variables. Duke Mathematical Journal, 1938, 4(2), 300-322

work page 1938
[4]

The Radon transform on Euclidean spaces, compact two point homogeneous spaces and Grassmann manifolds

Helgason S. The Radon transform on Euclidean spaces, compact two point homogeneous spaces and Grassmann manifolds. Acta Mathematica, 1965, 113(1): 153-180

work page 1965
[5]

The Radon transform on Euclidean space

Ludwig D. The Radon transform on Euclidean space. Commun. Pure Appl. Math., 1966, XIX: 49-81

work page 1966
[6]

Integral Invariants for Computed Tomography

Wei Y, Yu H, Wang G. Integral Invariants for Computed Tomography. IEEE Signal Processing Letters, 2006, 13(9): 549-552

work page 2006
[7]

Data Consistency Conditions for Cone - Beam Projections on a Circular Trajectory

Clackdoyle R, Desbat L, Lesaint J, et al. Data Consistency Conditions for Cone - Beam Projections on a Circular Trajectory. IEEE Signal Processing Letters, 2016, 23(12): 1746-1750

work page 2016
[8]

The mathematical equivalence of consistency conditions in the divergent -beam computed tomography

Tang S, Xu Q, Mou X, et al. The mathematical equivalence of consistency conditions in the divergent -beam computed tomography. Journal of X -Ray Science and Technology, 2012, 20(1): 45-68

work page 2012
[9]

Computation of Unmeasured Third -Generation VCT Views From Measured Views

Patch S K. Computation of Unmeasured Third -Generation VCT Views From Measured Views. IEEE Transactions on Medical Imaging, 2002, 21(7): 801-813

work page 2002
[10]

A general cone -beam reconstruction algorithm

Wang G, Lin TH, Cheng P C, Shinozaki D M. A general cone -beam reconstruction algorithm. IEEE Trans. on Medical Imaging, 1993, 12(3): 486-495

work page 1993
[11]

Theoretically exact filtered backprojection -type inversion algorithm for spiral CT

Katsevich, A. Theoretically exact filtered backprojection -type inversion algorithm for spiral CT. SIAM Journal on Applied Mathematics, 2002, 62(6), 2012-2026

work page 2002
[12]

Exact image reconstruction on PI -lines from minimum data in helical cone-beam CT

Zou, Y., Pan, X. Exact image reconstruction on PI -lines from minimum data in helical cone-beam CT. Physics in Medicine & Biology, 2004, 49(6), 941

work page 2004
[13]

Advanced single -slice rebinning in cone - beam spiral CT

Kachelrieß M, Schaller S, Kalender W. Advanced single -slice rebinning in cone - beam spiral CT. Med. Phys., 2000, 27: 754-772

work page 2000
[14]

A three -dimensional weighted cone beam filtered backprojection (CB -FBP) algorithm for image reconstruction in volumetric CT under a circular source trajectory

Tang, X, Hsiech J, Hagiwara A, et al. A three -dimensional weighted cone beam filtered backprojection (CB -FBP) algorithm for image reconstruction in volumetric CT under a circular source trajectory. Phys. Med. Biol., 2005, 50(16): 3889-3905

work page 2005
[15]

A three-dimensional weighted cone beam filtered backprojection (CB -FBP) algorithm for image reconstruction in volumetric CT - helical scanning

Tang X, Hsieh J, Nilsen R A, et al. A three-dimensional weighted cone beam filtered backprojection (CB -FBP) algorithm for image reconstruction in volumetric CT - helical scanning. Phys. Med. Biol., 2006, 51(4): 855-874

work page 2006
[16]

Mou X and Duan J. Exploring the redundancy of Radon transform using a set of partial derivative equations: could we precisely reconstruct the image from a sparse - view projection without any image prior? Phys. Med. Biol. 70: 115011, 2025

work page 2025
[17]

Statistical Projection Completion in X -ray CT Using Consistency Conditions

Xu J, Taguchi K, Tsui B M W. Statistical Projection Completion in X -ray CT Using Consistency Conditions. IEEE Transactions on Medical Imaging, 2010, 29(8): 1528- 1540

work page 2010
[18]

Sinogram recovery with the method of convex projections for limited-data reconstruction in computed tomography

Kudo H, Saito T. Sinogram recovery with the method of convex projections for limited-data reconstruction in computed tomography. J. Opt. Soc. Am., A, 1991, 8(7): 1148-1160

work page 1991
[19]

Data Consistency Based Rigid Motion Artifact Reduction in Fan - Beam CT

Yu H, Wang G. Data Consistency Based Rigid Motion Artifact Reduction in Fan - Beam CT. IEEE Trans. Medical Imaging, 2007, 26(2): 249-260

work page 2007
[20]

Data consistency condition –based beam -hardening correction

Tang S, Mou X, Xu Q, et al. Data consistency condition –based beam -hardening correction. Optical Engineering, 50(7): 076501, 2011

work page 2011
[21]

Optimization Based Beam -hardening Correction in CT under Data Integral Invariant Constraint

Tang S, Huang K, Cheng Y, et al. Optimization Based Beam -hardening Correction in CT under Data Integral Invariant Constraint. Phys. Med. Biol., 2018, 63(13): 135015

work page 2018
[22]

Non -convex Optimization based Optimal Bone Correction for Various Beam -hardening Artifacts in CT Imaging

Tang S, Huang T, Qiao Z, et al. Non -convex Optimization based Optimal Bone Correction for Various Beam -hardening Artifacts in CT Imaging. Journal of X -Ray Science and Technology, 2022, 30(4): 805-822

work page 2022
[23]

Projection correlation based view interpolation for cone beam CT: primary fluence restoration in scatter measurement with a moving beam stop array

Yan H, Mou X, Tang S, et al. Projection correlation based view interpolation for cone beam CT: primary fluence restoration in scatter measurement with a moving beam stop array. Phys. Med. Biol., 2010, 55: 6353-6375

work page 2010
[24]

Computed Tomography (From Photon Statistics to Modern Cone -Beam CT)

Buzug T M. Computed Tomography (From Photon Statistics to Modern Cone -Beam CT) . Springer-Verlag Berlin Heidelberg, 2008

work page 2008
[25]

Tang X. Spectral Multi -Detector Computed Tomography (sMDCT): Data Acquisition, Image Formation, Quality Assessment and Contrast Enhancement (Series in Medical Physics and Biomedical Engineering). CRC Press, 2023

work page 2023
[26]

Calibration for Circular Cone-Beam CT Based on Consistency Conditions

Jé rô me Lesaint; Simon Rit ; Rolf Clackdoyle ; Laurent Desbat . Calibration for Circular Cone-Beam CT Based on Consistency Conditions. IEEE Transactions on Radiation and Plasma Medical Sciences, 2017, 1(6): 517-526. Appendices Appendix A: Derivation of the John Equation Based on Ray Direction From the definition d = b - a, we get b = a + d, and g(a, d) =...

work page 2017