sinogram.m

几个图象处理程序

ellipse=[  0 0      0.92    0.69     90   2.0
           0 -0.0184  0.874   0.6624   90  -0.98  
           0.22 0      0.31   0.11     72   -0.02
           -0.22 0     0.41   0.16     108  -0.02
           0 0.35      0.25   0.21     90     0.01
           0 0.1       0.046   0.046     0     0.01
           0 -0.1       0.046  0.046    0    0.01
           -0.08 -0.605 0.046  0.023    0    0.01
           0 -0.605     0.023  0.023    0    0.01
           0.06 -0.605  0.046  0.023    90   0.01];
 width=200;
length=200;    
[number,para]=size(ellipse) ;
slmodel=zeros(width,length);
theta=linspace(0,pi,length);
 for i=1:number
        x0=ellipse(i,1);y0=ellipse(i,2); A=ellipse(i,3);
        B=ellipse(i,4); Alpha=ellipse(i,5)*pi/180; row=ellipse(i,6);        
           for j=1:width
                alpha=theta(j)-Alpha
                 ra=A*A*cos(alpha)*cos(alpha)+B*B*sin(alpha)*sin(alpha);
                 R=linspace(-1,1,length);                    
                     for k=1:length                 
                           Ra=R(k)-x0*cos(theta(j))-y0*sin(theta(j));                 
                  if ra-Ra*Ra<0
                     slmodel(k,j)=slmodel(k,j);
                 else PQ=2*A*B*sqrt(ra-Ra*Ra)/ra;
                      slmodel(k,j)=slmodel(k,j)+row*PQ;
             end     
         end 
     end
 end
   
     imagesc(slmodel);