Smooth Hard-Thresholding for Singular Values with Stein's Unbiased Risk Estimate

Low-rank matrix denoising is a central primitive in patch-based image restoration and many other inverse problems. Classical SVD-based image denoising methods often choose a truncation rank by matching residual singular-value energy with an estimated noise energy, but this rule is not a finite-sampl