k=3.141592653/180;
N=256;
fid?=?fopen（‘sample.raw‘）;
I?=?fread（fid[256256]）;
%?imshow（I[]）;
%?title（‘image?of?the?sample.raw?file‘）;
theta3=1:1:180;%?to?calculate?radon?transform?of?R3?which?includes?90?angles
[R3xp]=radon（Itheta3）;
%?figure?imagesc（theta3xpR3）;
%?title（‘radon?transform?of?90?angles?for?Shepp-logan?head?phanton‘）;
%?xlabel（‘\theta（degree）‘）;
%?ylabel（‘X\prime‘）;
%?set（gca‘Xtick‘0:20:180）;?%?to?set?the?properties?of?axis
%?colormap（hot）;
%?colorbar;
A=zeros（180NN）;
P=zeros（767180）;
[ab]=size（R3）;
for?i=1:a%行367
????for?j=1:b%列180
????????P（ij）=R3（ij）;
????end
end
for?j=1:b
????for?i=256:767
????????P（ij）=（P（255j）+P（256j））/2;
????end
end?

for?m=1:180
????m
%????C=127.5*（1-cos（m*k）-sin（m*k））;
????M=zeros（1256）;
????for?i=1:256
????????for?j=1:511;
????????????M（i）=M（i）+P（i+jm）*f2（j1N）;%求卷積
????????end
????end
????for?i=1:256
????????for?j=1:256
????????????x=128.5+（i-1）*cos（m*k）+（j-1）*sin（m*k）;
????????????if（x>1&&x???????????????xn=floor（x）;
???????????????thita=x-xn;
???????????????if（m+1>1）
?????????????????A（m+1ij）=A（mij）+（1-thita）*M（xn）+thita*M（xn+1）;
???????????????end
????????????end
????????end
????end
end
B=zeros（256256）;
for?i=1:256
????for?j=1:256
????????B（ij）=A（180ij）;
????end
end
figureimshow（B[]）
????????

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2013-04-15?17:01??卷積反投影\
?????目錄???????????0??2013-04-15?17:01??卷積反投影\CT圖像重建\
?????目錄???????????0??2013-04-15?17:01??卷積反投影\CT圖像重建\CT圖像重建\
?????文件?????????291??2010-03-01?13:18??卷積反投影\CT圖像重建\CT圖像重建\g_myfunction.m
?????文件?????????204??2010-03-01?13:17??卷積反投影\CT圖像重建\CT圖像重建\h_myfunction.m
?????文件????????5638??2010-03-02?11:02??卷積反投影\CT圖像重建\CT圖像重建\main.m
?????文件?????????228??2010-03-01?21:27??卷積反投影\CT圖像重建\CT圖像重建\my_conv.m
?????文件????????1090??2010-03-01?18:06??卷積反投影\CT圖像重建\CT圖像重建\my_ellipse.m
?????文件?????????690??2010-02-28?19:30??卷積反投影\CT圖像重建\CT圖像重建\my_projection.m
?????文件?????????181??2010-03-01?16:23??卷積反投影\CT圖像重建\CT圖像重建\trans.m
?????文件????????3518??2012-12-02?13:04??卷積反投影\CT圖像重建\FBP.m
?????文件????????1356??2011-06-09?01:56??卷積反投影\elec385.m
?????目錄???????????0??2013-04-15?17:01??卷積反投影\矩卷積反投影\
?????文件?????????688??2011-06-23?00:42??卷積反投影\矩卷積反投影\ct.m
?????文件?????????435??2011-06-22?22:25??卷積反投影\矩卷積反投影\recover.m