clear
sl=phantom;
the?=?0:179;
[Rxp]?=?radon（slthe）;
[widthlength]=size（R）;??????
width1=width;
length1=width;
rec=zeros（width1length1）;??
flou=zeros（widthlength）;
for?i=1:length
flou（:i）?=fftshift（R（:i））;?
flou（:i）=fft（flou（:i））;??
end
o=round（width1/2）;
?for?i=1:width1
????for?j=1:length1
????????x=j-o;
????????y=o-i;????????
????????????if（x==0）???????
??????????????????if（y>=0）
?????????????????????theta=90;
?????????????????else
?????????????????????theta=-90;
??????????????????end
????????????else
????????????????theta=（atan（y/x））*180/pi;
????????????end????????????
????????????if（x<=0）
????????????????theta=180+theta;
????????????elseif（y<0）
????????????????theta=2*180+theta;
????????????end????????????
???????????r=sqrt（x^2+y^2）;
????????????if（theta<180）
??