% Symmetric rectangle function for use with DFTs
% Ty, Tx are the widths of the rect in the y and x directions
% Tx and Ty must be odd values to lead to symmetry properties
% My,Nx are the rows (y) and columns (x) of Y, respectively
function Y=irect(Ty,Tx,My,Nx)
temp=floor(Tx/2); temp=temp*2; % test for odd or even
if temp==Tx
Tx2=Tx/2; % even
else 
Tx2=(Tx-1)/2; % odd
end;
temp=floor(Ty/2); temp=temp*2; % test for odd or even
if temp==Ty
Ty2=Ty/2; % even
else 
Ty2=(Ty-1)/2; % odd
end;
rectx=zeros(1,Nx);
rectx(1:(Tx2+1))=1;
rectx(((Nx-Tx2)+1):Nx)=1;
if (My>1) & (Nx>1)
	recty=zeros(1,My);
	recty(1:(Ty2+1))=1;
	recty(((My-Ty2)+1):My)=1;
   Y=recty'*rectx;
else
	Y=rectx;
end