REVIEW 3 major objections 4 minor 23 references
An Inversion Formula for Horizontal Conical Radon Transform
T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read A function in 3D can be reconstructed exactly from integrals over all cones with horizontal axes and vertices on a line, for every integer weight.
desk verdict A useful new inversion formula for a practical cone transform, but the proof of Lemma 2.5 drops boundary terms that appear to break the identity. 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 proof is carried by three transforms and one connecting identity. The vertical slice transform $\Gamma g(\phi,t)$ integrates a function on the sphere over the circle where $e_\phi\cdot\omega=t$, and its spherical inversion (Theorem 2.3) is the first reconstruction step. The weighted X-ray transform $\chi_k f(z,\omega)$ integrates $f$ along the ray from the vertex $b(z)$ in direction $\omega$ with weight $r^k$. Lemma 2.4 connects them: $\Gamma(\chi_k f)(z,\beta,\cos\psi)=T_k(f)(z,\beta,\psi)/(2\pi\sin^2\psi)$, so cone data give slice data of the weighted X-ray transform. Because slice circles are symmetric about the horizontal plane, only the even part $(\chi_k f)_e$ survives, and an identity expresses the unweighted even X-ray data $(\chi_0 f)_e$ as an integral of $z$-derivatives of $(\chi_k f)_e$, reducing any weight $k$ to the unweighted case. The even X-ray data are then read as V-line integrals in each vertical half-plane, and two-dimensional Radon inversion finishes the reconstruction.
What would settle it
Take a compactly supported smooth $f$ whose conical transform $T_k$ can be computed to high accuracy, evaluate the right-hand side of Theorem 2.6 at a point outside the support of $f$, and check that the principal-value integral is zero; a nonzero value, or divergence as the regularization is removed, would falsify the claimed inversion.
Extended reading notes
Core claim
The central claim, Theorem 2.6, is that for every integer $k\ge 1$ and every compactly supported smooth $f$, the value $f(x)$ at $x=(r\,e_\varphi,x_3)$ is reproduced by the displayed integral of the conical data $T_k(f)(p/\sin\eta,\beta,\psi)$ and its derivatives, with the singular kernels $1/(\cos\psi-\sin\gamma\cos(\beta-\varphi))$ and $1/(r\cos\eta+x_3\sin\eta-p)$. The integrations run over the opening angle, the detector direction, and the cone-vertex position. Thus the three-dimensional family of conical projections with line vertices and horizontal axes determines $f$ exactly, not merely approximately. The proof converts cone integrals into vertical-slice data on the sphere, recovers the even weighted X-ray transform, and inverts the resulting V-line data by the two-dimensional Radon transform.
Load-bearing premise
The derivation assumes that the singular integrals in the inversion formula are interpreted as principal values and that data for obtuse cone opening angles are known through the symmetry relation; if either assumption fails, the equality in Theorem 2.6 need not hold.
Editorial extensions
If this is right
- For every positive integer $k$, a compactly supported smooth function in $\mathbb{R}^3$ is determined exactly by its horizontal conical Radon transform on the three-parameter family of cones with vertices on a line.
- Compton-camera setups with detectors on a line acquire a direct analytic reconstruction route: recover even weighted X-ray data, then invert a planar Radon transform slice by slice.
- The inversion separates over the azimuth $\varphi$, so the reconstruction can be carried out independently for each vertical half-plane.
- The formula handles every integer weight $k$ with the same structure, so stronger radial weighting needs no new inversion theory.
Reading between the lines
- The $\psi$-integral in Theorem 2.6 runs to $\pi$ although $T_k$ is defined only for $\psi\in(0,\pi/2)$; the paper leaves implicit that obtuse-angle data must come from the symmetry $T_k(z,\beta,\pi-\psi)=T_k(z,\beta+\pi,\psi)$, so a detector recording only acute cones would need additional measurements or an extrapolation step.
- The two kernel denominators are singular inside the integration domain, and the paper states no principal-value prescription; any numerical implementation must choose and validate a regularization before the formula can handle measured data.
- Because the proof factors through the vertical slice transform, the same strategy could plausibly invert other three-parameter cone manifolds with horizontal axes, such as vertices lying on a closed curve; the paper does not pursue those geometries.
- The inversion differentiates the data $k+1$ times in $p$ and once in $\psi$, so noise sensitivity is a likely practical obstacle; simulated-phantom tests of the formula's stability would be a natural next step.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper considers the weighted conical Radon transform T_k(f)(z,beta,psi) over all one-sided circular cones whose vertices lie on the z-axis and whose symmetry axes are horizontal, with weight ||x-b(z)||^{k-1}. It claims an explicit inversion formula: for compactly supported smooth f in R^3, f(x) is recovered by a triple integral involving T_k and its derivatives. The strategy is to relate T_k to the vertical slice transform of the weighted X-ray transform (Lemma 2.4), to invert the vertical slice transform using Gindikin's theorem (Lemma 2.5), and then to invert the resulting V-line transform via the two-dimensional Radon transform (Theorem 2.6).
Significance. If the main theorem were correct, it would provide an exact inversion formula for a three-dimensional family of cones, corresponding to Compton camera data with detectors on a line, and would complement existing results for other cone manifolds. The paper is concise and the overall reduction to established inversions is appealing. The derivation has no fitted parameters, and the reliance on known inversion theorems (Gindikin's vertical slice inversion and the Radon inversion) is transparent. However, the central derivation contains a serious omitted-boundary-term issue in Lemma 2.5, and the statement of Theorem 2.6 also leaves the domain of psi and the interpretation of singular integrals unspecified. These issues are load-bearing for the claimed formula.
major comments (3)
- [§2, Lemma 2.5, proof of (4)] The passage from the q-integral to the psi-integral silently discards boundary terms. The function Gamma(chi_k f)_e is extended by zero for |t|>1, so its distributional derivative in q contains delta contributions at q = 1 - sin(eta) cos(beta-phi) and q = -1 - sin(eta) cos(beta-phi), corresponding to t = 1 and t = -1. The change of variables cos(psi) = sin(eta) cos(beta-phi) + q with psi in [0,pi] captures only the classical interior derivative and omits these endpoint deltas. The omitted contribution is of the form involving chi_k f(z, e_beta) and chi_k f(z, -e_beta), i.e. weighted X-ray data along horizontal rays. Such data are not present as an additive term in the displayed integrand of (4), and for a generic compactly supported f (for example, a positive bump on the x-axis) the omitted term is nonzero. Unless a cancellation is proved, or the boundary term is explicitly expressed through T_k (for instance via the limit of T_k/sin(psi) as psi -> 0+), formula (4) and hence Theorem 2.6 are incomplete.
- [§2, Theorem 2.6] The inversion formula integrates psi over [0,pi], whereas T_k is defined only for psi in (0,pi/2). The manuscript implicitly relies on an extension of T_k to obtuse cones, presumably through the symmetry T_k(z, beta, pi-psi) = T_k(z, beta+pi, psi), which follows from the relation e_{beta+pi} = -e_beta and the definition of the cone. This extension is never stated or proved. Since the formula uses values of T_k at psi > pi/2, the statement of Theorem 2.6 is currently not well-defined and needs an explicit definition or a symmetry argument.
- [§2, Theorem 2.6 and Lemma 2.5] The singular integrals in the formula are not given a precise interpretation. The inner psi-integral has the kernel 1/(cos(psi) - sin(gamma) cos(beta-phi)), which vanishes for generic gamma, beta, phi, psi, and the outer p-integral has the kernel 1/(r cos(eta) + x_3 sin(eta) - p), which is a Hilbert-type singular integral after the derivatives act on T_k. The paper neither states that these are principal-value integrals nor proves that the displayed iterated integral converges in the ordinary or principal-value sense for f in C_0^infinity(R^3). This is essential for the formula to be meaningful as an inversion formula.
minor comments (4)
- [§2, Theorem 2.3] The statement of Gindikin's inversion formula mixes notation: the integrand contains 1/t and derivatives with respect to t, while the proof uses 1/q and derivatives with respect to q. The variables should be made consistent.
- [§1, Definition of T_k] The domain of T_k is given as R x [0,2pi) x (0,pi/2), but Lemma 2.4 states an identity for psi in (0,pi). Please reconcile the domain of definition with the range of psi used throughout the paper.
- [§2, Lemma 2.5] The notation phi in [0,2pi] should presumably be phi in [0,2pi), and the interval notation in the statement should be made uniform.
- [§2, proof of Theorem 2.6] The sentence 'The inversion such transform can be reduced to that of the X-ray transform' contains a grammatical error and should read 'The inversion of such a transform...'.
Circularity Check
No significant circularity: the inversion is a composition of independent external inversion theorems and transform identities, with no fitted inputs or self-citation chain.
full rationale
The derivation chain is not circular. T_k(f) is defined directly as a weighted cone-surface integral (Section 1). Lemma 2.4 derives an exact identity relating T_k(f) to the vertical slice transform of the weighted X-ray transform chi_k f, using only spherical-coordinate substitution and integration over the fiber. Lemma 2.5 then applies Gindikin's vertical slice inversion theorem (Theorem 2.3) to recover (chi_0 f)_e from T_k(f), and uses an independent relation from Moon-Haltmeier [17] connecting chi_0 f to chi_k f. Theorem 2.6 finishes by passing from (chi_0 f)_e to the 2D Radon transform of a reflected restriction and applying the standard Radon inversion formula [13]. Each load-bearing step invokes an established external inversion theorem or a transform identity; no target quantity f or T_k(f) is fed back into the derivation as an assumption, no fitted parameter is renamed as a prediction, and no load-bearing step is justified solely by a self-citation. The concerns noted by the reader and skeptic about the psi-range extension, principal-value interpretation, and possible omitted boundary terms in Lemma 2.5 are potential rigor or correctness issues, not circularity: even if the formula were incomplete, that would not make the argument equivalent to its inputs by construction. Thus the circularity score is 0.
Assumptions & free parameters
assumptions (5)
- standard math Gindikin's vertical slice inversion formula (Theorem 2.3) for even continuous functions on S^2.
- standard math Standard 2D Radon inversion formula f = (1/2 pi^2) integral integral partial_p R f / (x dot theta - p) dp dtheta.
- domain assumption Relation (chi_0 f)_e(z, omega) = 1/(k-1)! integral (omega_3 - s)^(k-1) partial_z^k (chi_k f)_e(z, omega_1, omega_2, s) ds.
- domain assumption The conical transform data is known for all half-angles psi in (0, pi), or can be obtained from acute-angle data by the unstated symmetry T_k(z, beta, pi-psi) = T_k(z, beta+pi, psi).
- domain assumption Singular integrals in the inversion formula are interpreted as principal values.
Cite this review
Pith. "Pith review of An Inversion Formula for Horizontal Conical Radon Transform." pith.science (2026). https://pith.science/paper/AWHKEZRT
@misc{pith2026190808255,
author = {Pith},
title = {Pith review of: An Inversion Formula for Horizontal Conical Radon Transform},
year = {2026},
howpublished = {\url{https://pith.science/paper/AWHKEZRT}},
note = {Machine review of arXiv:1908.08255}
}
read the original abstract
In this paper, we consider the conical Radon transform on all cones with horizontal central axis whose vertices are on a straight line. We derive an explicit inversion formula for such transform. The inversion makes use of the vertical slice transform on a sphere and V-line transform on a plane.
Figures
Reference graph
Works this paper leans on
-
[1]
Detecting small low emission radiating sources
Moritz Allmaras, David P Darrow, Yulia Hristova, Guido Kanschat, and Peter Kuchment. Detecting small low emission radiating sources. arXiv preprint arXiv:1012.3373 , 2010
arXiv 2010
-
[2]
Inversion of the v-line radon transform in a disc and its applications in imaging
Gaik Ambartsoumian. Inversion of the v-line radon transform in a disc and its applications in imaging. Computers & Mathematics with Applications , 64(3):260–265, 2012
work page 2012
-
[3]
Analytical reconstruction formula for one- dimensional compton camera
Roman Basko, Gengsheng L Zeng, and Grant T Gullberg. Analytical reconstruction formula for one- dimensional compton camera. IEEE Transactions on Nuclear Science , 44(3):1342–1346, 1997
work page 1997
-
[4]
Application of spherical harmonics to image reconstruction for the compton camera
Roman Basko, Gengsheng L Zeng, and Grant T Gullberg. Application of spherical harmonics to image reconstruction for the compton camera. Physics in Medicine & Biology , 43(4):887, 1998
work page 1998
-
[5]
Towards direct reconstruction from a gamma camera based on compton scattering
Michael J Cree and Philip J Bones. Towards direct reconstruction from a gamma camera based on compton scattering. IEEE transactions on medical imaging , 13(2):398–407, 1994
work page 1994
-
[6]
Gamma-radiation imaging system based on the compton effect
DB Everett, JS Fleming, RW Todd, and JM Nightingale. Gamma-radiation imaging system based on the compton effect. In Proceedings of the Institution of Electrical Engineers , volume 124, pages 995–1000. IET, 1977. 5
work page 1977
-
[7]
Inversion formulas for the broken-ray radon transform
Lucia Florescu, Vadim A Markel, and John C Schotland. Inversion formulas for the broken-ray radon transform. Inverse Problems, 27(2):025002, 2011
work page 2011
-
[8]
Single-scattering optical tomography
Lucia Florescu, John C Schotland, and Vadim A Markel. Single-scattering optical tomography. Physical Review E, 79(3):036607, 2009
work page 2009
Show all 23 references
-
[9]
Spherical tomography and spherical integral geometry
S Gindikin, J Reeds, and L Shepp. Spherical tomography and spherical integral geometry. Tomography, Impedance Imaging, and Integral Geometry (South Hadley, MA, 1993), Lectures in Appl. Math , 30:83– 92, 1994
1993
-
[10]
Analytical reconstruction formula for n-dimensional conical radon transform
Rim Gouia-Zarrad. Analytical reconstruction formula for n-dimensional conical radon transform. Com- puters & Mathematics with Applications , 68(9):1016–1023, 2014
2014
-
[11]
Exact inversion of the conical radon transform with a fixed opening angle
Rim Gouia-Zarrad and Gaik Ambartsoumian. Exact inversion of the conical radon transform with a fixed opening angle. Inverse Problems, 30(4):045007, 2014
2014
-
[12]
Exact reconstruction formulas for a radon transform over cones
Markus Haltmeier. Exact reconstruction formulas for a radon transform over cones. Inverse Problems, 30(3):035001, 2014
2014
-
[13]
The radon transform, volume 2
Sigurdur Helgason and S Helgason. The radon transform, volume 2. Springer, 1999
1999
-
[14]
Inversion of a v-line transform arising in emission tomography
Yulia Hristova. Inversion of a v-line transform arising in emission tomography. Journal of Coupled Systems and Multiscale Dynamics , 3(3):272–277, 2015
2015
-
[15]
Inversion formulas for cone transforms arising in application of compton cameras
Chang-Yeol Jung and Sunghwan Moon. Inversion formulas for cone transforms arising in application of compton cameras. Inverse Problems, 31(1):015006, 2015
2015
-
[16]
Three-dimensional image reconstruction from compton camera data
Peter Kuchment and Fatma Terzioglu. Three-dimensional image reconstruction from compton camera data. SIAM Journal on Imaging Sciences , 9(4):1708–1725, 2016
2016
-
[17]
Analytic inversion of a conical radon transform arising in application of compton cameras on the cylinder
Sunghwan Moon and Markus Haltmeier. Analytic inversion of a conical radon transform arising in application of compton cameras on the cylinder. SIAM Journal on imaging sciences , 10(2):535–557, 2017
2017
-
[18]
On the v-line radon transform and its imaging applications
Marcela Morvidone, Ma¨ ı Khuong Nguyen, Tuong T Truong, and Habib Zaidi. On the v-line radon transform and its imaging applications. Journal of Biomedical Imaging , 2010:11, 2010
2010
-
[19]
Radon transforms on a class of cones with fixed axis direction
Mai K Nguyen, Tuong T Truong, and Pierre Grangeat. Radon transforms on a class of cones with fixed axis direction. Journal of Physics A: Mathematical and General , 38(37):8003, 2005
2005
-
[20]
The radon transform over cones with vertices on the sphere and orthogonal axes
Daniela Schiefeneder and Markus Haltmeier. The radon transform over cones with vertices on the sphere and orthogonal axes. SIAM Journal on Applied Mathematics , 77(4):1335–1351, 2017
2017
-
[21]
An electronically collimated gamma camera for single photon emission computed tomog- raphy
Manbir Singh. An electronically collimated gamma camera for single photon emission computed tomog- raphy. part i: Theoretical considerations and design criteria. Medical Physics, 10(4):421–427, 1983
1983
-
[22]
Some inversion formulas for the cone transform
Fatma Terzioglu. Some inversion formulas for the cone transform. Inverse Problems, 31(11):115010, 2015
2015
-
[23]
On new v-line radon transforms in R2 and their inversion
Tuong T Truong and Mai K Nguyen. On new v-line radon transforms in R2 and their inversion. Journal of Physics A: Mathematical and Theoretical , 44(7):075206, 2011. 6
2011
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.