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）
   
   asq = ellipse（2）^2;       % a^2
   bsq = ellipse（3）^2;       % b^2
   phi = ellipse（6）*pi/180;  % rotation angle in radians
   x0 = ellipse（4）;          % x offset
   y0 = ellipse（5）;          % y offset
   A = ellipse（1）;           % 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;

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

     文件       2590  2012-07-25 19:37  xiugai.m

     文件        700  2012-07-25 19:12  creat_ellipse.m

     文件        500  2012-07-25 16:38  gtheta.m

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

                 3790                    3