Large-scale linear systems coming from suitable discretizations of linear inverse problems are challenging to solve. Indeed, since they are inherently ill-posed, appropriate regularization should be applied; since they are large-scale, well-established direct regularization methods (such as Tikhonov regularization) cannot often be straightforwardly employed, and iterative linear solvers should be exploited. Moreover, every regularization method crucially depends on the choice of one or more regularization parameters, which should be suitably tuned. The aim of this paper is twofold: (a) survey some well-established regularizing projection methods based on Krylov subspace methods (with a particular emphasis on methods based on the Golub-Kahan bidiagonalization algorithm), and the so-called hybrid approaches (which combine Tikhonov regularization and projection onto Krylov subspaces of increasing dimension); (b) introduce a new principled and adaptive algorithmic approach for regularization similar to specific instances of hybrid methods. In particular, the new strategy provides reliable parameter choice rules by leveraging the framework of bilevel optimization, and the links between Gauss quadrature and Golub-Kahan bidiagonalization. Numerical tests modeling inverse problems in imaging illustrate the performance of existing regularizing Krylov methods, and validate the new algorithms. K E Y W O R D S hybrid methods, imaging problems, Krylov subspace methods, large-scale linear inverse problems, regularization parameter choice rules, Tikhonov regularization 1 INTRODUCTION This paper considers linear, large-scale, discrete ill-posed problems of the form Ax true + e = b true + e = b , (1) where the matrix A ∈ R m×n represents a forward mapping that is ill-conditioned with ill-determined rank (ie, the singular values of A decay and cluster at zero without an evident gap between two consecutive ones), x true ∈ R n is the desired solution, and e ∈ R m is some unknown noise that affects the data b ∈ R m. Systems like (1) typically stem from the discretization of first-kind Fredholm integral equations, and model inverse problems arising in a variety of applications, such This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.