shepplogan.m

几个图象处理程序

%仿真头模型
ra=[0.92 0.874 0.31 0.41 0.25 0.046 0.046 0.046 0.023 0.046];
rb=[0.69 0.6624 0.11 0.16 0.21 0.046 0.046 0.023 0.023 0.023];
x0=[0 0 0.22 -0.22 0 0 0 -0.08 0 0.06];
y0=[0 -0.0184 0 0 0.35 0.1 -0.1 -0.605 -0.605 -0.605];
ang=[90 90 72 108 90 0 0 0 0 90]*pi/180;
ref=[2.0 -0.98 -0.02 -0.02 0.01 0.01 0.01 0.01 0.01 0.01];
ref0=[1.0 -0.8 -0.2 -0.2 0.1 0.1 0.1 0.1 0.1 0.1]
clf
line=256;row=256;
sl=zeros(line,row);
sll=zeros(line,row);
len=length(ra);
n=1:-2/255:-1;
m=1:-2/255:-1;
for l=1:256;
     for r=1:256;
         for i=1:len;
            m1=(m(r)-x0(i))*cos(ang(i))+(n(l)-y0(i))*sin(ang(i));
            n1=-(m(r)-x0(i))*sin(ang(i))+(n(l)-y0(i))*cos(ang(i));
            if m1.^2/ra(i).^2+n1.^2/rb(i).^2<=1;
              sl(l,r)=ref(i)+sl(l,r);
              sll(l,r)=ref0(i)+sll(l,r);
            end
       end
    end
 end
%sl采用的是课本的折射系数数据，产生的图像不清楚        
imagesc(sl)
colormap(gray)

%sll采用的phantom函数中的一组折射系数数据
figure
imagesc(sll)
colormap(gray)