function [B,Theta] = radon0(R,tau,sigma,I);

%                  [B,Theta]=radon0(R,tau,sigma,I) 
%
% Calcule la transformée de Radon (sinogramme) d'une image
% I=I(-R:tau:R-tau,-R:tau:R-tau) apres l'avoir tronqué pour qu'elle soit
% localisée dans le disque D(0,R). La sortie B(theta,p) est une image du
% type B(theta,p)= B(-R:p<R-tau,0:sigma:2*pi). Cette routine n'utilise pas
% imrotate. Theta désigne la liste des angles. L'angle est figuré en
% abscisse (de 0 a 2*pi). 

D=disque(R,tau);
I=D.*I;
x=-R:tau:R-tau;
MM=length(x); 
x=x';
a=ones(1,MM); 
XX=x*a;
YY=a'*x';
theta=0:sigma:pi-sigma; 
M=length(theta); 
B1=zeros(M,MM); 

for i=1:M
   for j=1:MM
      Q=cos(theta(i))*XX+ sin(theta(i))*YY+(R-(j-1)*tau)*ones(MM);
      [II,JJ]=find(abs(Q)<tau/2);
      B1(i,j)=tau*sum(diag(I(II,JJ)));
   end
end   
B1=B1';
B2=flipud(B1);
B=[B1,B2]; 
theta1=pi:sigma:2*pi-sigma;
Theta=[theta theta1]; 