MATLAB的血管三维重建源代码.docVIP

图片下载2001数学建模A题 附录1：图像二值矩阵的0-1互换的matlab程序代码（zhuanhua.m） function b0=zhuanhua(b0) %图像二值矩阵的0-1互换 for i=1:512 for j=1:512 if b0(i,j)==1 b0(i,j)=0; else b0(i,j)=1; end end end 附录2：求各切片的最大内切圆的半径及圆心坐标matlab程序代码（ff.m） function [r, zhongxindian]=ff %输出各切片最大内切圆半径及圆心坐标 a=zeros(512,512); b=zeros(512,512); for i=1:512 for j=1:512 a(i,j)=i-257; %横坐标的对应 b(i,j)=j-257; %纵坐标的对应 end end %图像在xyz面上的x轴、y轴坐标 zhongxindian=zeros(100,2); r=zeros(100,1); for k=0:99 t=strcat(f:/,int2str(i),.bmp); b=imread(t); b=zhuanhua(b);%将01互换 blunkuo=edge(b,sobel);%提取轮廓 bgujia=bwmorph(b,skel,inf);%提取骨架 %寻找内切圆 [x0,y0,v0]=find(b0lunkuo); [a0,b0,c0]=find(b0gujia); m=length(a0); n=length(x0); juli=zeros(m,n); cunfang=zeros(m,2); for i=1:m for j=1:n p1=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(juli(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） function j=pczx(z,t) %根据不同次数的多项式拟合与原图数据偏差平方和的大小来确定多项式拟和的次数 delta=zeros(10,1); for k=1:10 [p,s]=polyfit(z,t,k); delta(k)=s.normr end [i,j]=min(delta); 附录4：根据轮廓画出血管的三维图像的matlab程序代码 for b=0:99 %提取原图的轮廓，根据轮廓画出血管的三维图像 m1=imread([int2str(b),.bmp]); m(:,:,b+1)=edge(m1,sobel); end for k=0:99 for i=1:512 for j=1:512 if (m(i,j,k+1)==1) plot3(i,j,k+1,r-.); hold on end end end end grid on title(血管三维图) rotate3d hold off 附录5：绘制中轴线及在各平面的投影图matlab程序代码 format long px=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); %画中心轴线图 plot3(x1,y1,z) grid on xlabel(X轴); ylabel(Y轴); zlabel(Z轴); title(血管中轴线图); figure(2); %画中心轴线在xoz平面上的投影 plot(z,x1,-r) ylabel