Regularization of Linear Systems with Sparsity Constraints with Applications to Large Scale Inverse Problems