function fout = schmidFilter(t,sigma)
    sz=fix(4*sigma);
    [x,y] = meshgrid(-sz:sz,sz:-1:-sz);
    r=sqrt(x.^2+y.^2);
    fout = filterNorm(cos((pi*t/sigma)*r).*exp(-(r.^2)/(2*sigma^2)));