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);