This paper investigates the use of cascadic multiresolution methods for image deblurring. Iterations with a conjugate gradient-type method are carried out on each level, and terminated by a stopping rule based on the discrepancy principle. Prolongation is carried out by nonlinear edge-preserving operators, which are defined via PDEs associated with Perona