Fast solvers for Toeplitz systems with applications to image restoration