by Joonas Ilmavirta, Olli Koskela, Jesse Railo
Abstract:
We present a new computed tomography (CT) method for inverting the Radon transform in 2D. The idea relies on the geometry of the flat torus, hence we call the new method Torus CT. We prove new inversion formulas for integrable functions, solve a minimization problem associated to Tikhonov regularization in Sobolev spaces and prove that the solution operator provides an admissible regularization strategy with a quantitative stability estimate. This regularization is a simple post-processing low-pass filter for the Fourier series of a phantom. We also study the adjoint and the normal operator of the X-ray transform on the flat torus. The X-ray transform is unitary on the flat torus. We have implemented the Torus CT method using Matlab and tested it with simulated data with promising results. The inversion method is meshless in the sense that it gives out a closed form function that can be evaluated at any point of interest.
Reference:
Torus computed tomography (Joonas Ilmavirta, Olli Koskela, Jesse Railo), SIAM Journal on Applied Mathematics, volume 80, number 4, 2020.
[show abstract]
[hide abstract]
We present a new computed tomography (CT) method for inverting the Radon transform in 2D. The idea relies on the geometry of the flat torus, hence we call the new method Torus CT. We prove new inversion formulas for integrable functions, solve a minimization problem associated to Tikhonov regularization in Sobolev spaces and prove that the solution operator provides an admissible regularization strategy with a quantitative stability estimate. This regularization is a simple post-processing low-pass filter for the Fourier series of a phantom. We also study the adjoint and the normal operator of the X-ray transform on the flat torus. The X-ray transform is unitary on the flat torus. We have implemented the Torus CT method using Matlab and tested it with simulated data with promising results. The inversion method is meshless in the sense that it gives out a closed form function that can be evaluated at any point of interest.
[arXiv]
Bibtex Entry:
@article{torus-ct,
author = {Joonas Ilmavirta and Olli Koskela and Jesse Railo},
title = {{Torus computed tomography}},
journal = {SIAM Journal on Applied Mathematics},
volume = 80,
number = 4,
month = aug,
year = {2020},
arxiv = {1906.05046},
url={http://users.jyu.fi/~jojapeil/pub/torus-ct.pdf},
gsid = {16426561243474227616},
doi = {10.1137/19M1268070},
abstract = {We present a new computed tomography (CT) method for inverting the Radon transform in 2D. The idea relies on the geometry of the flat torus, hence we call the new method Torus CT. We prove new inversion formulas for integrable functions, solve a minimization problem associated to Tikhonov regularization in Sobolev spaces and prove that the solution operator provides an admissible regularization strategy with a quantitative stability estimate. This regularization is a simple post-processing low-pass filter for the Fourier series of a phantom. We also study the adjoint and the normal operator of the X-ray transform on the flat torus. The X-ray transform is unitary on the flat torus. We have implemented the Torus CT method using Matlab and tested it with simulated data with promising results. The inversion method is meshless in the sense that it gives out a closed form function that can be evaluated at any point of interest.}
}