Although image restoration methods based on spectral filtering techniques are very efficient, they can be applied only to problems with fairly simple spatially invariant blurring operators. Iterative methods, however, are much more flexible; they can be very efficient for spatially invariant as well as spatially variant blurs, they can incorporate a variety of regularization techniques and boundary conditions, and they can more easily incorporate additional constraints, such as nonnegativity. This chapter describes a variety of iterative methods used in image restoration, with a particular emphasis on efficiency, convergence behavior, and implementation. Discussion of MATLAB software implementing the methods is also provided.