X-CT投影重建

清华大学，研究生医学成像系统，期末大作业~·

% 计算任意位置的椭圆的投影函数
function y=gtheta（ellipseRtheta1n）
y=zeros（1n）;
alpha=pi*ellipse（5）/180;
theta=pi*theta1/180;
r_alpha=sqrt（ellipse（3）^2*（cos（theta-alpha））^2+ellipse（4）^2*（sin（theta-alpha））^2）;
D_alpha=R-ellipse（1）*cos（theta）-ellipse（2）*sin（theta）;
for i=1:n;
    if abs（D_alpha（i）） <= r_alpha
        y（i）=2*ellipse（6）*ellipse（3）*ellipse（4）*sqrt（r_alpha^2-D_alpha（i）^2）/r_alpha^2;
    else
        y（i）=0;
    end
end

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

     文件        449  2011-12-25 23:04  Final_Project\gtheta.m

     文件        249  2011-12-25 22:47  Final_Project\hRl.m

     文件        133  2008-12-25 21:55  Final_Project\hSl.m

     文件       1625  2011-12-26 19:40  Final_Project\HWW.m

     文件       1003  2011-12-26 19:38  Final_Project\MLEM.m

     文件       1026  2011-12-26 19:39  Final_Project\RART.m

     文件       2202  2011-12-25 23:09  Final_Project\sheeplogan.m

     目录          0  2011-12-26 21:08  Final_Project

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

                 6687                    8