clear?all
close?all
tic
N?=?256;
n?=?0;%迭代次數
it0?=?ones（4*N）;%數據矩陣擴增為原來四倍，保證對旋轉過程的靈敏度
it1=it0;
sl=（0:1:100）;
slf=（0:1:99）;
angle=（0:180）‘;
nangle?=?length（angle）;
I?=?phantom（‘Modified?Shepp-Logan‘256）;
%I=imread（‘C:\Users\adminraden\Desktop\醫學成像系統大作業\CT圖像\800px-CT_ScoutView.jpg‘）;
%I=rgb2gray（I）;
%I=double（I）;
%I?=?imresize（I[256?256]）;
%I=I./255;
P?=?radon?（I?angle）;
R0?=?iradon?（P?angle）;
R0?=?imresize（R0[256?256]）;
P1?=?radon?（I?0:10:180）;
%P1=P;
%P1（2001:181）=100;
%P1（2001:181）=（P1（1991:181）+P1（2011:181））/2;
R1?=?iradon?（P1?0:10:180）;
%R1?=?iradon?（P1?angle）;
for?a?=?1:N;
????for?b?=?1:N;
????????it0?（（4*a-3）:4*a?（4*b-3）:4*b）?=?R1?（a?b）/16;
????????I1?（（4*a-3）:4*a?（4*b-3）:4*b）?=?I?（a?b）/16;
????end
end
while?（n?????n?=?n?+?1;
????for?ii?=?1:N
????????pj?=?zeros?（4*N）;
????????pj?（:（4*ii-3）:4*ii）?=?1;
????????pj1?=?imrotate?（pj?angle?（mod（n-1?nangle）+1）?‘crop‘）;%旋轉采樣
???????
????????P_pj?=?sum（sum（I1.*pj1））;
????????R_pj?=?sum（sum（it0.*pj1））;

????????it1?=?it0?+?pj1.*（P_pj-R_pj）/sum（sum（pj1））;
????????it0?=?it1;????
????end
%??????for?ii?=?N:2
%?????????pj?=?zeros?（4*N）;
%?????????pj?（:（4*ii-3）:4*ii）?=?1;
%?????????pj1?=?imrotate?（pj?angle?（mod（n-1?nangle）+1）?‘crop‘）;%旋轉采樣
%????????
%?????????P_pj?=?sum（sum（I1.*pj1））;
%?????????R_pj?=?sum（sum（it0.*pj1））;
%?
%?????????it1?=?it0?+?pj1.*（P_pj-R_pj）/sum（sum（pj1））;
%?????????it0?=?it1;????
%?????end

??for?a?=?1:N;
??????for?b?=?1:N;
????????????????R?（a?b）?=?sum?（sum?（it0?（?（4*a-3）:4*a?（4*b-3）:4*b）））;
??????end
??end
??dA=R-R0;
%?slf（n）=max（max（abs（dA）））;
slf（n）=norm（dA1）;
end

%?for?a?=?1:N;
%?????for?b?=?1:N;
%?????????????????R?（a?b）?=?sum?（sum?（it0?（?（4*a-3）:4*a?（4*b-3）:4*b）））;
%?????end
%?end
figure
subplot（221）;imshow（I）?title（‘原始圖像‘）;stem（slf）;
subplot（222）;imshow（R0）?title（‘完整的采集信號重構圖像‘）;
subplot（223）;imshow（R1）?title（‘不完整信號重構圖像‘）;
subplot（224）;imshow（R）?title（‘迭代重建圖像‘）;

toc

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2017-06-15?19:39??大作業MATLAB程序\
?????目錄???????????0??2017-06-15?19:39??大作業MATLAB程序\art算法\
?????文件????????2071??2017-01-10?22:11??大作業MATLAB程序\art算法\main.m
?????文件??????????45??2017-01-10?21:53??大作業MATLAB程序\art算法\Untitled111.m
?????文件?????????684??2016-12-22?00:51??大作業MATLAB程序\art算法\Untitled2.m
?????文件?????????202??2016-12-28?00:22??大作業MATLAB程序\art算法\Untitled3.m
?????目錄???????????0??2017-06-15?19:39??大作業MATLAB程序\sart算法\
?????文件?????????719??2016-12-24?23:29??大作業MATLAB程序\sart算法\CWM.m
?????文件?????????609??2017-01-10?22:29??大作業MATLAB程序\sart算法\main.m
?????文件?????????445??2016-12-24?23:23??大作業MATLAB程序\sart算法\sart.m
?????文件?????????737??2016-12-24?23:17??大作業MATLAB程序\sart算法\W.m