clear;

f = double(imread('jon.TIF'));
%g(:,1:101) = f(:,75:175);
%h(1:101,:) = g(100:200,:);
f = discretize(f,70);
figure(1)
imagesc(f)
axis('square')
colormap(gray)


%f = normalize(f);
%subplot(2,1,2)
%imagesc(f)
%axis('square')
%colormap(gray)


%f = cropImg(f,1,16);

sze = size(f);
sze = sze(1);

fDir = dirImg(f);

fAve = tileImg(fDir, 2);
%fAve = averImg(fAve, 3);
%fAve = averImg(fAve, 3);
%fAve = averImg(fAve, 3);

fAngle = angleImg(fAve);


fClass = Poincare(fAngle, 2);


