Primal-dual first order methods for total variation image restoration in presence of Poisson noise