MATLAB的血管三维重建源代码

图片下载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('Z轴');

xlabel('X轴')

title('血管中轴线XOZ平面投影图');

grid on

figure(3);%画中心轴线在yoz平面上的投影

plot(z,y1,'-b')

xlabel('Z轴');

ylabel('Y轴');

title('血管中轴线YOZ平面投影图');

grid on

figure(4);%画中心轴线在xoy平面上的投影

plot(x1,y1,'-g')

xlabel('X轴');

ylabel('Y轴');