The Eigenspace Spectral Regularization Method for solving Discrete Ill-Posed Systems