function B = radonatt(I,Theta,A);

%                  B=radonatt(I,Theta,A) 
% Calcule la transformée de Radon atténuée (sinogramme) d'une image carrée
% I de lorsque l'atténuation est donnée par l'image A (I et A sont des
% images carrées de même taille, de support inclus dans le disque de rayon
% maximal inclus dans le cadre de l'image). Les angles sont exprimés en
% degrés entiers (entre 0 et 359); les valeurs des angles sont données par
% le vecteur Theta.   

[M,M]=size(I); 
l=length(Theta);
B=[]; 
for k=1:l
    Arot=imrotate(A,-Theta(k),'bilinear','crop');
    Irot=imrotate(I,-Theta(k),'bilinear','crop');
    S=radonatt_aux(Arot)*2/M ;
    Irot=Irot.*exp(-S);
    B=[B,(sum(Irot)*2/M)']; 
end