High Resolution Radon Transform In the T-x Domain Using "intelligent" Prioritization of the Gauss-Seidel Estimation Sequence

This paper describes a technique for computing a fast, high resolution Radon transform. The technique entails first computing a preliminary forward Radon transform using the Gauss-Seidel method. Next, the algorithm incorporates a priori information into the current Radon transform estimate by first measuring the energy of each p trace in the preliminary transform (denoting the Radon transform as u(τ, p), a p trace is defined as the collection of transform points at constant p for all τ ), then using this measure to establish a p trace estimation sequence based on the descending order of the p trace energy. Finally, the current transform is estimated via Gauss-Seidel iteration using this new p trace estimation sequence. Synthetic tests show that the method achieves excellent τ -p localization and produces very good multiple attenuation.