function mask=boxmask(p,xx,yy);
[Nx Ny]=size(xx);

mask=((yy*(p(2)-p(1))/Nx+p(1))>xx)+((yy*(p(4)-p(3))/Nx+p(3))<xx);
mask=(1-mask).*(yy>p(5)).*(yy<p(6));
