% test of 2d projective
% to illustrate 2d chirplet by using datwarp on a planewave

[X Y] = meshdom (0:739, 0:479);
whos
[Xp Yp] = project2(X,Y,pi/6,853,853);
whos

Xp = Xp.';  % datfiles TV order, Matlab is Fortran order
Yp = Yp.';

save /hdtv/steve/projx Xp
save /hdtv/steve/projy Yp

disp('shelling out for a little reducing')
!/u/steve/reduce/reduce /hdtv/steve/projx.mat 23 /hdtv/steve/projx/data
!/u/steve/reduce/reduce /hdtv/steve/projy.mat 23 /hdtv/steve/projy/data