Image restoration using sparse approximations of spatially varying blur operators in the wavelet domain