We introduce a simple algorithm to reduce artifacts that often appear in image restoration techniques such as Wiener filtering. The algorithm starts with the inverse filter solution and iteratively calculates the correction term. At each iteration we use an entropy gradient and an analytically calculated step size. The algorithm uses two Fourier transforms per iteration. We show both 1 -D and 2-D examples to illustrate the algorithm.