The method of projections onto convex sets can be used to solve many problems in image restoration, e.g., restoration from phase, spectral extrapolation, and signal recovery in computer-aided tomography. However, image-restoration problems involving nonconvex constraints cannot be handled by the method of projection onto convex sets in a fashion that ensures convergence. The restoration-from-magnitude (RFM) problem is such a case. To handle the RFM as well as other nonconvex constraints, we describe an algorithm known as generalized projections and discuss its properties. When sets are nonconvex, it is possible for the algorithm to exhibit pathological behavior that is never manifest in convex projections. We introduce an error criterion called the summed-distance error (SDE) and show under what circumstances the SDE is a monotonically decreasing function of the number of iterations. Near-optimum performance of the algorithm is achieved by relaxation parameters. Comparisons with other RFM methods are furnished for synthetic imagery.
© 1984 Optical Society of America
Original Manuscript: July 15, 1983
Manuscript Accepted: May 5, 1984
Published: September 1, 1984
Aharon Levi and Henry Stark, "Image restoration by the method of generalized projections with application to restoration from magnitude," J. Opt. Soc. Am. A 1, 932-943 (1984)