x = imread('lena.bmp');x = double(x) / 256;[r, m, ind] = bfrit(x, 16, 3, [1,3]);  % block-16 ridgelet transformimagesc(r);rc = compress(r, 0.95);  % zero-out 95% of coefficientsxr = ibfrit(rc, m, ind, 3, [1,3]);SNR(x, xr)figure; imshow(xr);