Pyr_L=half_octave(Left,7,0.5); Pyr_R=half_octave(Right,7,0.5);
patchsize=5; sigma=5/6;
Lap_Pyr_L=laplacian(Pyr_L,patchsize,sigma); Lap_Pyr_R=laplacian(Pyr_R,patchsize,sigma);

G=mask(patchsize,sigma,0,0,'gau',sigma);
G=G/sum(sum(G));

[Dx, Dy]=disparity(Lap_Pyr_L{4},Lap_Pyr_R{4},G,patchsize,sigma);