The l 1 -norm regularization has attracted attention for image reconstruction in computed tomography. The l 0 -norm of the gradients of an image provides a measure of the sparsity of gradients of the image. In this paper, we present a new combined l 1 -norm and l 0 -norm regularization model for image reconstruction from limited projection data in computed tomography. We also propose an algorithm in the algebraic framework to solve the optimization effectively using the nonmonotone alternating direction algorithm with hard thresholding method. Numerical experiments indicate that this new algorithm makes much improvement by involving l 0 -norm regularization.