Now we want to study the range of the Radon transform. This is important for our inverse problem: when we want to reconstruct a function from its Radon transform, all we really have to work with is a finite number of error-prone samples of the Radon transform. We need to analyze this statistically to find the function “most likely” to have produced the observed data. We can only do this if we know what sort of functions can be Radon transforms of other functions.
In other words, we can only do this if we know the Range of .
A Warm-up: The Range of the Fourier Transform
For the Fourier Transform we know
But what is ? For , is analytic (we can write down its power series by expanding the series of the exponential in the definition of the Fourier transform). But an analytic function with compact support must be zero everywhere, so if has compact support, it (and ) must be identically zero.
So . We will talk about this more later when we discuss the Paley-Wiener theorem.
The Range of Acting on
Now let’s look at the Radon transform, . We know that the Radon transform of a function must be even, i.e. , so the range of must be somewhat restricted. Let’s start with .
Proposition The range of acting on is
In other words, the range of is everything it can be.
Proof We know that . First we will show that must be even (even though we’ve already seen this for ).
Now changing to
So the range of is a subset of the even functions in . Now we will go the other way and show that every even function in is in the range.
Let be an even function in . We need to find an such that .
Note that if exists, then
so we define with the Fourier transform
We need top be careful, this is a polar change of coordinates which is singular.
We need to check something first: , so for to be well defined we must have
To see this, note that
by applying the change of variables . Next we need to check that .
So is in the domain of , and is an isometric isomorphism from to .
The Range of on the Schwartz Class
Now let’s look at how acts on another space, . We have already seen that is linear and continuous.
So what is ? Is it the set of all even Schwartz functions?
Here is our first surprise. The range of is quite a bit more restricted.
In the 1960’s Helgason and Ludwig found the following. Let and consider the moments of in :
Then observe that
But for fixed , is a homogeneous polynomial of degree in . Our integral is just summing coefficients of homogeneous polynomials, and must be homogeneous of degree as well.
So if is in the range or , then must be a homogeneous polynomial of degree ! Functions in don’t generally have finite moments, so we did not see this restriction there.
It turns out that this is the only real restriction on the range
Theorem [Helgason, Ludwig]
We have already seen the necessity of these moment conditions. We’ll get an idea of the sufficiency proof before treating it rigorously.
Say (We will add the moment conditions later). We want to find an with , so with the Fourier Slice Theorem we write
It is easy to check that this function is in . We must check that is smooth at zero.
What is ?
It is . But is this independent of ?
But by our assumption, the 0-th moment is a homogeneous polynomial of degree zero in . This is a fancy way of saying it is a constant. So by our definition of , is well defined.
So what about ?
First we should ask “what should it be?”
is differentiable at zero if and only if we can approximate it with a linear function:
Writing this with we see
If we differentiate this with respect to , we would see that
In other words, should be a homogeneous polynomial of degree one in .
Now what is ?
Which, by our moment condition, is a homogeneous polynomial of degree 1 in , exactly what we were looking for.
We can continue like this, expanding the power series of in the Fourier transform:
Where is a homogenous polynomial in of degree and is the exponential function minus the first terms of its Taylor series. .We have used the moment conditions to get this expansion.
So using expression we derived above for we can write a Taylor series for :
This is not a proof of the theorem yet, and we need more information to make it rigorous (e.g. if we knew had compact support, we coul;d show that the power series converged and was analytic) . But it does show how the moment conditions appear naturally in this setting.
We will pursue a more technical, but straightforward proof in the next post.