In this paper, we study the low-rank matrix completion problem with Poisson observations, where only partial entries are available and the observations are in the presence of Poisson noise. We propose a novel model composed of the Kullback-Leibler (KL) divergence by using the maximum likelihood estimation of Poisson noise, and total variation (TV) and nuclear norm constraints. Here the nuclear norm and TV constraints are utilized to explore the approximate low-rankness and piecewise smoothness of the underlying matrix, respectively. The advantage of these two constraints in the proposed model is that the low-rankness and piecewise smoothness of the underlying matrix can be exploited simultaneously, and they can be regularized for many real-world image data. An upper error bound of the estimator of the proposed model is established with high probability, which is not larger than that of only TV or nuclear norm constraint. To the best of our knowledge, this is the first work to utilize both low-rank and TV constraints with theoretical error bounds for matrix completion under Poisson observations. Extensive numerical examples on both synthetic data and real-world images are reported to corroborate the superiority of the proposed approach.