We present two Matlab-based routines, Radon_inverse and Radon_forward, that perform discrete inverse and forward Radon transforms. These frequency-based methods offer the user flexible choices among multiple regularization methods and path functions. They are designed for simplicity and flexibility, and the Matlab foundation ensures portability between platforms. Both Radon_inverse and Radon_forward could be readily adapted or modified for a broad range of geophysical applications.
For sample use see "example.m"
References:
Schultz, R., Gu, Y. J., 2012.
Flexible, inversion-based Matlab implementation of the Radon transform
Computers and Geosciences 52, 437-442.
doi: 10.1016/j.cageo.2012.08.013.
An, Y., Gu, Y. J., Sacchi, M., 2007.
Imaging mantle discontinuities using least-squares Radon transform
Journal of Geophysical Research 112, B10303.
doi: 10.1029/2007JB005009.
As a stipulation for the use of these codes, in part or full, we request that our efforts are acknowledged with a proper citation of the previous two papers.
Author: R. Schultz, 2012
This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or any later version.
This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details: http://www.gnu.org/licenses/