CT的卷积反投影

卷积反投影的matlab程序，能够对matlab自己生成的头部幻影图进行卷积反投影

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[]）
        

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           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