医学成像系统实验报告

仿真头模型 仿真投影的实现 仿真反投影的实现

医学成像系统实验报告

姓名：__陆想想_______

学号：___030840110____

指导老师：___陶玲____

南京航空航天大学 2011年11月

仿真头模型 仿真投影的实现 仿真反投影的实现

实验一 仿真头模型

一、 实验目的

1、掌握MATLAB制图； 2、绘制S-L头模型。 二、实验内容

1、用MATLAB绘制S-L头模型，该模型由十个位置、大小、方向、密度各异的椭圆组成。 三、编程实现

程序源代码： h=zeros(512,512); %a

x0=256; y0=256; a=236; b=177;

w=90/180*pi; for i=1:512 for j=1:512

c=b*b*((j-x0)*cos(w)+(i-y0)*sin(w))^2+a*a*(-(j-x0)*sin(w)-(i-y0)*cos(w))^2; d=a*a*b*b; if(c<=d)

h(i,j)=1.0; end end End %b x1=256; y1=261; a=224; b=170;

w=90/180*pi; for i=1:512

for j=1:512

c=b*b*((j-x1)*cos(w)+(i-y1)*sin(w))^2+a*a*(-(j-x1)*sin(w)-(i-y1)*cos(w))^2; d=a*a*b*b; if(c<=d)

h(i,j)=0.7; end end end %c

x2=312; y2=256; a=79; b=28;

w=72/180*pi; for i=1:512 for j=1:512

c=b*b*((j-x2)*cos(w)+(i-y2)*sin(w))^2+a*a*(-(j-x2)*sin(w)-(i-y2)*cos(w))^2; d=a*a*b*b; if(c<=d)

h(i,j)=0.5; end end

仿真头模型 仿真投影的实现 仿真反投影的实现

end %d x3=200; y3=256; a=105; b=41;

w=108/180*pi; for i=1:512 for j=1:512

c=b*b*((j-x3)*cos(w)+(i-y3)*sin(w))^2+a*a*(-(j-x3)*sin(w)-(i-y3)*cos(w))^2; d=a*a*b*b; if(c<=d)

h(i,j)=0.5; end end end %e

x4=256; y4=166; a=64; b=54;

w=90/180*pi; for i=1:512 for j=1:512

c=b*b*((j-x4)*cos(w)+(i-y4)*sin(w))^2+a*a*(-(j-x4)*sin(w)-(i-y4)*cos(w))^2; d=a*a*b*b; if(c<=d)

h(i,j)=0.8; end end end

%重叠部分 x3=200; y3=256; a1=105; b1=41;

w1=108/180*pi; for i=1:512 for j=1:512

c1=b1*b1*((j-x3)*cos(w1)+(i-y3)*sin(w1))^2+a1*a1*(-(j-x3)*sin(w1)-(i-y3)*cos(w1))^2;

c=b*b*((j-x4)*cos(w)+(i-y4)*sin(w))^2+a*a*(-(j-x4)*sin(w)-(i-y4)*cos(w))^2; d=a*a*b*b;

d1=a1*a1*b1*b1; if(c<=d&c1<=d1) h(i,j)=0.6; end end end %f

x5=256; y5=230; a=12; b=12;

w=0/180*pi; for i=1:512 for j=1:512

c=b*b*((j-x5)*cos(w)+(i-y5)*sin(w))^2+a*a*(-(j-x5)*sin(w)-(i-y5)*cos(w))^2; d=a*a*b*b; if(c<=d)

h(i,j)=0.8; end end end

%重叠部分 x4=256; y4=166; a1=64; b1=54;

w1=90/180*pi; for i=1:512 for j=1:512

仿真头模型 仿真投影的实现 仿真反投影的实现

c1=b1*b1*((j-x4)*cos(w1)+(i-y4)*sin(w1))^2+a1*a1*(-(j-x4)*sin(w1)-(i-y4)*cos(w1))^2;

c=b*b*((j-x5)*cos(w)+(i-y5)*sin(w))^2+a*a*(-(j-x5)*sin(w)-(i-y5)*cos(w))^2; d=a*a*b*b;

d1=a1*a1*b1*b1; if(c<=d&c1<=d1) h(i,j)=0.9; end end end %g x6=256; y6=282; a=12; b=12;

w=0/180*pi; for i=1:512 for j=1:512

c=b*b*((j-x6)*cos(w)+(i-y6)*sin(w))^2+a*a*(-(j-x6)*sin(w)-(i-y6)*cos(w))^2; d=a*a*b*b; if(c<=d)

h(i,j)=0.9; end end end

%重叠部分 x3=200; y3=256; a1=105; b1=41;

w1=108/180*pi; for i=1:512 for j=1:512

c1=b1*b1*((j-x3)*cos(w1)+(i-y3)*sin(w1))^2+a1*a1*(-(j-x3)*sin(w1)-(i-y3)*cos(w1))^2;

c=b*b*((j-x6)*cos(w)+(i-y6)*sin(w))^2+a*a*(-(j-x6)*sin(w)-(i-y6)*cos(w))^2; d=a*a*b*b;

d1=a1*a1*b1*b1; if(c<=d&c1<=d1) h(i,j)=0.7; end end end %h x7=236; y7=411; a=12; b=6;

w=0/180*pi; for i=1:512 for j=1:512

c=b*b*((j-x7)*cos(w)+(i-y7)*sin(w))^2+a*a*(-(j-x7)*sin(w)-(i-y7)*cos(w))^2; d=a*a*b*b; if(c<=d)

h(i,j)=0.8; end end end %i

x8=256; y8=411; a=6; b=6;

w=0/180*pi; for i=1:512 for j=1:512

c=b*b*((j-x8)*cos(w)+(i-y8)*sin(w))^2+a*a*(-(j-x8)*sin(w)-(i-y8)*cos(w))^2; d=a*a*b*b; if(c<=d)

h(i,j)=0.8; end

仿真头模型 仿真投影的实现 仿真反投影的实现

end end %j

x9=271; y9=411; a=12; b=6;

w=90/180*pi; for i=1:512

for j=1:512

c=b*b*((j-x9)*cos(w)+(i-y9)*sin(w))^2+a*a*(-(j-x9)*sin(w)-(i-y9)*cos(w))^2; d=a*a*b*b; if(c<=d)

h(i,j)=0.8; end end end

四、实验结果

仿真头模型 仿真投影的实现 仿真反投影的实现

实验二 仿真投影的实现

一、 实验目的

将头模型进行投影得头模型投影图。 二、 实验内容

要求每隔5度进行进行一次投影。 三、 编程实现

程序源代码

%存放0~90°投影数据 g1=zeros(18,512);

%存放90~180°投影数据 g2=zeros(18,512);

%存放0~180°投影数据 g=zeros(36,512);

%存放0~90°时，像素点在投影轴各个点的个数

c1=zeros(18,512);

%存放90~180°时，像素点在投影轴各个点的个数

c2=zeros(18,512); %反投影矩阵 h1=zeros(512,512);

%反投影矩阵上得到反投影像素点在的个数

c=zeros(512,512);

%0~90°投影 for w1=1:18

w=(w1-1)*5/180*pi; for i=1:512 for j=1:512

if(w1<=10)

r1=(j-256.5)*cos(w)^2+(256.5-i)*sin(w)*cos(w)+256.500000000001;

r=round(r1); else

r1=(j-256.5)*cos(w)*sin(w)+(256.5-i)*sin(w)

^2+256.500000000001;

r=round(r1); end

if(r>0&r<513)

g1(w1,r)=g1(w1,r)+h(i,j);

c1(w1,r)=c1(w1,r)+1; end end end end

%90~180°投影 for w1=1:18

w=( …… 此处隐藏：3093字，全部文档内容请下载后查看。喜欢就下载吧 ……