ct-fbp

ct扇束等角滤波反投影重建算法matlab源码 分成源图像创建 投影 滤波反投影三个文件

function tuoyuan=creat_ellipse（ellipsen）
p = zeros（n）;

xax =  （ （0:n-1）-（n-1）/2 ） / （（n-1）/2）; 
xg = repmat（xax n 1）;   % x coordinates the y coordinates are rot90（xg）

for k = 1:size（ellipse1）    
   asq = ellipse（k2）^2;       % a^2
   bsq = ellipse（k3）^2;       % b^2
   phi = ellipse（k6）*pi/180;  % rotation angle in radians
   x0 = ellipse（k4）;          % x offset
   y0 = ellipse（k5）;          % y offset
   A = ellipse（k1）;           % Amplitude change for this ellipse
   x=xg-x0;                    % Center the ellipse
   y=rot90（xg）-y0;  
   cosp = cos（phi）; 
   sinp = sin（phi）;
   idx=find（（（x.*cosp + y.*sinp）.^2）./asq + （（y.*cosp - x.*sinp）.^2）./bsq <= 1）; 
   p（idx） = p（idx） + A;
   tuoyuan=p;
end

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

     文件       3722  2012-07-28 13:35  shanshu2.m

     文件        745  2012-07-24 17:35  creat_ellipse.m

     文件        506  2012-07-25 23:03  gtheta.m

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

                 4973                    3