matlab SIRT CT重建算法

基本的SIRT重建算法，已调试，可以运行。 做的比较简单，有兴趣的话可以在这个的基础上加点东西。欢迎交流讨论~

%SIRT 算法，可以用于CT重建
%做得比较简单，只用了四个像素，六条投影路径，具体情况见附图
%像素的实际值是A=[97 50 ;66 90]
W=[0 0.1 0.1 0;0.1 0.1 0 0;0 0 0.1 0.1;0.1 0 0 0.1;0 0.1 0 0.1;0.1 0 0.1 0];
%W是6*4的参数矩阵，元素代表某条路径在某个像素中的长度，这里做的很简略，把所有经过的长度都记为0.1，实际上应该好好算算的
P=[12 ;15 ;17 ;20; 15 ;17];
%TOF是投影值
count=100;
ng=4;%像素数
delta=1e-9;

%以下是SIRT迭代
X=ones（ng1）;%slowness
Xn=ones（ng1）;
Lie=sum（W）;%W的每列之和，是1*ng的行向量
k=count;
while k>0
    P1=W*X;
    for j=1:ng
      A=W（:j）.*（P-P1）;
      B=sum（A）;
      C=0.0001+sum（Lie2）;
      Xn（j）=B/（C*Lie（j））;
      X（j）=X（j）+Xn（j）;
    end
    e=sqrt（sum（Xn.^2））;
    if e        k=0;
    else
        k=k-1;
    end
end
X

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        778  2014-11-05 11:10  SIRT4grids.m

     文件     972054  2014-11-05 11:07  投影示意图.bmp

----------- ---------  ---------- -----  ----

               972832                    2