function gBank= GaborPyramid()
	Fmax  = 8;
	dF = 1.4;		% frequency bandwidth = 1.4 octave
	dOmega = 40;	%orientation bandwidth = 40 degree
	Ka = (2^dF - 1) / (2^dF + 1);
	Kb = tand(0.5*dOmega);
	dTheta = 22.5;
	szx=96;
	szy=96;
	nLayers = 5;
	nDir = 180/dTheta;
	gBank = struct();
	%Create the Base Layer of the Pyramid
	for f = 1:nLayers
		for d= 1:nDir 
			omega = dTheta*(d-1);
			F = Fmax/2^(f-1);
			rawFilter=spGabor(omega, F, Ka, Kb, szx, szy);
			gBank(f, d).rawFilter=imresize(rawFilter, 0.5^(f-1));
		end
	end
	return
