![Arxiv:1810.09017V1 [Math.FA] 21 Oct 2018 Omadsmaitiitgasadivrinformulas](http://data.docslib.org/img/058b70f0401c7e8963d3dab992be890b-1.webp)
arXiv:1810.09017v1 [math.FA] 21 Oct 2018 omadSmaitiitgasadivrinformulas. inversion and integrals Semyanistyi and form om[,4 htge pfo h oko ikwk 1] nthe In trans- Funk [13]. the Minkowski is of example work typical the A from up 35]. grew [2, 27, that 18, n e.g., 4] 17, see, [3, to- 12, areas; form in 11, other arise 9, many int kind 7, and this their of geometry, from Problems convex sphere mography, submanifolds. Euclidean prescribed the over on grals functions reconstructing to peeisl.Adsiciefaueo uhtasom ncomparis in transforms if that such is of case feature Funk distinctive the with A itself. sphere rnfrsoe peia etosb yepae asn thro passing point hyperplanes fixed by arbitrary sections spherical over transforms sphere asn hog h origin the through passing oictoso hs rnfrsato ucin hc r defin are which functions on invest act their transforms for these of tools Modifications special need and transformations thogonal dmninlstig hstasomitgae ucino the on function a integrates transform this setting, -dimensional e od n phrases. and words Key 2010 peia nerlgoer oesawd ag fpolm relat problems of range wide a covers geometry integral Spherical ntepeetpprw bannwrslsrltdt h similar the to related results new obtain we paper present the In EOSRCINO UCIN NTHE ON FUNCTIONS OF RECONSTRUCTION ahmtc ujc Classification. Subject Mathematics r ls ominimal. to close are lnspsigtruhacmo point hype common by a sections through spherical to passing associated planes kinds two of transforms Abstract. rpi rjcinwt h oeat of hyper- pole families the the analytic to with are projection reduced tools graphic be main can The and λ transform. bowls, Radon or Transforms subspher plane truncated caps over transform. spherical integration Funk perform classical like kind the second subspheres to the complete reduced over of integration be by can defined and are kind first the n S PEEFO HI NERL OVER INTEGRALS THEIR FROM SPHERE csn rnfrs eynsisitgas n oie stereor- modified and integrals, Semyanisyi’s transforms, -cosine dmninlui peeo nteshr tef rnfrsof Transforms itself. sphere the on or sphere unit -dimensional n in R n +1 eoti e neso omlsfrteFn type Funk the for formulas inversion new obtain We YEPAESECTIONS HYPERPLANE vrcosscin fti peeb hyperplanes by sphere this of cross-sections over A ao rnfr n uktasomadcsn trans- cosine and transform Funk and transform Radon hc isete nieteshr ro the on or sphere the inside either lies which 1. o Introduction of A .RUBIN B. R 6= n o +1 1 hyd o omt with commute not do they , e 24 below. (2.4) see ; A rmr 41;Scnay44A15. Secondary 44A12; Primary supin o functions for Assumptions . A hc isisd the inside lies which es, r- igation. g an ugh all don ed unit or- on 5– ed e- 2 B. RUBIN truncated spheres, like spherical caps or bowls. An important problem is to reconstruct these functions from their integrals over the corre- sponding planar slices. Setting of the Problem, Brief History, and Main Results. n+1 Given a point A in R , let τA be a hyperplane passing through A. We assume A = (0,..., 0, a), 1 < a 1, and set − ≤ Sn = η =(η ,...,η ) Sn : η < a . a { 1 n+1 ∈ n+1 } Our main concern is the Radon-like transforms ( f)(τ )= f(η) d η, ( f)(τ )= f(η) d η, (1.1) Fa A τA Sa A τA n n τAZS τAZS ∩ ∩ a associated with spherical slices by the hyperplanes τA. These trans- forms coincide when a = 1, but have different features if a < 1. The case a = 0 for a gives the classical Funk transform. In the case, when a = 1, n = 2,F and f is smooth, these operators were introduced by Abouelaz and Daher [1] and studied by Helgason [7, p. 145]. The case a = 1 for all n 2 was studied in [27, Section 7.2], [28, Section 5] in the general context≥ of Lebesgue integrable functions. Operators a for 0 a < 1 and n = 2 were studied by Quellmalz [19], who suggestedF a nice≤ change of variable in the meridional direction and reduced the inversion problem for af to the similar one for the Funk transform. The function f wasF assumed either continuous or 2 2 1 belonging to L (S ). The injectivity of a with a = 1 was obtained in [19] as a limiting case under certain smoothnessF assumption for f. Another point of view on the inversion problem for a and more general operators, related to planar sections of S2, was suggestedF by Gindikin et al. [6] in terms of the kappa-operator. The integrals af were studied by Salman [30, 31] when 0 a< 1, n 2, and f isS a continuously differentiable function with≤ compact ≥ Sn support strictly inside a . He used the stereorgraphic projection to reduce inversion of af to the similar problem for a certain Radon-like transform over spheresS in Rn. The latter was explicitly inverted by making use some technique borrowed from John [10]. The case a = 1 was covered in [30, 31] by passing to the limit as a 1. The aim of the present paper is twofold. Firstly,→ we extend Quell- malz’s approach from [19] to all n 2 and all a ( 1, 1), and per- form this extension under minimal assumptions≥ for∈f.− Specifically, we
1After the manuscript had been completed, the author became aware of the recent preprint by Quellmalz [20], where the results of [19] were extended to all n 2. The method of [20] essentially differs from ours. ≥ RECONSTRUCTION OF FUNCTIONS ON THE SPHERE 3 assume f L1(Sn), as in the Funk case. To circumvent the main tech- ∈ nical difficulty related to the change of variables, we consider a as a member of the analytic family of the relevant λ-cosine transformsF for which the corresponding Jacobians can be easily computed. Secondly, we suggest a new approach to the study of the operators a. Instead of the stereographic projection with the pole (0,..., 0, 1), Sas in [31], we choose the pole at (0,..., 0, a). The case a = 1 agrees with the usual stereographic projection, as in [7] for n = 2. Asa result, a expresses through the classical hyperplane Radon transform Rf of functionsS on Rn and all known properties of the latter, including inversion formulas, transfer to a. To perform this transition, we invoke analytic families of the λ-cosineS transforms and Semyanistyi’s fractional integrals [33], [27, Section 4.9]. These remarkable integrals include the Radon transform as a particular case. As an interesting by-product, we obtain a formula Cλ = Rλ , (1.2) a Va,λ Ua,λ Cλ establishing connection between the λ-cosine transform a and the Se- λ myanistyi transform R via certain bijections a,λ and a,λ (see Lemma 4.72). Similar formula was known before (see,V e.g., [27,U Section 5.2]) only for λ = 1 and a = 0, that corresponds to the Funk and Radon transforms. The− limiting case λ = 1 in (1.2) gives the desired con- − nection between a and R. It is worth notingS that the λ-cosine transforms, Semyanistyi’s inte- grals, and their modifications have proved to be important objects in integral geometry and harmonic analysis (see, e.g., [15, 16, 21–27]) be- cause they enable one to obtain diverse results for particular values of λ by making use of analytic continuation. We believe that (1.2) paves the way to new investigations in the area. For possible applications and further purposes, we separately men- tion particular cases a = 0, a = 1, and the case of zonal functions, when all formulas look much simpler. The paper is organized as follows. Section 2 contains necessary pre- liminaries. Section 3 and 4 are devoted to the operators a and a, respectively. The Appendix contains technical results relatedF to limSit- ing properties of some integrals depending on a parameter.
2. Preliminaries 2.1. Notation. In the following Sn, n 2, denotes the unit sphere n+1 ≥ in the real Euclidean space R ; e1,...,en+1 are the coordinate unit n vectors; the space R = Re1 Ren is considered as a coordinate hyperplane of Rn+1. For η =⊕···⊕ (η ,...,η ) Sn, dη stands for the 1 n+1 ∈ 4 B. RUBIN
Sn (n+1)/2 Riemannian measure on ; σn Sn dη =2π Γ((n +1)/2). We ≡ n write d η = dη/σn for the normalized surface element on S ; η⊥ is the ∗ R linear subspace orthogonal to η. The notation
p = (0,..., 0, 1), p = (0,..., 0, 1) N S − will be used for the North Pole and the South Pole of Sn, respectively. Sn Sn Sn Sn We set + = η : ηn+1 0 , = η : ηn+1 < 0 , { ∈ ≥ } − { ∈ } A = (0,..., 0, a), 1 < a 1. − ≤ p n n n The notation for the function spaces L (S ), C(S ), C∞(S ) and the Rn Rn Rn similar spaces on is standard. We write Cc( ) and Cc∞( ) for the corresponding spaces of compactly supported functions. Dealing with integrals, we say that the integral exists in the Lebesgue sense if it is finite when the expression under the sign of integration is replaced by its absolute value. The letter c (sometimes with sub- scripts) stands for an inessential positive constant that may vary at each occurrence.
2.2. Riesz potentials. The Riesz potential of a function g on R is defined by the formula
α 1/2 α 1 g(y) dy 2 π Γ(α/2) (I g)(x)= 1 α , γ1(α)= . (2.1) γ1(α) x y − Γ((1 α)/2) ZR | − | − We will need this operator only with 0 <α< 1, though other values of α can be considered. The following fact is a particular case of the more general statement in [27, Lemma 3.2].
Lemma 2.1. Let g be a continuous function in some neighborhood Ω of x R, and let ∈ yg(y) dy < . (2.2) | | ∞ RnZ Ω \ Then lim (Iαg)(x)= g(x). (2.3) α 0 → A generalization of this statement to the case when the function under the sign of integration depends also on α (see Lemma 5.2 in Appendix) will be used in the present paper. RECONSTRUCTION OF FUNCTIONS ON THE SPHERE 5
2.3. The Funk Transform and the Spherical Slice Transform. The Funk transform of a continuous function f on Sn has the form
1 n ( f)(ξ)= f(η) dξη, ξ S , (2.4) F σn 1 ∈ − SnZξ⊥ ∩ n where dξη stands for the Riemannian surface measure on S ξ⊥. If f L1(Sn), then ( f)(ξ) is finite for almost all ξ Sn and the∩ operator ∈ extends as a linearF bounded operator from L1(S∈n) to L1(Sn). PropertiesF of the Funk transform are described in numerous publications; see, e.g., [5, 7, 17, 27]. In particular, if η = (η ,...,η ) Sn and f(η) 1 n+1 ∈ ≡ f0(ηn+1), i.e. f is a zonal function, then f is zonal too. If, for f F 2 zonal, we set ( f)(ξ)= F0(ρ), where ρ = 1 ξn+1 is the sine of the geodesic distanceF between ξ and the North Pole− p , then p N ρ 2 n 2 2 (n 3)/2 2σn 2 0(ρ)= c ρ − f0(t)(ρ t ) − dt, c = − ; (2.5) F − σn 1 Z0 − cf. [27, Theorem 5.21]. Abel type integrals of the form (2.5) were studied in many sources; see, e.g., [27, 32]. We shall consider modifications of (2.4) for the non-central cross- sections of Sn by the hyperplanes passing through an arbitrary fixed point inside the sphere Sn or on the sphere itself. Without loss of generality, we choose such a point to be A = (0,..., 0, a), 1 < a 1. − ≤ Replacing ξ⊥ in (2.4) by the shifted hyperplane ξ⊥ + aen+1, we obtain the collection of spherical slices γ (ξ) = η Sn : ξ (η ae )=0 a { ∈ · − n+1 } = η Sn : ξ η = aξ , (2.6) { ∈ · n+1} labeled by the point ξ. The corresponding spherical slice transform is defined by ( f)(ξ)= f(η) d η, ξ Sn, (2.7) Fa γa ∈ γaZ(ξ) where dγa η stands for the surface area measure on γa(ξ). Clearly, γ ( ξ) = γ (ξ) and ( f)( ξ)=( f)(ξ). Every slice γ (ξ) is an a − a Fa − Fa a (n 1)-dimensional Euclidean sphere with center at (aξn+1)ξ and ra- dius− (1 a2ξ2 )1/2. Hence − n+1 2 2 (n 1)/2 ( a1)(ξ)= σn 1(1 a ξn+1) − . (2.8) F − − This simple equality implies the following statement. 6 B. RUBIN
Proposition 2.2. Let 1 < a 1, n 2. The slice transform a is − ≤ n ≥ F a linear bounded operator on L∞(S ). If 1 1 2 (Mξf)(t) = f(√1 t η + tξ) dξη (2.9) σn 1 − − SnZξ⊥ ∩ 2 (1 n)/2 (1 t ) − = − f(η) dξ,tη, 1 1 2 (n 2)/2 h(ξ η) f(η) dη = σn 1 h(t)(Mξf)(t)(1 t ) − dt (2.12) · − − SZn Z1 − provided that at least one of these integrals exists in the Lebesgue sense. In the case h 1 (2.12) has a simpler form ≡ 1 2 (n 2)/2 2 f(η) dη = (1 t ) − dt f √1 t ψ + te dψ; (2.13) − − n+1 SZn Z1 SnZ−1 − see [27, Lemma 1.34]. 2.4. The λ-Cosine Transform. In many occurrences it is convenient to consider the Funk transform (2.4) as a member of the meromorphic family of normalized cosine (or λ-cosine) transforms λ λ ( f)(ξ)= γn,λ f(η) ξ η d η, (2.14) C | · | ∗ SZn where d η = dη/σn, Reλ > 1, λ =0, 2, 4,..., ∗ − 6 π1/2 Γ( λ/2) γ = − . (2.15) n,λ Γ((n+1)/2) Γ((1+λ)/2) RECONSTRUCTION OF FUNCTIONS ON THE SPHERE 7 This transform has a long history; see [27, pp. 283, 361, 527], [15] and references therein. We restrict to the case 1 <λ< 1, which is sufficient for our purposes. − Lemma 2.3. (cf. [27, p. 531]) If f is a continuous function on Sn, then for every ξ Sn, ∈ 1/2 λ π lim ( f)(ξ)= cn ( f)(ξ), cn = . (2.16) λ 1 →− C F Γ(n/2) Both transforms (2.4) and (2.14) annihilate odd functions that form the kernel (the null space) of the corresponding operators. We introduce the “shifted” cosine transform ( λf)(ξ)= γ f(η) ξ (η ae ) λ dη, 1 < a 1, (2.17) Ca n,λ | · − n+1 | − ≤ SZn where λ and γn,λ are the same as in (2.14). Lemma 2.4. If f is a continuous function on Sn, then for every ξ Sn and aξ < 1, ∈ | n+1| λ dn π lim ( a f)(ξ)= ( af)(ξ), dn = . (2.18) λ 1 C 1 a2ξ2 F Γ((n+1)/2) →− − n+1 Proof. By (2.12), p 1 λ λ 2 (n 2)/2 ( a f)(ξ)= γn,λσn 1 t aξn+1 (Mξf)(t)(1 t ) − dt. C − | − | − Z1 − For every fixed ξ, this integral is a constant multiple of the Riesz po- tential 1 α 1/2 α 1 g(t) dt 2 π Γ(α/2) (I g)(x)= 1 α , γ1(α)= (2.19) γ1(α) x t − Γ((1 α)/2) Z1 | − | − − (cf. (5.7)) with 2 (n 2)/2 α = λ +1, x = aξ , g(t)=(M f)(t)(1 t ) − . n+1 ξ − Hence, by Lemma 2.1 and (2.11), a simple calculation gives λ 2 2 (n 2)/2 lim ( a f)(ξ) = dnσn 1(Mξf)(aξn+1)(1 a ξn+1) − λ 1 − →− C − dn = ( af)(ξ). 1 a2ξ2 F − n+1 p 8 B. RUBIN 3. Spherical Slice Transforms via the Funk Transform In this section we extend Quellmalz’s method for the slice transform a (see [19] for n = 2) to all n 2. Everywhere in this section we Fassume 1 1 u a µ− :η ˜(ψ, v) η(ψ,u)=˜η ψ, − (3.4) → 1 au − and moves the “equator” v = 0 to the “parallel” u = a while keeping the poles pN and pS fixed. For the sake of convenience, we abuse the notation and write f(ψ,u) in place of f(η(ψ,u)). We will need one more bijective map Sn ξ(ϕ,s) ν ξ˜(ϕ, t) Sn, (3.5) ∋ −→ ∈ where √1 a2 t = s − . (3.6) √1 a2s2 − Then 2 1 √1 a ν− : ξ˜(ϕ, t) ξ(ϕ,s)= ξ˜ ϕ,s − . (3.7) → √1 a2s2 − Lemma 3.1. If f and Φ belong to L1(Sn), then ∈ 1 2 n/2 (f µ− )(˜η) f(η) dη = (1 a ) ◦ n dη,˜ (3.8) − (1 + aη˜n+1) SZn SZn (Φ ν)(ξ) ˜ ˜ √ 2 Φ(ξ) dξ = 1 a 2 ◦2 (n+1)/2 dξ. (3.9) − (1 a ξn+1) SZn SZn − RECONSTRUCTION OF FUNCTIONS ON THE SPHERE 9 Proof. By (2.13), 1 2 (n 2)/2 2 f(η) dη = (1 u ) − du f(√1 u ψ + ue ) dψ. − − n+1 SZn Z1 SnZ−1 − The result follows if we change variables according to (3.3), taking into account that du 1 a2 (1 a2)(1 v2) = − , 1 u2 = − − . dv (1 + av)2 − (1 + av)2 The proof of the second statement is similar. Our next aim is to convert the slice transform (2.7) into the Funk transform (2.4) by making use of the map (3.2). To circumvent tech- nicalities related to Jacobians, we invoke Lemma 2.4 for the shifted cosine transform. Given 1 Proposition 3.2. The maps a,λ, a,λ, and their inverses preserve n M N n p n the space C∞(S ) and are linear bounded operators on C(S ) and L (S ), 1 p . ≤ ≤∞ Proof. The statements become obvious if we observe that 2 1+ aη˜ 1 a > 0, 1 1 a2ξ2 1 a2 > 0, ≥ n+1 ≥ −| | ≥ − n+1 ≥ − because 0 a<1. The result for Lp functions follows from Lemma 3.1. ≤ λ λ Lemma 3.3. Let and a be the λ-cosine transforms (2.14) and (2.17), respectively.C If f LC1(Sn), then ∈ λf = λ f. (3.12) Ca Na,λC Ma,λ λ Sn Moreover, a preserves the space C∞( ) and is a linear bounded op- erator on CC(Sn) and Lp(Sn), 1 p . ≤ ≤∞ 10 B. RUBIN Proof. We transform the left-hand side by making use of Lemma 3.1. Changing variables, for ξ = ξ(ϕ,s) and η = η(ψ,u) we have ξ (η ae )= √1 s2 √1 u2 (ϕ ψ)+ s(u a) · − n+1 − − · − √1 s2 √1 a2 √1 v2 sv(1 a2) = − − − (ϕ ψ)+ − 1+ av · 1+ av √1 a2 √1 a2s2 √1 v2 √1 s2 sv√1 a2 = − − − − (ϕ ψ)+ − 1+ av √1 a2s2 · √1 a2s2 − − √1 a2 √1 a2s2 s2(1 a2) s√1 a2 = − − √1 v2 1 − (ϕ ψ)+ v − 1+ av − − 1 a2s2 · √1 a2s2 " r − − # √1 a2 √1 a2s2 √1 a2 1 a2ξ2 = − − (ξ˜ η˜)= − − n+1 (ξ˜ η˜). 1+ av · 1+paη˜n+1 · Hence, by Lemma 3.1, f(η) ξ (η ae ) λ dη | · − n+1 | SZn 1 2 (n+λ)/2 2 2 λ/2 (f µ− )(˜η) ˜ λ = (1 a ) (1 a ξn+1) ◦ n+λ ξ η˜ dη.˜ − − (1 + aη˜n+1) | · | SZn This gives (3.12). The last statement follows from the similar facts for λ [27, Section 5.1] in view of Proposition 3.2. C Our next aim is to pass to the limit (3.12) as λ 1. This will give → − us the desired relation between the slice transform a and the Funk transform . We denote F F 1 (f µ− )(˜η) af)(˜η) = ◦ n 1 ( a,λf)(˜η) λ= 1, (3.13) M (1 + aη˜n+1) − ≡ M − 2 (n 1)/2 ( aΦ)(ξ) = σn 1 (1 a ) − (Φ ν)(ξ) (3.14) N − − ◦ 2 2 σn 1 1 a ξn+1 ( a,λΦ)(ξ) λ= 1. ≡ − − N − q Theorem 3.4. If f L1(Sn), then ∈ f = f. (3.15) Fa NaFMa Proof. In view of Proposition 3.2 and the L1-boundedness of the Funk transform, it suffices to prove (3.15) for f C(Sn). We write (3.12) as ∈ ( λf)(ξ)= h (ξ)( λg )(ξ˜), where Ca λ C λ 1 2 (n+λ)/2 2 2 λ/2 (f µ− )(˜η) hλ(ξ)=(1 a ) (1 a ξn+1) , gλ(˜η)= ◦ n+λ . − − (1 + aη˜n+1) RECONSTRUCTION OF FUNCTIONS ON THE SPHERE 11 Hereη ˜ and ξ˜ are defined by (3.2) and (3.5), respectively. Note that n gλ(˜η) is a continuous function of (λ, η˜) on [ 1, 0] S (we recall that 1+ aη˜ > 0 because 0 a< 1). If λ −1, then× n+1 ≤ → − 2 (n 1)/2 1 (1 a ) − (f µ− )(˜η) hλ(ξ) − , gλ(˜η) ◦ =( af)(˜η). 2 2 n 1 → 1 a ξ → (1 + aη˜n+1) − M − n+1 Hence, by Lemmap 5.3 in Appendix, 1/2 λ π lim ( gλ)(ξ˜)= cn ( af)(ξ˜), cn = , λ 1 Γ(n/2) →− C FM where F is the Funk transform (2.4). Further, by Lemma 2.4, λ dn π lim ( a f)(ξ)= ( af)(ξ), dn = . λ 1 C 1 a2ξ2 F Γ((n+1)/2) →− − n+1 Finally, (3.12) yields p 2 (n 1)/2 ( af)(ξ)= σn 1 (1 a ) − ( af)(ξ˜)=( a af)(ξ). F − − FM N FM Because the Funk transform is well-investigated [5, 7, 17, 27], Theo- rem 3.4 enables us to establish the corresponding properties of a. In particular, Proposition 3.2 yields the following corollary. F n Corollary 3.5. The slice transform a preserves the space C∞(S ) and is a linear bounded operator on C(FSn) and Lp(Sn), 1 p . ≤ ≤∞ It is known [7, 27] that the Funk transform , acting on L1(Sn) or C(Sn), is injective on even functions. Odd functionsF form the kernel of . Because the operators a and a in (3.15) are bijective, it follows Fthat is injective on theM class of functionsN f satisfying Fa ( f)(˜η)=( f)( η˜) (3.16) Ma Ma − for all (or almost all)η ˜ Sn. Similarly, the functions f for which ∈ ( f)(˜η)= ( f)( η˜) (3.17) Ma − Ma − form the kernel of Fa. Below we characterize (3.16) and (3.17) in terms of f itself. Let, as (3.1) and (3.2), η η(ψ,u),η ˜ = µ(η) η˜(ψ, v), where ≡ ≡ η(ψ,u)= √1 u2 ψ + ue , η˜(ψ, v)= √1 v2 ψ + ve , − n+1 − n+1 v + a u a u = , v = − . (3.18) 1+ av 1 ua − 12 B. RUBIN By (3.13) 1 1 2 − √ ( af)(˜η) = n 1 (f µ )( 1 v ψ + ven+1), M (1 + av) − ◦ − 1 1 2 ( f)( η˜) = (f µ− )( √1 v ψ ve ). Ma − (1 av)n 1 ◦ − − − n+1 − − Using (3.18), it is straightforward to show that (3.16) and (3.17) can be equivalently written as f(η)= ρ(η ) f(η ) (or f(η )= ρ(η) f(η)) (3.19) ∗ ∗ ∗ and f(η)= ρ(η ) f(η ) (or f(η )= ρ(η) f(η)), (3.20) − ∗ ∗ ∗ − respectively. Here ρ(η) ρ(ψ,u), ρ(η ) ρ(ψ ,u ), and η η (ψ ,u ) have the following definition≡ in spherical∗ ≡ coordinates:∗ ∗ ∗ ≡ ∗ ∗ ∗ n 1 n 1 1+ a2 2ua − 1+ a2 2u a − ρ(η)= − , ρ(η )= − ∗ , 1 a2 ∗ 1 a2 − − 2 2 2a u(1 + a ) η (ψ ,u )= 1 u ψ + u en+1, ψ = ψ, u = − . ∗ ∗ ∗ − ∗ ∗ ∗ ∗ − ∗ 1+ a2 2ua − One can showp (we leave this calculation to the reader) that the points η, η , and A = (0,..., 0, a) are located on the same line and the corre- sponding∗ reflection map Sn η η Sn is an involution. The weighted equality (3.19)∋ can→ be∗ ∈ regarded as a certain symmetry condition which is necessary and sufficient for injectivity of a. In the Funk transform case a = 0, (3.19) boils down to the classicalF evenness condition f(η)= f( η). The above reasoning− can be summarized as follows. 1 n n Theorem 3.6. The slice transform a, acting on L (S ) or C(S ), is injective if and only if its action is restrictedF to functions f satisfying (3.19). The functions f satisfying (3.20) form the kernel of . Fa Theorem 3.7. An integrable function f satisfying the symmetry condi- tion (3.19) can be recovered from the slice transform f by the formula Fa 1 1 1 f = − − − f (3.21) Ma F Na Fa 1 where − is the inverse Funk transform, F 1 ˜ 1 ˜ ( a− f)(η) ( a− f)(ψ,u) M ≡n 1M 1 a2 − √1 a2 √1 u2 ψ +(u a)e = − f˜ − − − n+1 , (3.22) 1 ua 1 ua − − RECONSTRUCTION OF FUNCTIONS ON THE SPHERE 13 1 1 ( − Φ)(˜ ξ˜) ( − Φ)(˜ ϕ, t) Na ≡ Na 2 (1 n)/2 2 2 (1 a ) − √1 a √1 t ϕ+te = − Φ˜ − − n+1 . (3.23) σn 1 √1 a2 +a2t2 − − Proof. The formula (3.22) follows by routine calculation from (3.1)- (3.4), (3.13). Similarly, (3.23) is a consequence of (3.5)-(3.7) and (3.14). Remark 3.8. The class of all functions f that fall into the scope of 1 1 Sn Theorem 3.7 is non-empty and has the form a− (Leven( )), where 1 M 1 1 Sn the operator a− is defined by (3.22). Similarly, a− (Lodd( )) is the kernel of Min L1(Sn). M Fa 1 Remark 3.9. Numerous formulas for the inverse Funk transform − in smooth and nonsmooth settings can be found in the literature;F see, e.g., [27, Section 5.1] and references therein. For example, the following statement in terms of spherical means (2.9) is known (cf. [27, Theorem 5.40]). p Sn 2 Theorem 3.10. Let f Leven( ), 1 p< , Φξ(s)=(Mξ f)(√1 s ). Then ∈ ≤ ∞ F − t n 1 1 ∂ − 2 2 2 (n 3)/2 n 1 f(ξ)=lim (t s ) − Φξ(s) s − ds . t 1 → 2t ∂t (n 2)! − − Z0 (3.24) In particular, for n odd, 1/2 (n 1)/2 π 1 ∂ − n 2 f(ξ)=lim [t − Φξ(t)]. (3.25) t 1 Γ(n/2) 2t ∂t → The limit in these formulas is understood in the Lp-norm. If f n ∈ Ceven(S ), it can be interpreted in the sup-norm. 4. Spherical Slice Transforms via the Radon Transform In this section the setting of the problem is slightly changed and the tools are different. Given 1 < a 1, let − ≤ Sn Sn Sn Sn Sn a = η : ηn+1 < a , 0 = η : ηn+1 < 0 . (4.1) { ∈ } − ≡ { ∈ } We consider a modified spherical slice transform that takes a function f on Sn to a function f on the set of all spherical slices a Sa γ (ξ) = η Sn : ξ (η ae )=0 a { ∈ a · − n+1 } = η Sn : ξ η = aξ , (4.2) { ∈ a · n+1} 14 B. RUBIN labeled by ξ Sn . The value of f at γ (ξ) is defined by the integral ∈ + Sa a ( f)(ξ)= f(η) d η, (4.3) Sa γa γaZ(ξ) where dγa η is the surface measure on γa(ξ). In the cases a = 0 and a = 1, we write af as 0f and 1f, respectively. The slices (4.2) are truncated versionsS of thoseS in theS previous section. They resemble multidimensional spherical “arcs” that become complete spheres only if a =1. If ξ Sn , then ∈ + ( f)(ξ)=( f˜)(ξ), (4.4) Sa Fa where f˜(η)= f(η) if ηn+1 < a and f˜(η) = 0, otherwise. Proposition 4.1. If 1 f p Sn = f˜ p Sn c f˜ p Sn = c f p Sn . ||Sa ||L ( +) ||Fa ||L ( ) ≤ || ||L ( ) || ||L ( a ) 4.1. Stereographic Projection for Truncated Sphere. We recall Rn that the space is identified with the coordinate hyperplane en⊥+1 in Rn+1. Consider the bijective map Rn x π η Sn, (4.5) ∋ −→ ∈ a where the point η = π(x) is located on the line ℓ passing through the point A = (0,..., 0, a) and the point X = x en+1 lies in the tangent n − plane to S at the South Pole pS = (0,..., 0, 1). To find η = π(x) explicitly, we make use of the parametric representation− of ℓ, so that ℓ = y Rn+1 : y = tA + (1 t)X; t R , { ∈ − ∈ } and find the value of t for which y Sn. Let ∈ a n 1 η = ψ sin ω + e cos ω, ψ S − , 0 ω π. n+1 ∈ ≤ ≤ n+1 Then, for y =(y′,y ) R we have n+1 ∈ ℓ = y : y′ = (1 t)x, y = ta + t 1; t R . { − n+1 − ∈ } The condition y = 1 with y Sn yields the quadratic equation | | ∈ a (1 t)2 x 2 +(ta + t 1)2 =1, t< 0. − | | − The solution to this equation has the form x 2 + a +1 D(x) t = | | − , x 2 +(a + 1)2 | | RECONSTRUCTION OF FUNCTIONS ON THE SPHERE 15 where D(x)= √a +1 x 2(1 a)+ a +1. (4.6) | | − This gives p η π(x)= P (x) x + Q(x) e Sn, (4.7) ≡ n+1 ∈ a where a(a +1)+ D(x) a x 2 (a + 1)D(x) P (x)= , Q(x)= | | − . (4.8) x 2 +(a + 1)2 x 2 +(a + 1)2 | | | | In particular, 1 for a = 0: D(x)= x 2 +1, P (x)= ; (4.9) | | x 2 +1 | | p 4 for a = 1: D(x)=2, P (x)= p. (4.10) x 2 +4 | | The inverse map −1 p Sn η π x Rn, (4.11) a ∋ −−−→ ∈ has much simpler analytic expression that follows from elementary ge- ometrical consideration. Specifically, if η =(η′, ηn+1), η′ =(η1,...,ηn), n 1 or η = ψ sin ω + e cos ω, ψ S − , r = x , then n+1 ∈ | | (a + 1) 1 η2 (a + 1) sin ω x = − n+1 = , (4.12) | | a η a cos ω −p n+1 − and therefore 1 (a + 1) η′ x = π− (η)= . (4.13) a η − n+1 In particular, (4.7) and (4.12) yield 1 a ηn+1 1 P (π− (η)) = − , Q(π− (η)) = η . (4.14) a +1 n+1 Remark 4.2. In the cases a = 0 and a = 1, of slices through the origin and the North Pole pN , we have the following simple formulas: For a = 0: x e x e η = − n+1 = − n+1 ; x Rn, (4.15) x 2 +1 x en+1 ∈ | | | − | η p x = ′ , η Sn . (4.16) η ∈ − | n+1| For a = 1: 4x +( x 2 4)e η = | | − n+1 , x Rn; (4.17) x 2 +4 ∈ | | 2η x = ′ , η Sn p . (4.18) 1 η ∈ \ N − n+1 16 B. RUBIN Lemma 4.3. Let 1 < a 1, − ≤ (a + 1)n(1 aη ) w(η)= − n+1 , η Sn, (4.19) (a η )n+1 ∈ a − n+1 so that D(x) (w π)(x)= , x Rn, (4.20) ◦ (a + 1) P n(x) ∈ where D(x) and P (x) are defined by (4.6) and (4.8), respectively. Then for any g L1(Rn) and f L1(Sn), ∈ ∈ a 1 g(x) dx = (g π− )(η) w(η) dη, (4.21) ◦ n n RZ SZa dx f(η) dη = (f π)(x) . (4.22) ◦ (w π)(x) n n SZa RZ ◦ Proof. We observe that (4.20) follows from (4.19) by virtue of (4.14). Further, passing to polar coordinates and using (4.12), we obtain ∞ n 1 g(x) dx = r − dr g(rθ) dθ RZn Z0 SnZ−1 π n 1 (a+1) sin ω − (a+1)(1 cos ω) = − dω a cos ω (a cos ω)2 cosZ−1 a − − (a+1) sin ω g θ dθ, × a cos ω SnZ−1 − which is the right-hand side of (4.21) in spherical coordinates. The 1 equality (4.22) follows from (4.21) if we set f(η)=(g π− )(η) w(η). ◦ Example 4.4. In the simplest cases, a = 0 and a = 1, the formulas (4.21) and (4.22) have the following form. For a = 0: η′ dη g(x) dx = g n+1 , (4.23) ηn+1 ηn+1 n n RZ SZ− | | | | x e dx f(η) dη = f − n+1 . (4.24) x 2 +1 ( x 2 + 1)(n+1)/2 n n ! SZ− RZ | | | | For a = 1: p RECONSTRUCTION OF FUNCTIONS ON THE SPHERE 17 n 2η′ dη g(x) dx =2 g n , (4.25) 1 ηn+1 (1 ηn+1) RZn SZn − − 4x +( x 2 4)e dx f(η) dη =4n f | | − n+1 . (4.26) x 2 +4 ( x 2 + 4)n SZn RZn | | | | 4.2. Semyanistyi’s Integrals, Cosine, and Radon Transforms. 1. The map (4.5) extends to the corresponding relationship between affine hyperplanes in Rn and spherical slices (4.2). Specifically, to every n (n 1)-plane τ in R = e⊥ we associate a spherical slice in the n- − n+1 dimensional plane spanned by the shifted (n 1)-plane τ en+1 and the point A = (0,..., 0, a). This relationship gives− rise to the− correspond- ing connection between the hyperplane Radon transform of functions on Rn and the spherical slice transform (4.3). Before we describe this connection analytically, it is instructive to recall some basic facts from related to the Radon transform and the corresponding analytic family of Semyanistyi’s integrals; see, e.g., [7, 27]. The hyperplane Radon transform R takes a sufficiently good contin- uous function g on Rn to integrals of g over affine hyperplanes in Rn. Every such hyperplane τ can be parameterized as n n 1 τ τ(θ, t)= x R : x θ = t , θ S − , t R, ≡ { ∈ · } ∈ ∈ so that τ(θ, t)= τ( θ, t). In this notation, − − (Rg)(τ) (Rg)(θ, t)= g(tθ + u) d u, (4.27) ≡ θ θZ⊥ where dθu denotes the Euclidean measure on θ⊥. Rn If Cc∞( ) is the space of compactly supported infinitely differen- Rn Sn 1 R tiable functions on and Cc∞( − ) is the similar space on Sn 1 R Rn × Sn 1 R − , then R acts from Cc∞( ) to Cc∞( − ). A similar × n n 1 × statement holds for the spaces Cc(R ) and Cc(S − R) of compactly supported continuous functions. × n 1 The integral (4.27) is well-defined for almost all (θ, t) S − R whenever g Lp(Rn), 1 p (Rλg)(θ, t)= γ g(x) x θ t λ dx, Reλ> 1, (4.30) n,λ | · − | − RZn where γn,λ is the constant (2.15). If g is good enough, then λ π lim (R g)(θ, t)= dn(Rg)(θ, t), dn = . (4.31) λ 1 →− Γ((n+1)/2) This equality (in different notation) is well known (see, e.g., [27, Section 4.9]). It will be proved in Appendix (see Lemma 5.4) for more general integrals, when g depends also on λ. 2. We introduce an analytic family of cosine transforms associated Sn with slices in a by the formula (Cλf)(ξ)= γ f(η) ξ η aξ λ dη, ξ Sn , (4.32) a n,λ | · − n+1| ∈ + n SZa Rn 2 Sn (cf. (2.17)). For x and ξ ξ(θ,s)= √1 s θ+sen+1 + pN , n 1 ∈ ≡ − ∈ \ θ S − , s [0, 1), we set ∈ ∈ (f π)(x) ( f)(x) = ◦ P λ(x), (4.33) Ua,λ (w π)(x) ◦ 2 λ/2 (a + 1)s ( a,λΦ)(ξ) ( a,λΦ)(θ,s)=(1 s ) Φ θ, , (4.34) V ≡ V − √1 s2 − where P (x) is defined by (4.8), and (w π)(x) has the form (4.20). ◦ RECONSTRUCTION OF FUNCTIONS ON THE SPHERE 19 Example 4.5. In the cases a = 0 and a = 1, (4.33) and (4.34) have the following form (cf. Remark 4.2). For a = 0: 1 x en+1 ( 0,λf)(x) = 2 (λ+n+1)/2 f − , (4.35) U ( x + 1) x 2 +1! | | | | 2 λ/p2 s ( 0,λΦ)(ξ) ( 0,λΦ)(θ,s)=(1 s ) Φ θ, . (4.36) V ≡ V − √1 s2 − For a = 1: 4 n+λ 4x +( x 2 4)e ( f)(x) = f | | − n+1 , (4.37) U1,λ x 2 +4 x 2 +4 | | | | 2 λ/2 2s ( 1,λΦ)(ξ) ( 1,λΦ)(θ,s)=(1 s ) Φ θ, . (4.38) V ≡ V − √1 s2 − Sn Sn Let Cc( a ) be the space of continuous functions on having com- Sn pact support in a . Such functions vanish identically near the parallel ηn+1 = a. The corresponding space of infinitely differentiable functions Sn Sn Sn is denoted by Cc∞( a ). The notations Cc( + pN ) and Cc∞( + pN ) \ Sn \ will be used for the similar spaces of functions on + pN . The definitions (4.33) and (4.34) imply the following\ statements. Proposition 4.6. Sn Sn (i) The operator a,λ acts bijectively from Cc( a ) (from Cc∞( a )) to Rn RUn 1 Cc( ) (to Cc∞( )). The inverse operator a,λ− acts similarly in the opposite direction. U Sn 1 R Sn 1 (ii) The operator a,λ acts bijectively from Cc( − )(from Cc∞( − R Sn V Sn × 1 × )) to Cc( + pN ) (to Cc∞( + pN )). The inverse operator a,λ− acts similarly in the\ opposite direction.\ V Lemma 4.7. Let f C (Sn), 1 < a 1, ∈ c a − ≤ 2 n n 1 ξ = √1 s θ + se S p , θ S − , s [0, 1). − n+1 ∈ + \ N ∈ ∈ Then (Cλf)(ξ) (Cλf)(θ,s)=( Rλ f)(θ,s), (4.39) a ≡ a Va,λ Ua,λ Proof. By (4.22), dx (Cλf)(ξ)= γ (f π)(x) ξ π(x) aξ λ . a n,λ ◦ | · − n+1| (w π)(x) RZn ◦ 20 B. RUBIN Further, by (4.7), ξ π(x) aξ = ξ (π(x) ae ) · − n+1 · − n+1 =(√1 s2θ + se ) [P (x) x +(Q(x) a)e ] − n+1 · − n+1 (Q(x) a)s = P (x)√1 s2 x θ + − = P (x) √1 s2 (x θ t), − · P (x) √1 s2 − · − − where (a Q(x))s (a + 1)s t = − = . (4.40) P (x) √1 s2 √1 s2 − − The last expression follows from (4.8). Thus, dx (Cλf)(ξ) = γ (1 s2)λ/2 (f π)(x) P (x)(x θ t) λ a n,λ − ◦ | · − | (w π)(x) RZn ◦ = γ (1 s2)λ/2 ( f)(x) x θ t) λ dx, n,λ − Ua,λ | · − | RZn as desired; cf. (4.34). Our next aim is to pass to the limit in (4.39) as λ 1. For x Rn and ξ ξ(θ,s) Sn , we denote → − ∈ ≡ ∈ + n 1 (f π)(x) (a + 1) P − (x) ( f)(x) = ◦ = (f π)(x), (4.41) Ua (w π)(x) P (x) D(x) ◦ ◦ 1 a2s2 (a + 1)s ( aΦ)(ξ) = − Φ θ, , (4.42) V 1 s2 √1 s2 r − − where P (x) and D(x) are defined by (4.8) and (4.6), respectively, and (w π)(x) has the form (4.20). In particular, ◦ for a = 0: 1 x en+1 ( 0f)(x) = 2 n/2 f − , (4.43) U ( x + 1) x 2 +1! | | | | 1 ps ( 0Φ)(ξ) = Φ θ, ; (4.44) V 1 s2 √1 s2 r − − for a = 1: n 1 4 − 4x +( x 2 4)e ( f)(x) = f | | − n+1 , (4.45) U1 x 2 +4 x 2 +4 | | | | 2s ( 1Φ)(ξ) = Φ θ, . (4.46) V √1 s2 − RECONSTRUCTION OF FUNCTIONS ON THE SPHERE 21 The analytic expressions for the inverse operators can be easily ob- tained by making use of (4.13), (4.14), (4.19), and (4.40). Specifically, n 1 1 (a + 1) − (1 aηn+1) (a + 1) η′ ( − f˜)(η)= − f˜ , (4.47) Ua (a η )n a η − n+1 − n+1 1 a +1 (a + 1) θ + ten+1 ( − Φ)(˜ θ, t)= Φ˜ . (4.48) Va t2(1 a)+ a +1 t2 +(a + 1)2 s − ! In particular, p for a = 0: 1 1 η′ ( − f˜)(η) = f˜ , (4.49) U0 η n η | n+1| | n+1| 1 1 θ + ten+1 ( − Φ)(˜ θ, t) = Φ˜ ; (4.50) V0 √t2 +1 √t2 +1 for a = 1: n 1 1 2 − 2η′ ( − f˜)(η) = f˜ , (4.51) U1 1 η 1 η − n+1 − n+1 1 2θ + ten+1 ( − Φ)(˜ θ, t) = Φ˜ . (4.52) V1 √t2 +4 Lemma 4.8. If f C (Sn), 1 < a 1, then ∈ c a − ≤ ( f)(ξ)=( R f)(ξ) ξ Sn p . (4.53) Sa Va Ua ∀ ∈ + \ N Proof. We first show that d ( f)(ξ) π Cλ n a lim ( af)(ξ)= S , dn = . (4.54) λ 1 1 a2ξ2 Γ((n+1)/2) →− − n+1 To this end, we make usep of Lemma 2.4, according to which 1 λ Cλ 2 2 (n 2)/2 lim ( a f)(ξ)= lim ( a f)(ξ)= cn ( af)(ξ)(1 a ξn+1) − . λ 1 λ 1 →− C σn →− F − By (2.11) and (2.10), 2 2 (1 n)/2 (1 a ξn+1) − ( af)(ξ)=(Mξf)(aξn+1)= − f(η) dσ(η). F σn 1 − ξ η=Zaξn+1 · Hence c σ Cλ n n 2 2 1/2 lim ( af)(ξ)= (1 a ξn+1)− ( af)(ξ), λ 1 σn 1 − S →− − which coincides with (4.54). 22 B. RUBIN To find the limit of the right-hand side of (4.39) as λ 1, we make use of the properties of the Semyanistyi type integral→ in − Lemma 5.4. The function ρ(α, x)=( f)(x) satisfies the conditions of this Ua,λ λ=α 1 lemma, and we get − λ lim (R a,λf)(θ, t)= dn(R a, 1f)(θ, t). (4.55) λ 1 − →− U U Hence, for the right-hand side of (4.39) we have λ dn (a + 1)s lim ( a,λR a,λf)(θ,s)= (R a, 1f) θ, . (4.56) λ 1 V U √1 s2 U − √1 s2 →− − − Comparing (4.54) with (4.56), we obtain the result. To extend the factorization (4.67) to Lebesgue integrable functions and obtain the corresponding inversion formula, some auxiliary esti- mates will be needed. We recall the notation from (4.22): a(a+1)+D(x) P (x)= ,D(x)=√a+1 x 2(1 a)+a+1. (4.57) x 2 +(a+1)2 | | − | | p Lemma 4.9. Let 1 p< , 1 < a 1, λ R. Then ≤ ∞ − ≤ ∈ ( f)(x) p P λ(x) dx (4.58) | Ua | RZn p λ 1+n(p 1) p λ+n(1 p) f(η) (a ηn+1) − − − − =(a+1) | | − p 1 dη (1 aηn+1) − n SZa − provided that the integral in either side exists in the Lebesgue sense. Proof. By (4.41) and (4.21), p p (f π)(x) λ f(η) l.h.s. = ◦ P (x) dx = 1 | p λ| p 1 dη. (w π)(x) P (x) [P (π− (η)] − [w(η)] − RZn ◦ SZn a Hence, by (4.14) and (4.19), λ p n 1 p p a ηn+1 − (a + 1) (1 aηn+1) − l.h.s. = f(η) − − n+1 dη = r.h.s. | | a +1 (a ηn+1) n SZa − The following corollary contains important particular cases of (4.58) that will be used in the following. Here we make use of the formulas (4.9) and (4.10). As above, we assume that the integral in either side of the corresponding equality exists in the Lebesgue sense. RECONSTRUCTION OF FUNCTIONS ON THE SPHERE 23 Corollary 4.10. λ =0, 1 p< , 1 < a 1: • ≤ ∞ − ≤ p p n+p(1 n) f(η) dη − ( af)(x) dx=(a+1) n|+1 np| p 1 ; | U | (a ηn+1) − (1 aηn+1) − n n RZ SZa − − (4.59) λ =0, 1 p< , a = 1: • ≤ ∞ p p n+p(1 n) f(η) dη − ( 1f)(x) dx =2 | | n p(n 1) ; (4.60) | U | (1 ηn+1) − − RZn SZn − λ =0, 1 p< , a = 0: • ≤ ∞ p p f(η) dη ( 0f)(x) dx = | n|+1 np ; (4.61) | U | ηn+1 − n n RZ SZ− | | λ =1, p =1, 1 < a 1: • − ≤ ( f)(x) P (x) dx = f(η) dη; (4.62) Ua n n RZ SZa λ =1, p =1, a = 1: • ( f)(x) 1 U1 dx = f(η) dη; (4.63) x 2 +4 4 RZn | | SZn λ =1/2, p =1, a = 1: • ( f)(x) 1 f(η) U1 dx = dη. (4.64) x 2 +4 √2 √1 ηn+1 RZn | | SZn − p Corollary 4.11. If 1 Theorem 4.12. If 1 f 1 Sn,ε c f 1 Sn . (4.68) ||Sa ||L ( + ) ≤ || ||L ( a ) holds by Proposition 4.1. Further, we observe that there exist positive constants cε and aε, depending only on ε and a, such that aε I ( Φ)(ξ) dξ c dt Φ(θ, t) dθ (4.69) ε ≡ | Va | ≤ ε | | SZn,ε Z SnZ−1 + 0 n 1 for any function Φ(θ, t) on S − R whenever the right-hand side of this inequality is finite. Indeed, by× (2.13) and (4.42), 1 ε − 2 (n 2)/2 2 I = (1 s ) − ds ( Φ)(√1 s θ + se ) dθ ε − | Va − n+1 | Z0 SnZ−1 1 ε − 2 2 2 (n 2)/2 1 a s (a + 1)s = (1 s ) − − ds Φ θ, dθ. − 1 s2 √1 s2 Z r SnZ−1 0 − − Changing variable t =(a +1)s/√1 s2, so that s = t/ t2 +(a + 1)2, and setting − p (a + 1)(1 ε) aε = − , √ε √2 ε − we continue a ε (n 2)/2 (a + 1)2 − t2(1 a)+ a +1 I = − ε t2 +(a + 1)2 a +1 Z0 r (a + 1)2 dt Φ(θ, t) dθ. × [t2 +(a + 1)2]3/2 | | SnZ−1 Now the routine estimate gives (4.69). RECONSTRUCTION OF FUNCTIONS ON THE SPHERE 25 If Φ= R f, then (4.69) yields Ua aε aε (R f )(θ, t) I c dt (R f)(θ, t) dθ c˜ dt Ua| | dθ ε ≤ ε | Ua | ≤ ε (1 + t2)n/2 Z0 SnZ−1 Z0 SnZ−1 (R f )(θ, t) c˜ dt Ua| | dθ. ≤ ε (1 + t2)n/2 ZR SnZ−1 Hence, by (4.28) and (4.66), ( f )(x) I c˜ Ua| | dx c˜ c f(η) dη. ε ≤ ε (1 + x 2)1/2 ≤ ε a | | Rn Sn Z | | Za Thus there is a constant c(ε, a)=˜cε ca > 0 such that R f 1 Sn,ε c(ε, a) f 1 Sn . (4.70) ||Va Ua ||L ( + ) ≤ || ||L ( a ) To complete the proof, we observe that the estimates (4.70) and Sn 1 Sn (4.68) in conjunction with (4.67) on the dense subset Cc( a ) L ( a ) imply (4.67) for all f L1(Sn). ⊂ ∈ a In the case a = 1, that was excluded in the previous theorem, it is easier to prove an analogue of (4.67) directly, by making use of the stereographic projection. Theorem 4.13. If f(η) | | dη < , (4.71) √1 ηn+1 ∞ SZn − then, for almost all ξ Sn , ( f)(ξ) is finite and ∈ + S1 ( f)(ξ)=( R f)(ξ). (4.72) S1 V1 U1 Proof. We first show that ( R f )(ξ) < for almost all ξ Sn . V1 U1| | ∞ ∈ + Then we prove (4.72). This will give us the a.e. finiteness of ( 1f)(ξ), and the proof will be complete. S By (4.28) and (4.64), (R 1f)(θ, t) ( 1f)(x) f(η) U 2 n/2 d θdt= U 2 1/2 dx c | | dη < . (1+t ) ∗ (1+ x ) ≤ √1 ηn+1 ∞ Sn−Z1 R RZn | | SZn − × It follows that (R 1f)(θ, t) < for almost all (θ, t), and therefore U ∞ Sn ( 1R 1 f )(ξ) < for almost all ξ +. The latter makes all subse- quentV U operations| | ∞ with integrals well-justified.∈ n 1 Let ξ = θ sin ψ + en+1 cos ψ θ S − , 0 < ψ π/2. Because both and R commute with rotations∈ about the x≤ -axis, it suffices to S1 n+1 26 B. RUBIN assume θ = en = (0,..., 0, 1, 0). Let τγ be the hyperplane containing the slice γ = γ(en, ψ), and let 2 o′ = e cos ψ sin ψ + e cos ψ τ (4.73) n n+1 ∈ γ be the center of the sphere γ. We translate τγ so that o′ moves to the origin o = (0,..., 0). Then we rotate the translated plane τ o′, γ − making it coincident with the coordinate plane en⊥ and keeping the Rn 1 R R subspace − = e1 ... en fixed. Letγ ˜ en⊥ be the image of γ under this transformation.⊕ ⊕ We stretchγ ˜ up to⊂ the unit sphere Sn in Sn Rn 1 en⊥ and project stereographically onto − . Thus, can write In 1 0 sin ψ cos ψ γ = o′ + ργ,˜ ρ = − , ρ = , (4.74) 0 ρ ψ cos ψ −sin ψ ψ n 1 ( f)(e , ψ)= f(o′ +ρσ) d σ =r − f(o′ +ρrσ) dσ, r =sin ψ. S1 n γ˜ Zγ˜ SZn By (4.26) (with n replaced by n 1), we obtain − n 1 dy − ′ ( 1f)(θ, ψ)=(4r) f(o + ρrπ˜(y)) 2 n 1 , S ( y + 4) − RnZ−1 | | 4y +( y 2 4) e π˜(y)= | | − n+1 . y 2 +4 | | By (4.73) and (4.74), 2 4y sin ψ +4en sin 2ψ + en+1( y +4cos2ψ) o′ +ρrν˜(y)= | | . y 2 +4 | | Hence, n 1 ( 1f)(en, ψ) = (4 sin ψ) − S 2 4y sin ψ +4en sin 2ψ + en+1( y +4cos2ψ) dy f 2 | | 2 n 1 . × y +4 ( y +4) − RnZ−1 | | | | Changing variable y = u sin ψ and setting t = cot ψ, we represent this expression as 2 n 1 4(u+2ten)+( u+2ten 4) en+1 du − 4 f | 2 | − 2 n 1 . u +2ten +4 ( u+2ten +4) − RnZ−1 | | | | The latter is the Radon transform (R f)(e , 2t); cf. (4.45). U1 n RECONSTRUCTION OF FUNCTIONS ON THE SPHERE 27 Theorems 4.13 and 4.12 enable us to investigate the slice transform a using properties of the hyperplane Radon transform R. Some results Sof this kind are presented below. Others, like range characterization, support theorem, etc. are left to the reader. Owing to (4.67) and (4.72), inversion of a reduces to inversion of R by the formula S 1 1 1 1 a− = a− R− a− , (4.75) 1 1 S U V where a− and a− have the form (4.47) and (4.48), respectively. Nu- U V 1 merous formulas for R− are available in the literature; see [5, 7, 27] and references therein. Applicability of these formulas essentially de- pends on the class of functions. If the action of a is considered on the 1 Sn S entire space L ( a ) (for 0 a< 1) or on functions f satisfying (4.71) (for a = 1), then R must be≤ inverted on the weighted L1 space f˜(x) L˜1(Rn)= f˜ : | | dx < . (4.76) (1 + x 2)1/2 ∞ n RZn | | o Pointwise inversion formula for R on this entire space need additional technical work. Luckily, there is a plenty of inversion formulas when f˜ Lp(Rn) with 1 p Theorem 4.14. A function f can be uniquely reconstructed from af by the formulas (4.75), (4.47), and (4.48) under the following assump-S tions. In the case 1 1 1 η′ f(η)= (R− Ψ ) , (4.79) η n 0 η | n+1| | n+1| 1 θ + ten+1 Ψ0(θ, t)= ( 0f) ; √t2 +1 S √t2 +1 28 B. RUBIN for a = 1: n 1 2 − 1 2η′ f(η)= (R− Ψ ) , (4.80) 1 η 1 1 η − n+1 − n+1 2θ + ten+1 Ψ1(θ, t)=( 1f) . S √t2 +4 Here the inversion formulas (4.79) and (4.80) follow from (4.49), (4.50), (4.51), and (4.52). The following statement gives an example of the inversion formula for the Radon transform; cf. [27, Theorem 4.66]. Theorem 4.15. A function f Lp(Rn) with 1 p < n/(n 1) and n odd, can be recovered from its∈ Radon transform≤ by the formula− (n 1)/2 (1 n)/2 1 d − f(x) = lim π − (Rf)(θ, t + x θ) d θ, (4.81) t 0 ∗ → −2t dt · SnZ−1 where the limit is understood in the Lp-norm. Many inversion formulas that work for diverse classes of functions and any parity of n can be found in [27, Section 4.8]. 4.3. Slice Transforms of Zonal Functions. The case when f is a zonal function and af can be expressed as a one-dimensional Abel type integral, is worthS being mentioned separately. Looking at the Sn formula af = aR af, we observe that a takes zonal function on a to radialS functionsV onU Rn, R takes radialU functions to single variable functions on R+, and the latter are transformed by a to zonal function Sn V on +. To make things precise, we set ◦ f(η) f0(ηn+1), ( af)(x) ( a f0)( x ); (4.82) ≡ U ≡ U | | ◦ g(x) g0( x ), (Rg)(θ, t) (R g0)(t); (4.83) ≡ | | ≡ ◦ Φ(θ, t) Φ0(t), t> 0; ( aΦ)(ξ)=( a Φ0)(ξn+1). (4.84) ≡ V V Here (cf. (4.42), (4.41), (4.20)), n 1 ◦ (a + 1) P0 − (r) ( a f0)(r)= (f0 π0)(r), (4.85) U D0(r) ◦ a(a +1)+ D (r) D (r)= √a +1 r2(1 a)+ a +1, P (r)= 0 , 0 − 0 r2 +(a + 1)2 p ar2 (a + 1)D (r) π (r) Q (r)= − 0 ; 0 ≡ 0 r2 +(a + 1)2 RECONSTRUCTION OF FUNCTIONS ON THE SPHERE 29 ∞ ◦ 2 2 (n 3)/2 (R g0)(t)= σn 2 g0(r)(r t ) − rdr, t> 0; (4.86) − − Zt 2 2 ◦ 1 a s (a + 1)s ( a Φ0)(s)= − Φ0 , s> 0. (4.87) V 1 s2 √1 s2 r − − Thus we have the following statement. Sn Proposition 4.16. Let f be a zonal function on a , 1 < a 1. Suppose that f(η) f (η ), so that − ≤ ≡ 0 n+1 a f (t)(1 t2)δ dt < , (4.88) 0 − ∞ Z1 − where δ =(n 2)/2 if 1 ◦ ◦ ◦ ◦ ( a f0)(s)=( aR a f0)(s), 0 Corollary 4.17. Let f(η) f (η ). ≡ 0 n+1 (i) If a =0, then ( f)(ξ)= F ( 1 ξ2 ), where S0 0 − n+1 ρ p σn 2 2 2 (n 3)/2 − − F0(ρ)= n 2 f0( t)(ρ t ) dt. (4.90) ρ − − − Z0 (ii) If a =1, then ( f)(ξ)= F (2ξ2 1), where S1 1 n+1 − 1 (3 n)/2 1 ρ − (n 3)/2 n 3 F1(ρ)= σn 2 − f0(t)(t ρ) − (1 t) − dt. (4.91) − 2 − − Zρ 30 B. RUBIN Proof. (i) Owing to (4.43), (4.86), and (4.44), from (4.89) we have ∞ ◦ σn 2 2 2 (n 3)/2 ( 0 f0)(s)= − g0(r)(r t ) − rdr S √1 s2 − − Zt s 1 1 here t= , g0(r)= f0 √1 s2 (r2 +1)n/2 −√r2 +1 − ∞ 2 (n 3)/2 σn 2 1 2 s − rdr = − f0 r . √1 s2 −√r2 +1 − 1 s2 (r2 +1)n/2 − s/√Z1 s2 − − Changing variables, we obtain (4.90). (ii) We proceed as above, by making use of (4.45), (4.46), (4.86), and (4.89): ∞ ◦ 2 2 (n 3)/2 ( 1 f0)(s)= σn 2 g0(r)(r t ) − rdr S − − Zt n 1 2s 4 − r2 4 here t= , g0(r)= f0 − √1 s2 r2 +4 r2 +4 − ! ∞ n 1 2 2 (n 3)/2 4 − r 4 2 4s − = σn 2 f0 − r rdr. − r2 +4 r2 +4 − 1 s2 2s/√Z1 s2 − − Changing variables, we obtain (4.91). 5. Appendix It is known that if f is integrable on [0, ) and continuous from the right at t = 0, then ∞ ∞ 1 α 1 lim t − f(t) dt = f(0); (5.1) α 0+ → Γ(α) Z0 see, e.g., [27, Theorem 2.41]. We need the following generalization of this statement. Lemma 5.1. Let 0 <α0 < 1, and let u(α, t) be a measurable function on (0,α ] (0, ), satisfying the following conditions. 0 × ∞ (a) lim u(α, t)= u0. (α,t) (0,0) → RECONSTRUCTION OF FUNCTIONS ON THE SPHERE 31 (b) There exists c> 0 such that ∞ u(α, t) dt c for all α (0,α ]. (5.2) | | ≤ ∈ 0 Z0 Then ∞ 1 α 1 lim t − u(α, t) dt = u0. (5.3) α 0 → Γ(α) Z0 Proof. Choose any ε> 0 and let δ δ(ε) < 1 be such that ≡ u(α,s) u <ε whenever 0 <α + s < δ. (5.4) | − 0| We set t t 1 ϕ (t)= u(α,s) ds, ψ (t)= [u(α,s)) u )] ds, α α t − 0 Z0 Z0 so that ϕα(t)= t [ψα(t)+ u0] (5.5) and, by (5.4), ψ (t) <ε whenever 0 <α + t < δ. (5.6) | α | We have δ ∞ ∞ 1 α 1 1 α 1 1 α 1 t − u(α, t) dt = t − u(α, t)dt + t − u(α, t)dt Γ(α) Γ(α) Γ(α) Z0 Z0 Zδ = I1(α)+ I2(α). By (5.2), ∞ α 1 ∞ α 1 1 α 1 δ − c δ − I (α) t − u(α, t) dt u(α, t) dt 0 | 2 | ≤ Γ(α) | | ≤ Γ(α) | | ≤ Γ(α) → Zδ Zδ 32 B. RUBIN as α 0. To estimate I (α), we proceed as follows. → 1 δ δ δ α 1 1 α 1 t − ϕα(t) 1 α α 2 I1(α) = t − dϕα(t)= + − t − ϕα(t) dt Γ(α) Γ(α) Γ(α) Z0 0 Z0 δ α 1 α 1 δ − ϕα(δ) t − ϕα(t) 1 α α 1 = lim + − t − [ψα(t)+ u0] dt t 0 Γ(α) − → Γ(α) Γ(α) Z0 δ α 1 α δ − ϕα(δ) t [ψα(t)+ u0] 1 α α 1 = lim + − t − [ψα(t)+ u0] dt. t 0 Γ(α) − → Γ(α) Γ(α) Z0 By (5.6), the limit in this expression is zero, and therefore δ α 1 α δ − ϕα(δ) 1 α α 1 δ (1 α) I (α)= + − t − ψ (t) dt + u − . 1 Γ(α) Γ(α) α 0 Γ(α + 1) Z0 Hence δ α α 1 δ (1 α) δ − ϕα(δ) 1 α α 1 I (α) u = u − 1 + + − t − ψ (t) dt 1 − 0 0 Γ(α + 1) − Γ(α) Γ(α) α Z0 = A1(α)+ A2(α)+ A3(α). Clearly, A (α) 0, as α 0. Further, by (5.2), 1 → → δ α 1 α 1 δ − c δ − A2(α) = u(α,s) ds 0 | | Γ(α) ≤ Γ(α) → Z0 as α 0. It remains to note that by (5.6), → ε (1 α) δα A (α) − ε, | 3 | ≤ Γ(α + 1) → as α 0. Finally, → lim I1(α) f(0) ε. α 0 → | − | ≤ Since ε> 0 is arbitrarily, the proof is complete. The next statement is a consequence of Lemma 5.1. For the classical Riesz potential α n/2 α 1 g(y) dy 2 π Γ(α/2) (I g)(x)= n α , γn(α)= , (5.7) γn(α) x y − Γ((n α)/2) RZn | − | − RECONSTRUCTION OF FUNCTIONS ON THE SPHERE 33 it is known that lim (Iαg)(x)= g(x) if g is good enough; see, e.g., [27, α 0 Lemma 3.2]. A generalization→ of this fact to more general integrals 1 v(α,y) Vα(x)= n α dy (5.8) γn(α) x y − RZn | − | is the following. Lemma 5.2. Let 0 < α0 < 1, n 1, and let v be a measurable n ≥ n function on (0,α0] R . Given a point x R , suppose that there is a ball B (r)= y R×n : x y r , r > 0∈, and nonnegative functions x { ∈ | − | ≤ } c1(x) and c2(x) such that (a) v(α,y) c (x) for all (α,y) (0,α ] B (r), | | ≤ 1 ∈ 0 × x (b) v(α,y) | n |2 dy c2(x) for all α (0,α0]. (5.9) x y − ≤ ∈ n R BZx(r) | − | \ If lim v(α,y)= v(x), then lim Vα(x)= v(x). (α,y) (0,x) α 0 → → Proof. Changing variables and passing to polar coordinates, we obtain ∞ σn 1 α 1 Vα(x) = − r − dr v(α, x rθ) d θ γn(α) − ∗ Z0 SZn ∞ σn 1 Γ((n α)/2) 1 α/2 1 = − − t − u (α, t) dt , 2α+1 πn/2 Γ(α/2) x Z0 where ux(α, t)= v(α, x √t θ) d θ. (5.10) − ∗ SnZ−1 Because the constant in front of the square bracket equals 1 when α = 0, it remains to show that the expression inside the brackets falls into the scope of Lemma 5.1 with u0 = v(x). By the assumption (a), we can pass to the limit under the sign of integration in (5.10) to get 34 B. RUBIN lim ux(α, t)= v(x). The validity of (5.2) is checked as follows. (α,t) (0,0) → ∞ 2 ∞ ux(α, t) dt rdr v(α, x rθ) dθ | | ≤ σn 1 | − | Z0 − Z0 SnZ−1 2 v(α,y) = + | n |2 dy σn 1 x y − n ! − BxZ(r) R BZx(r) | − | \ 2 dy v(α,y) c1(x) n 2 + | n |2 dy . ≤ σn 1 x y − x y − n ! − BxZ(r) | − | R ZBx(r) | − | \ By the assumptions (a) and (b), the last expression is finite and the proof is complete. Below we show two applications of Lemma 5.2 that have been used in the present paper. Lemma 5.3. Let λ (Wλ)(ξ)= γn,λ w(λ, η) ξ η d η, (5.11) | · | ∗ SZn π1/2 Γ( λ/2) γ = − , 1 <λ 0, n 2, n,λ Γ((n + 1)/2)Γ((1 + λ)/2) − ≤ ≥ where w is a continuous function on [ 1, 0] Sn. If lim w(λ, η) = − × λ 1 1/→−2 w0(η), then lim (Wλ)(ξ) = cn ( w0)(ξ), where cn = π /Γ(n/2) and λ 1 F is the Funk→− transform (2.4). F Proof. By (2.12) and (2.9), 1 γn,λ σn 1 λ Wλ(ξ)= − y vξ(λ,y) dy, (5.12) σn | | Z1 − where 2 (n 2)/2 2 vξ(λ,y)=(1 y ) − w(λ, 1 y ψ + yξ) d ψ. (5.13) − − ∗ SnZξ⊥ ∩ p We can now apply Lemma 5.2 with α = λ + 1, n = 1, x = 0, to the integral (5.12), assuming that v (λ,y) extends by zero outside [ 1, 1]. ξ − RECONSTRUCTION OF FUNCTIONS ON THE SPHERE 35 By the continuity of w, lim vξ(λ,y)= w0(ψ) d ψ =( w0)(ξ) (λ,y) ( 1,0) ∗ F → − SnZξ⊥ ∩ and vξ(λ,y) const uniformly in (λ,y) [ 1, 0] [ 1, 1]. The latter| gives (a)| ≤ in Lemma 5.2. The validity∈ of− (b) in× that− lemma is trivial because v ( , ) is uniformly bounded and compactly supported ξ · · for every ξ. The constant cn can be easily calculated. Hence the result follows. Our next example is related to the Semyanistyi type integral 1 α 1 α(θ, t)= ρ(α, x) x θ t − dx, (5.14) R γ1(α) | · − | RZn 2απ1/2Γ(α/2) γ (α)= , 0 <α< 1, 1 Γ((1 α)/2) − and the relevant Radon transform (4.27). n Lemma 5.4. Let 0 <α0 < 1, Br = x R : x r , r> 0, and let ρ be a function on [0,α ] Rn satisfying{ ∈ the following| | ≤ } conditions. 0 × (a) For each α [0,α ], ρ(α, x)=0 whenever x / B ; ∈ 0 ∈ r (b) ρ is continuous on [0,α ] B ; 0 × r (c) lim ρ(α, x)= ρ0(x). α 0 → Then lim α(θ, t)=(Rρ0)(θ, t). α 0 → R Proof. By Fubini’s theorem, 1 α 1 α(θ, t)= y t − vθ(α,y) dy, (5.15) R γ1(α) | − | ZR where vθ(α,y)= ρ(α, x) dx = ρ(α,yθ + u) dθu. (5.16) x θZ=y θZ⊥ · Our aim is to show that the integral (5.15) falls into the scope of Lemma 5.2 with n = 1. If u > r, then yθ + u = u 2 + u 2 > r and therefore, by (a), ρ(α,yθ| | + u) = 0. Hence| | | | | | p vθ(α,y)= ρ(α,yθ + u) dθu. (5.17) ⊥ θ ZBr ∩ 36 B. RUBIN This equality, in conjunction with (b) and (c), allows us to pass to the limit in (5.16) and obtain lim vθ(α,y)= ρ0(tθ + u) dθu =(Rρ0)(θ, t). (5.18) (α,y) (0,t) → θZ⊥ The latter gives the limit assumption in Lemma 5.2. Further, by (5.17), vθ(α,y) ρ(α,yθ + u) dθu c vn 1(r), (5.19) | | ≤ | | ≤ − ⊥ θ ZBr ∩ where c = max ρ(α, x) and vn 1(r) is the volume of the (n 1)- (α,x) [0,α0] Br | | − − ∈ × dimensional ball of radius r. This implies (a) in Lemma 5.2, uniformly in t. Note also that by (5.19), v (α,y) y t dy r dy ρ(α,yθ + u) d u | θ || − | ≤ | | θ ⊥ R [t Zr,t+r] R [t Zr,t+r] θ ZBr \ − \ − ∩ 2 2cr vn 1(r). ≤ − Thus the integral in (5.15) satisfies all conditions of Lemma 5.2, and the result follows. References 1. Abouelaz, A. and Daher, R.: Sur la transformation de Radon de la sph`ere Sd. Bull. Soc. Math. France 121, 353-382 (1993) 2. Agranovsky, M., Kuchment, P., Kunyansky, L.: On reconstruction formulas and algorithms for the thermoacoustic tomography. In Photoacoustic imaging and spectroscopy, pp. 89–102, CRC Press (2009) 3. Funk, P. G.: Uber¨ Fl¨achen mit lauter geschlossenen geod¨atischen Linien”, Thesis, Georg-August-Universitt G¨ottingen, 1911. 4. Funk, P. G.: Uber¨ Fl¨achen mit lauter geschlossenen geod¨atschen Linen. Math. Ann. 74, 278–300 (1913) 5. Gel’fand, I.M., Gindikin, S.G., Graev, M.I.: Selected topics in in- tegral geometry. Translations of Mathematical Monographs, AMS, Providence, RI (2003) 6. Gindikin, S., Reeds, J., Shepp, L.: Spherical tomography and spherical integral geometry. In Tomography, impedance imaging, and integral geometry (South Hadley, MA, 1993), 83-92, Lectures in Appl. Math., 30, Amer. Math. Soc., Providence, RI (1994) 7. Helgason, S.: Integral geometry and Radon transform. Springer, New York-Dordrecht-Heidelberg-London (2011) RECONSTRUCTION OF FUNCTIONS ON THE SPHERE 37 8. Hielscher, R. and Quellmalz, M.: Reconstructing a function on the sphere from its means along vertical slices. Inverse Probl. Imaging 10(3), 711-739, (2016) 9. Hristova, Y., Moon, S. Steinhauer, D.: A Radon-type transform arising in photoacoustic tomography with circular detectors: spher- ical geometry. Inverse Probl. Sci. Eng. 24 , no. 6, 974-989 (2016) 10. John, F.: Plane Waves and Spherical Means Applied to Partial Differential Equations. Reprint of the 1955 original. Dover Publi- cations, Inc., Mineola, NY (2004) 11. Kazantsev, S. G.: Funk-Minkowski transform and spherical con- volution of Hilbert type in reconstructing functions on the sphere. Preprint, arXiv:1806.06672 (2018) 12. Kuchment, P.: The Radon transform and medical imaging. CBMS- NSF Regional Conference Series in Applied Mathematics 85. Soci- ety for Industrial and Applied Mathematics (SIAM), Philadelphia, PA (2014) 13. Minkowski, H.: Uber¨ die K¨orper konstanter Breite [in Russian]. Mat. Sbornik 25, 505–508 (1904); German translation in Gesam- melte Abhandlungen 2, Bd. Teubner, Leipzig, 277–279 (1911) 14. Oberlin, D. M. and Stein, E. M.: Mapping properties of the Radon transform. Indiana Univ. Math. J. 31, 641–650 (1982) 15. Olafsson,´ G., Pasquale, A., and Rubin, B.: Analytic and group- theoretic aspects of the cosine transform, Contemp. Math., 598, 167–188 (2013) 16. Ournycheva, E. and Rubin, B.: Semyanistyi’s integrals and Radon transforms on matrix spaces. The Journal of Fourier Analysis and Applications, 14, 60–88 (2008) 17. Palamodov, V.P.: Reconstructive Integral Geometry. Monographs in Mathematics, 98. Birkh¨auser Verlag, Basel (2004) 18. Palamodov, V.P.: Reconstruction from Integral Data. Monographs and Research Notes in Mathematics. CRC Press, Boca Raton, FL (2016) 19. Quellmalz, M.: A generalization of the Funk-Radon transform. In- verse Problems 33, no. 3, 035016, 26 pp. (2017) 20. Quellmalz, M.: The Funk-Radon transform for hyperplane sections through a common point, Preprint, arXiv:1810.08105 (2018) 21. Rubin, B.: Fractional calculus and wavelet transforms in integral geometry. Fractional Calculus and Applied Analysis 1, 193–219 (1998) 22. Rubin, B.: Inversion formulas for the spherical Radon transform and the generalized cosine transform. Advances in Appl. Math. 29, 38 B. RUBIN 471–497 (2002) 23. Rubin, B.: Radon, cosine, and sine transforms on real hyperbolic space. Advances in Math. 170, 206–223 (2002) 24. Rubin, B.: Intersection bodies and generalized cosine transforms. Adv. in Math. 218, 696–727 (2008) 25. Rubin, B.: Semyanistyi fractional integrals and Radon transforms. Contemp. Math. 598, 221–237 (2013) 26. Rubin, B.: Funk, Cosine, and Sine Transforms on Stiefel and Grass- mann manifolds, Journal of Geometric Analysis 23, Issue 3, 1441– 1497 (2013) 27. Rubin, B.: Introduction to Radon transforms: With elements of fractional calculus and harmonic analysis (Encyclopedia of Math- ematics and its Applications), Cambridge University Press (2015) 28. Rubin, B.: Radon Transforms and Gegenbauer-Chebyshev inte- grals, II; Examples, Analysis and Mathematical Physics, 7, Issue 4, 349–375 (2017) 29. Rubin, B.: The Vertical Slice Transform in Spherical Tomography. Preprint, arXiv:1807.07689 (2018) 30. Salman, Y.: An inversion formula for the spherical transform in S2 for a special family of circles of integration. Anal. Math. Phys. 6, no. 1, 43-58 (2016) 31. Salman, Y.: Recovering functions defined on the unit sphere by integration on a special family of sub-spheres. Anal. Math. Phys. 7, no. 2, 165-185 (2017) 32. Samko, S. G., Kilbas, A. A., Marichev, O. I.: Fractional Integrals and Derivatives. Theory and Applications. Gordon and Breach Sc. Publ., New York (1993) 33. Semyanistyi, V.I.: On some integral transformations in Euclidean space (In Russian). Dokl. Akad. Nauk SSSR 134, 536–539 (1960) 34. Stepanov, V.N.: The method of spherical harmonics for integral transforms on a sphere. Mathematical Structures and Modeling 2(42), 36–48 (2017) 35. Zangerl, G. and Scherzer, O.: Exact reconstruction in photoacous- tic tomography with circular integrating detectors II: spherical ge- ometry. Math. Methods Appl. Sci. 33(15): 1771-1782 (2010) Department of Mathematics, Louisiana State University, Baton Rouge, Louisiana 70803, USA E-mail address: [email protected]