function T = curlab(f,l)
% returns topographic label (Jain's version)
%      T = curlab(f,l)
%
%   f: image
%   l: window size to fit bicubic

[f1,f2,f3,f4,f5,f6,f7,f8,f9,f10] = bicubicm(f,l);
H = Mcurv(f1,f2,f3,f4,f5,f6,f7,f8,f9,f10);
K = Gcurv(f1,f2,f3,f4,f5,f6,f7,f8,f9,f10);

[m,n] = size(H);
onemat = ones(m,n);
T = onemat+3*onemat.*(onemat+sign(H)) + (onemat+sign(K));
