医学成像系统实验报告.doc

《医学成像系统实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、医学成像系统实验报告姓名：__陆想想_______学号：_______指导老师：___陶玲____南京航空航天大学2011年11月实验一仿真头模型一、实验目的1、掌握MATLAB制图；2、绘制S-L头模型。二、实验内容1、用MATLAB绘制S-L头模型，该模型由十个位置、大小、方向、密度各异的椭圆组成。三、编程实现程序源代码：h=zeros(512,512);%ax0=256;y0=256;a=236;b=177;w=90/180*pi;fori=1:512forj=1:512c=b*b*((j-x0)*cos(w)+(i-y0)*si

2、n(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;endendEnd%bx1=256;y1=261;a=224;b=170;w=90/180*pi;fori=1:512forj=1:512c=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;endendend%cx2=312;y2=256

3、;a=79;b=28;w=72/180*pi;fori=1:512forj=1:512c=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;endendend%dx3=200;y3=256;a=105;b=41;w=108/180*pi;fori=1:512forj=1:512c=b*b*((j-x3)*cos(w)+(i-y3)*sin(w))^2+a*a*(-(j-x3)*sin(

4、w)-(i-y3)*cos(w))^2;d=a*a*b*b;if(c<=d)h(i,j)=0.5;endendend%ex4=256;y4=166;a=64;b=54;w=90/180*pi;fori=1:512forj=1:512c=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;endendend%重叠部分x3=200;y3=256;a1=105;b1=41;w1=108/180

5、*pi;fori=1:512forj=1:512c1=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;endendend%fx5=256;y5=230;a=12;b=12;w=

6、0/180*pi;fori=1:512forj=1:512c=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;endendend%重叠部分x4=256;y4=166;a1=64;b1=54;w1=90/180*pi;fori=1:512forj=1:512c1=b1*b1*((j-x4)*cos(w1)+(i-y4)*sin(w1))^2+a1*a1*(-(j-x4)*sin(w1)

7、-(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;endendend%gx6=256;y6=282;a=12;b=12;w=0/180*pi;fori=1:512forj=1:512c=b*b*((j-x6)*cos(w)+(i-y6)*sin(w))^2+a*a*(-(j-x6)*sin(w)-(i

8、-y6)*cos(w))^2;d=a*a*b*b;if(c<=d)h(i,j)=0.9;endendend%重叠部分x3=200;y3=256;a1=105;b1=41;w1=108/180*pi;fori