A polynomial fit preconditioner for band Toeplitz matrices in image reconstruction