(en) In this paper, we propose new multilevel optimization methods for minimizing continuously differentiable functions obtained by discretizing models for image registration problems. These multilevel schemes rely on a novel two-step Gauss-Newton method, in which a second step is computed within each iteration by minimizing a quadratic approximation of the objective function over a certain two-dimensional subspace. Numerical results on image registration problems show that the proposed methods can outperform the standard multilevel Gauss-Newton method.