资源描述：
《MATLAB的血管三维重建源代码》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、图片下载2001数学建模A题附录1：图像二值矩阵的0-1互换的matlab程序代码（zhuanhua.m）functionb0=zhuanhua(b0)%图像二值矩阵的0-1互换fori=1:512forj=1:512ifb0(i,j)==1b0(i,j)=0;elseb0(i,j)=1;endendend附录2：求各切片的最大内切圆的半径及圆心坐标matlab程序代码（ff.m）function[r,zhongxindian]=ff%输出各切片最大内切圆半径及圆心坐标a=zeros(512,512);b=zeros(512
2、,512);fori=1:512forj=1:512a(i,j)=i-257;%横坐标的对应b(i,j)=j-257;%纵坐标的对应endend%图像在xyz面上的x轴、y轴坐标zhongxindian=zeros(100,2);r=zeros(100,1);fork=0:99t=strcat('f:/',int2str(i),'.bmp');b=imread(t);b=zhuanhua(b);%将01互换blunkuo=edge(b,'sobel');%提取轮廓bgujia=bwmorph(b,'skel',inf);%提
3、取骨架%寻找内切圆[x0,y0,v0]=find(b0lunkuo);[a0,b0,c0]=find(b0gujia);m=length(a0);n=length(x0);juli=zeros(m,n);cunfang=zeros(m,2);fori=1:mforj=1:np1=a0(i);q1=b0(i);p2=x0(j);q2=y0(j);juli(i,j)=sqrt((a(p1,q1)-a(p2,q2))^2+(b(p1,q1)-b(p2,q2))^2);%骨架上的各个点到轮廓的距离end[zx,zxxh]=min(j
4、uli(i,:));%骨架上一点到轮廓的最短距离即以骨架上各个点为圆心的内切园的半径cunfang(i,1)=zx;cunfang(i,2)=zxxh;end[zd,zdxh]=max(cunfang(:,1));%寻找半径中最大的半径和其对应的圆心坐标g=a0(zdxh);h=b0(zdxh);zhongxindian(k+1,1)=a(g,h);zhongxindian(k+1,2)=b(g,h);r(k+1)=zd;end附录3：通过计算不同次数多项式拟合的偏差平方和确定拟和次数的matlab程序代码（pczx.m）f
5、unctionj=pczx(z,t)%根据不同次数的多项式拟合与原图数据偏差平方和的大小来确定多项式拟和的次数delta=zeros(10,1);fork=1:10[p,s]=polyfit(z,t,k);delta(k)=s.normrend[i,j]=min(delta);附录4：根据轮廓画出血管的三维图像的matlab程序代码forb=0:99%提取原图的轮廓，根据轮廓画出血管的三维图像m1=imread([int2str(b),'.bmp']);m(:,:,b+1)=edge(m1,'sobel');endfork=
6、0:99fori=1:512forj=1:512if(m(i,j,k+1)==1)plot3(i,j,k+1,'r-.');holdonendendendendgridontitle('血管三维图')rotate3dholdoff附录5：绘制中轴线及在各平面的投影图matlab程序代码formatlongpx=polyfit(z,x,7);%x,z的7次多项式拟合x1=polyval(px,z);py=polyfit(z,y,5);%y,z的5次多项式拟合y1=polyval(py,z);figure(1);%画中心轴线图p
7、lot3(x1,y1,z)gridonxlabel('X轴');ylabel('Y轴');zlabel('Z轴');title('血管中轴线图');figure(2);%画中心轴线在xoz平面上的投影plot(z,x1,'-r')ylabel('Z轴');xlabel('X轴')title('血管中轴线XOZ平面投影图');gridonfigure(3);%画中心轴线在yoz平面上的投影plot(z,y1,'-b')xlabel('Z轴');ylabel('Y轴');title('血管中轴线YOZ平面投影图');gridon
8、figure(4);%画中心轴线在xoy平面上的投影plot(x1,y1,'-g')xlabel('X轴');ylabel('Y轴');title('血管中轴线XOY平面投影图');gridon附录6：求第pn张拟合图的轮廓的二值矩阵的matlab程序代码（dian.m）functionp
