Code related to the paper "Recovering wavelet coefficients from binary samples using fast transforms", by V. Antun.
This software implements fast matrix-vector multiplication with an NxM section of the infinite change-of-basis matrix between a Walsh sampling basis and an orthonormal wavelet reconstruction basis in one and two dimensions. It is accompanied by several examples.
- SPGL1. A solver for large-scale sparse reconstruction
- Fastwht. A fast implementation of MATLAB's
fwht
-function and Walsh functions. - WL - General purpose wavelet library. Wavelet library for boundary wavelet support. Note that the branch
new_interface
is used. - CIlib - A software library for compressive imaging. This library is used for one of the two-dimensional sampling patterns.
Start by installing the dependencies listed above. In the examples
folder, you will find a set of examples used to produce the figures in the paper. To run these, it is important to set all the default parameters used for plotting. This is done by running the script etc/cww_set_detaults.m
. This script will produce a file called var/cww_defaults.mat
, which is read by all of the scripts in the examples folder. Finally, add the following folders to your MATLAB path
addpath('/path/to/cww/var');
addpath('/path/to/cww/utils');
It is possible to speed up the computations of the two-dimensional algorithm by using MATLAB's Parallel Computing Toolbox. To do so, write parfor
, instead of for
, in the script utils/cww_handle_2d.m
. The relevant line numbers are 33, 39, 51 and 57.