Randomized matrix approximation to enhance regularized projection schemes in inverse problems