Code covered by the BSD License
analysis(f, J, alpha, tau) % signal is first prefiltered, and then a J-level wavelet decomposition
analysis(inImg, J, alpha,... % First, project the image into the four (tensor) approximation space.
autocorr(M, alpha) % computes the autocorrelation filter A(z) for the fractional B-spline.
autocorr(M, alpha) % computes the autocorrelation filter A(z) for the fractional B-spline.
displayFigures (inImg, w,... % Plots the original and the reconstructed images, and the modulus of the
displayResults(f, recon, w) % display the original and the reconstructed signal, the real and the
exception( f ) % Checks for inconsistencies.
exception(f, J) % check and correct for inconsistency.
filters(M, alpha, tau, flag) % sets up the multiresolution spline filters.
filters(M, alpha, tau, tag) % Computes the analysis and synthesis filters corresponding to the
generate.m
generateGaborWavelets.m
postfilter(Y, L, alpha, t... % Inverts the prefiltering operation.
prefilter(im0, alpha, tau... % projects the image onto the four different approximation spaces used in
projectionFilters(M, alph... % projects the signal onto the four 'matched' approximtion spaces.
projectionFilters(M, alph... % projects the signal onto the two 'matched' approximtion spaces.
synthesis(lowpass, w, J, ... % reconstructs the signal from the lowpass signal and the complex wavelet
synthesis(lowpass, w, J, ... % reconstructs the image from the lowpass signals, and the complex wavelet
demo.m
demo.m
View all files
Multiresolution Gabor-like transforms
by
Kunal Chaudhury
11 May 2012
(Updated
14 May 2012 )
Matlab implementation of the multiresolution Gabor filters in 1 and 2 dimensions.
exception(f, J)
function [ f , flag] = exception(f, J)
%
% check and correct for inconsistency.
%
flag = 1;
[p, q] = size(f);
if p ~= 1
fprintf('The input signal must be in the format 1xM!\n');
flag = 0;
return;
end
if q ~= 2^round( log(q)/log(2) )
fprintf('The size of the input signal must be dyadic!\n');
flag = 0;
return
end
if q /2^J < 4
fprintf('Either signal is too short, or the decomposition depth is too large \n!');
flag = 0;
return;
end
Contact us