clear; close all;
x = mandrill;
x = x(1:128,:);
% compute filter bank
opt_filters = struct();
opt_filters.J = 10;
opt_filters.min_margin = [25, 800];
filters = morlet_filter_bank_2d(size(x), opt_filters);
%%
opt_wavelet = struct();
tic;
[x_phi, x_psi, meta_phi, meta_psi] = wavelet_2d(x, filters, opt_wavelet);
toc;
