Optimal image restoration with the fractional Fourier transform

The classical Wiener filter, which can be implemented in O(N log N) time, is suited best for space-invariant degradation models and space-invariant signal and noise characteristics. For space-varying degradations and nonstationary processes, however, the optimal linear estimate requires O(N-2) time for implementation. Optimal filtering in fractional Fourier domains permits reduction of the error compared with ordinary Fourier domain Wiener filtering for certain types of degradation and noise while requiring only O(N log N) implementation time. The amount of reduction in error depends on the signal and noise statistics as well as on the degradation model. The largest improvements are typically obtained for chirplike degradations and noise, but other types of degradation and noise may also benefit substantially from the method (e.g., nonconstant velocity motion blur and degradation by inhomogeneous atmospheric turbulence). In any event, these reductions are achieved at no additional cost. (C) 1998 Optical Society of America.