Papers › Why POCS works, and how to make it better
Why POCS works, and how to make it better
Matteo Ravasi, Nick Luiken
The archive published only this paper's code-link row. Authors, date and abstract are from arXiv's metadata (CC0), read from the Kaggle arXiv metadata snapshot of 2026-09-12 where its title matched the archive's; the title is the archive's.
Projection Over Convex Sets (POCS) is one of the most widely used algorithms in geophysical data processing to interpolate seismic data. Whilst usually described as a modification of the Gerchberg-Saxton algorithm, a formal understanding of the underlying objective function and its implication for the associated optimization process is lacking to date in the literature. We show that the POCS algorithm can be interpreted as the application of the Half-Quadratic Splitting (HQS) method to the L₀ norm of the projected sought after data (where the projection can be any orthonormal transformation), constrained on the available traces. Similarly, the popular, apparently heuristic strategy of using a decaying threshold in POCS is revealed to be the result of the continuation strategy that HQS must employ to converge to a solution of the minimizer. Finally, we prove that for off-the-grid receivers, the POCS algorithms must operate in an inner-outer fashion (i.e., an inverse problem is solved at each outer iteration). In light of our new theoretical understanding, we suggest the use of a modern solver such as the Chambolle-Pock Primal-Dual algorithm and show that this can lead to a new POCS-like method with superior interpolation capabilities at nearly the same computational cost of the industry-standard POCS method.
Code
Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.
Code Syntology ran Syntology
Not run by Syntology. Nothing on this page verifies that the listed code works.
Results from the paper archive 2025-07-28
No leaderboard rows for this paper in the archive.
Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections