bfrit_demo.m

来自「脊波工具和相关实验以及有关实验数据」· M 代码 · 共 9 行

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