We consider a real-valued function on a plane of the form m(x,y,θ)=A(x,y)+Bc(x,y)cos(2θ)+Bs(x,y)sin(2θ)+Cc(x,y)cos(4θ)Cs(x,y)sin(4θ) that models anisotropic acoustic slowness (reciprocal velocity) perturbations. This “slowness function” depends on Cartesian coordinates and polar angle θ. The five anisotropic “component functions” A (x,y), Bc(x,y), Bs(x,y), Cc(x,y) and Cs(x,y) are assumed to be real-valued Schwartz functions. The “travel time” function d(u, θ) models the travel time perturbations on an indefinitely long straight-line observation path, where the line is parameterized by perpendicular distance u from the origin and polar angle θ; it is the Radon transform of m ( x, y, θ). We show that: 1) an A can always be found with the same d(u, θ) as an arbitrary (Bc,Bs) and/or an arbitrary (Cc,Cs) ; 2) a (Bc,Bs) can always be found with the same d(u, θ) as an arbitrary A, and furthermore, infinite families of them exist; 3) a (Cc,Cs) can always be found with the same d(u, θ) as an arbitrary A, and furthermore, infinite families of them exist; 4) a (Bc,Bs) can always be found with the same d(u, θ) as an arbitrary (Cc,Cs) , and vice versa; and furthermore, infinite families of them exist; and 5) given an arbitrary isotropic reference slowness function m0(x,y), “null coefficients” (Bc,Bs) can be constructed for which d(u, θ) is identically zero (and similarly for Cc,Cs ). We provide explicit methods of constructing each of these “equivalent functions”.
We consider a real-valued function on a plane of the form m ( x , y , θ ) = A ( x , y ) + B C ( x , y ) cos ( 2 θ ) + B S ( x , y ) sin ( 2 θ ) + C C ( x , y ) cos ( 4 θ ) + C S ( x , y ) sin ( 4 θ ) that models anisotropic acoustic slowness (reciprocal velocity) perturbations. This “slowness function” depends on Cartesian coordinates ( x , y ) and polar angle θ. The five anisotropic “component functions” A ( x , y ) , B C ( x , y ) , B S ( x , y ) , C C ( x , y ) and C S ( x , y ) are assumed to be real-valued Schwartz functions. The “travel time” function d ( u , θ ) models the travel time perturbations on an indefinitely long straight-line observation path, where the line is parameterized by perpendicular distance u from the origin and polar angle θ; it is the Radon transform of m ( x , y , θ ) . We show that: 1) an A can always be found with the same d ( u , θ ) as an arbitrary ( B C , B S ) and/or an arbitrary ( C C , C S ) ; 2) a ( B C , B S ) can always be found with the same d ( u , θ ) as an arbitrary A, and furthermore, infinite families of them exist; 3) a ( C C , C S ) can always be found with the same d ( u , θ ) as an arbitrary A, and furthermore, infinite families of them exist; 4) a ( B C , B S ) can always be found with the same d ( u , θ ) as an arbitrary ( C C , C S ) , and vice versa; and furthermore, infinite families of them exist; and 5) given an arbitrary isotropic reference slowness function m 0 ( x , y ) , “null coefficients” ( B C , B S ) can be constructed for which d ( u , θ ) is identically zero (and similarly for ( C C , C S ) ). We provide explicit methods of constructing each of these “equivalent functions”.
Keywords:
Geotomography, Radon Transform, Acoustic Anisotropy, Non-Uniqueness, Null Space
This paper addresses the non-uniqueness of the two-dimensional inverse Radon transform, when the real-valued function m ( x , y , θ ) being transformed is presumed to be anisotropic; that is, varying with polar angle θ as well as with position ( x , y ) . This is a common special case in geotomography, where m ( x , y , θ ) represents perturbations in anisotropic acoustic slowness (reciprocal velocity) on a plane in the Earth with respect to a homogeneous and isotropic background and where its Radon transform represents the corresponding travel time perturbations on indefinitely long straight-line observation paths [
As we will discuss further below, the presence of weak anisotropy results in slowness with three modes of angular variation: isotropic, 2θ and 4θ (with polar angle θ), which are described by a total of five spatially-varying “component functions” [
Assuming straight line observation paths of indefinite length, the travel time perturbation d ( u , θ ) associated with an isotropic slowness perturbation m ( x , y ) is computed via the Radon transform:
d ( u , θ ) = ∫ L ( u , θ ) m ( x , y ) d l ≡ R m (1)
The line integral is taken over the straight line L with arc length l , parameterized by its perpendicular distance u to the origin and the counter-clockwise angle θ that the perpendicular makes with the x-axis. We limit our discussions here to two-dimensional functions drawn from Schwatz space [
d ¯ ( k u , θ 0 ) = F u d ( u , θ 0 ) = m ¯ ( k u cos ( θ 0 ) , k u sin ( θ 0 ) ) with m ¯ ( k x , k y ) ≡ F x F y m (2)
Here the two-dimensional Fourier transform F x F y takes ( x , y ) into ( k x , k y ) . Consequently, a Radon transform has no null space; that is, there is no non-zero isotopic slowness function m ( x , y ) for which R m = 0 . (As we shall show below, the same is not true for anisotropic slowness functions).
Multiplying the Radon transform by a smooth real-valued function F ( θ ) of polar angle θ to yield F ( θ ) d ( r , θ ) only changes the overall scaling of a projection and, because of the Fourier Projection Theorem, corresponds to a fan filter in the ( k x , k y ) domain:
F ( θ ) m ¯ ¯ ( k x , k y ) with tan ( θ ) = k x / k y (3)
A fan filter satisfying | F ( θ ) | ≤ 1 does not increase the overall energy of the Fourier transform, so by Plancherel’s theorem [
The slowness function associated with weak anisotropy [
m ( x , y , θ ) = A ( x , y ) + B C ( x , y ) cos ( 2 θ ) + B S ( x , y ) sin ( 2 θ ) + C C ( x , y ) cos ( 4 θ ) + C S ( x , y ) sin ( 4 θ ) (4)
contains three “modes” of angular variation specified by a total of five spatially-varying “coefficient functions”. The coefficient function A specifies the isotropic mode; the pair of coefficient functions ( B C , B S ) specify the 2θ mode and the pair of coefficient functions ( C C , C S ) specify the 4θ mode. The travel time is the Radon transform of this slowness function:
d ( u , θ ) = R A + R [ B C cos ( 2 θ ) ] + R [ B S sin ( 2 θ ) ] + R [ C C cos ( 4 θ ) ] + R [ C S sin ( 4 θ ) ] (5)
The trigonometric functions can be moved outside the Radon transforms, since they do not vary along the transform’s straight line integration path:
d ( u , θ ) = R A + cos ( 2 θ ) R B C + sin ( 2 θ ) R B S + cos ( 4 θ ) R C C + sin ( 4 θ ) R C S (6)
Note that the term cos ( 2 θ ) R B C implies fan filtering of B ¯ ¯ C ( k x , k y ) (and similarly for the other terms containing trigonometric functions). Because of their 2θ and 4θ dependence, these fan filters obey the symmetry conditions for Fourier transforms of real functions.
Consider a case where only the isotropic component function A ( x , y ) is nonzero, and another case where only the 2θ anisotropic component functions, B c ( x , y ) and B s ( x , y ) , are nonzero. The two cases can be made to imply the same travel times with the choice:
A = R − 1 d B with d B ≡ cos ( 2 θ ) R B C + sin ( 2 θ ) R B S (7)
Since d B is just the sum of fan-filtered versions of functions whose Radon transforms are presumed to exist, we are assured that its inverse Radon transform exists, too. Hence, an A can always be found that is equivalent to an arbitrary ( B C , B S ) (in the sense of the two having the same travel time perturbation).
Similarly, for a case where only C c ( x , y ) and C s ( x , y ) are nonzero, the equivalent A is:
A = R − 1 d C with d C ≡ cos ( 4 θ ) R C C + sin ( 4 θ ) R C S (8)
Hence, an A can always be found that mimics a set of ( C C , C S ) . We demonstrate Equation (7) by proposing a ( B C , B S ) , computing the equivalent A, and noting that their Radom transforms are identical (
First, note the identity:
A ( x , y ) = A ( x , y ) cos 2 ( 2 θ ) + A ( x , y ) sin 2 ( 2 θ ) (9)
Now consider one case where only A ( x , y ) is nonzero:
m ( x , y , θ ) = A ( x , y ) cos 2 ( 2 θ ) + A ( x , y ) sin 2 ( 2 θ ) (10)
and another case where only ( B c , B s ) and are nonzero:
m ( x , y , θ ) = B C ( x , y ) cos ( 2 θ ) + B S ( x , y ) sin ( 2 θ ) (11)
By matching terms, two cases can be made to have the same travel time:
R B C = ( cos 2 ( 2 θ ) cos ( 2 θ ) R A ) or B C = R − 1 [ cos ( 2 θ ) R A ] (12)
and
R B S = ( sin 2 ( 2 θ ) sin ( 2 θ ) R A ) or B S = R − 1 [ sin ( 2 θ ) R A ] (13)
Note that this approach relies upon the zeros in the cos 2 ( 2 θ ) and sin 2 ( 2 θ ) in the numerators cancelling the zeros cos ( 2 θ ) and sin ( 2 θ ) in the denominators, so that the fractions are finite. If we could find a pair of functions c ( 2 θ ) and s ( 2 θ ) with the property c ( 2 θ ) + s ( 2 θ ) = 1 and with zeros in the appropriate locations, we could write A ( x , y ) = A ( x , y ) c ( 2 θ ) + A ( x , y ) s ( 2 θ ) , and:
B C = R − 1 ( c ( 2 θ ) cos ( 2 θ ) R A ) and B S = R − 1 ( s ( 2 θ ) sin ( 2 θ ) R A ) (14)
would follow. Evidentially, many such pairs of function exist, including:
c ( 2 θ ) = cos 2 n ( 2 θ ) and s ( 2 θ ) = 1 − cos 2 n ( 2 θ ) (15)
and
c ( 2 θ ) = 1 − sin 2 n ( 2 θ ) and s ( 2 θ ) = sin 2 n ( 2 θ ) (16)
with n ≥ 1 . The s ( 2 θ ) function in Equation (15) has sin 2 ( 2 θ ) as a factor, as can be demonstrated by writing the unity term as [ cos 2 ( 2 θ ) + sin 2 ( 2 θ ) ] n , expanding it with the Binomial Theorem [
s ( 2 θ ) = ( cos 2 ( 2 θ ) + sin 2 ( 2 θ ) ) n − cos 2 n ( 2 θ ) = ∑ k = 0 n ( k n ) cos 2 ( n − k ) ( 2 θ ) sin 2 k ( 2 θ ) − cos 2 n ( 2 θ ) = ∑ k = 1 n ( k n ) cos 2 ( n − k ) ( 2 θ ) sin 2 k ( 2 θ ) = sin 2 ( 2 θ ) ∑ k = 1 n ( k n ) cos 2 ( n − k ) ( 2 θ ) sin 2 ( k − 1 ) ( 2 θ ) (17)
Similarly, the c ( 2 θ ) in Equation (16) has cos 2 ( 2 θ ) as a factor. We conclude that a ( B C , B S ) can always be found that is equivalent to an arbitrary A, and furthermore, that infinite families of such pairs exist. By replacing 2θ with 4θ in the above argument, we conclude that a ( C C , C S ) can always be found that is equivalent to an arbitrary A, and furthermore, that infinite families of such pairs exist. We demonstrate Equations (14) and (16) by proposing a ( B C , B S ) , computing the equivalent A, and showing that their Radon transforms are identical (
We can find the ( C C , C S ) equivalent to a ( B C , B S ) by using the method of Section 4 to find the A equivalent to ( B C , B S ) and then using the method of Section 5 to find the ( C C , C S ) equivalent to that A. Similarly, we can find the ( B C , B S ) equivalent to an ( C C , C S ) by using the method of Section 4 to find the A equivalent to ( C C , C S ) and then using the method of Section 5 to find the ( B C , B S ) equivalent to that A. Recalling that the method of Section 5 is non-unique, we conclude that a ( B C , B S ) can always be found that mimics a ( C C , C S ) , and vice versa; and furthermore, infinite families of them exist.
As an example, we find a ( C C , C S ) equivalent to an exemplary ( B C , B S ) and a ( B C , B S ) equivalent to an exemplary ( C C , C S ) (
We have demonstrated in Section 5 that two slowness functions, say m 1 ( x , y , θ ) and m 2 ( x , y , θ ) , each with only a 2θ mode, can have the same travel time; that is R m 1 = R m 2 . Because the Radom transform is a linear operator, the “null slowness function” m n u l l = m 1 − m 2 must have zero Radon transform; that is, R m n u l l = 0 [
because the travel time perturbations of a 2θ image is the sum of two fan-filtered Radon transforms). The ( B C , B S ) of an exclusively 2θ null slowness function can be constructed from an isotropic reference slowness function m 0 ( x , y ) in the following way:
B C = R − 1 ( sin ( 2 θ ) d 0 ) and B S = − R − 1 ( cos ( 2 θ ) d 0 ) with d 0 = R m 0 (18)
The travel time d B associated with this ( B C , B S ) is identically zero:
d B = R [ cos ( 2 θ ) B C + sin ( 2 θ ) B s ] = cos ( 2 θ ) R B C + sin ( 2 θ ) R B s = cos ( 2 θ ) sin ( 2 θ ) d 0 − sin ( 2 θ ) cos ( 2 θ ) d 0 = 0 (19)
This argument can be extended to the 4θ mode simply by replacing 2θ with 4θ. Hence, given an isotropic reference slowness function m 0 ( x , y ) , a 2θ null slowness function, with coefficients ( B C , B S ) , can be constructed for which the travel time perturbation is identically zero (and similarly for 4θ null slowness function). We demonstrate Equation (18) by building ( B C , B S ) and ( C C , C S ) that correspond to null slowness functions from an exemplary isotropic reference slowness function (
A weakly-anisotropic acoustic slowness function on a plane has three modes of angular variability: isotropic, 2θ with 4θ (where θ is polar angle) that are described by a total of five spatially-varying coefficients. The Radon transform of the slowness determines travel time perturbations with respect to a homogeneous and isotropic background, observed in an idealized experiment with indefinitely long straight line observation paths. In this paper, we show that an arbitrary travel time perturbation can be matched by a slowness function with
only one of any of the three mode of angular variability. Consequently, anisotropic slowness functions cannot be uniquely determined from these travel times. Many slowness functions, equivalent in the sense of having the same travel times, exist and can be computed from one another by applying straightforward procedures. Furthermore, slowness functions with exclusively two-theta variability are non-unique; infinite families of equivalent coefficients exist. The same is true of exclusively four-theta slowness functions. Consequently, some exclusively 2θ and 4θ slowness functions have identically zero travel time and, consequentially, are inherently undetectable.
In geotomography, non-uniqueness is usually handled through the addition of prior information that regularizes the problem by singling out one “best” slowness function from among all slowness functions consistent with the travel time data [
I thank the students in the Columbia University 2018 Anisotropy Seminar that I taught for helpful discussion. This research was supported by the US National Science Foundation under awards OCE-0426369 and EAR 11-47742.
The author declares no conflicts of interest regarding the publication of this paper.