Приводится параллельный алгоритм для вычисления параболического преобразования Радона – суммирование двумерной функции вдоль парабол с разной кривизной и вертикальным смещением вершины. Основу вычислительного алгоритма составляет параллельная реализация быстрого преобразования Фурье на нерегулярных сетках, учитывающая особенности современных центральных процессоров. Дается описание оптимизации алгоритмов и сравнение с существующими аналогами. Приводится применение параболического преобразования Радона при обработке сейсмических данных.
The article presents parallel algorithm for computing parabolic Radon transformation – summation of 2D function along parabolas with varying curvature and vertical shift of its apex. The basis of the computational algorithm is a parallel implementation of the Unequally-Spaced Fast Fourier Transform optimized for architecture of modern CPUs. We give a description of optimizations applied to the algorithm and the results of testing and comparison with alternative existing libraries. We also present an application of the parabolic Radon transform to seismic data processing.
