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

欢迎阅读

图片下载2001数学建模A题

附录1:图像二值矩阵的0-1互换的matIab程序代码(zhuanhua.m)

functionbO=zhuanhua(bO)%图像二值矩阵的0-1互换

fori=l:512

forj=l:512

ifbO(i,j)==l

bO(i,j)=O;

lsbO(i,j)=l;

nd

nd

nd

附录2:求各切片的最内切圆的半径及圆心坐标matlab程序代码(ff.m)

function[r,zhongxindian]=ff%输出各切片最内切圆半径及圆心坐标

a=zros(512,512);

b=zros(512,512);

fori=l:512

forj=l:512

a(i,j)=i-257;%横坐标的对应

b(i,j)=j-257;%纵坐标的对应

nd

nd%图像在xyz面上的x轴、y轴坐标

欢迎阅读

zhongxindian=zros(100,2);

r=zros(100,l);

fork=0:99

,,

t=strcat(f:/,int2str(i),.bmp)；

b=imrad(t);

b=zhuanhua(b);%将01互换

blunkuo=dg(b,sobl);%提取轮廓

bgujia=bwmorph(bs「kl,inf);%提取骨架

%寻找内切圆

[x0,y0,v0]=find(b01unkuo);

[a0,b0,c0]=find(b0gujia);

m=lngth(aO);

n=lngth(xO);

juli=zros(in,n);

cunfang=zros(m,2);

fori=l:m

forj=l:n

pl=a0(i);ql=b0(i);

p2=x0(j);q2=y0(j);

juli(i,j)=sqrt((a(pl,ql)-a(p2,q2))八2+(b(pl,ql)-b(p2,q2))八2);%骨架上的各个点到

轮廓的距离

nd

[zx,zxxh]=min(juli(i,:));%骨架上一点到轮廓的最短距离即以骨架上各个点为圆心

的内切园的半径

欢迎阅读

cunfang(i,l)=zx;

cunfang(i,2)=zxxh;

nd

[zd,zdxh]=max(cunfang(:,l));%寻找半径中最的半径和其对应的圆心坐标

g=aO(zdxh);h=bO(zdxh);

zhongxindian(k+1,l)=a(g,h);

zhongxindian(k+1,2)=b(g,h);

r(k+l)=zd;

nd

附录3:通过计算不同次数多项式拟合的偏差平方和确定拟和次数的matlab程序代

码(pczx.m)

functionj=pczx(z,t)%根据不同次数的多项式拟合与原图数据偏差平方和的小来

确定多项式拟和的次数

dlta=zrob(10,l)；

fork=l:10

[p,s]=polyfit(z,t,k);

dlta(k)=s.normr

nd

[i,j]=min(dlta);

附录4:根据轮廓画出血管的三维图像的matlab程序代码

forb=0:99%提取原图的轮廓，根据轮廓画出血管的三维图像

ml=imrad([int2str(b),.bmp]);

m(:,:,b+l)=dg(ml,sobl)；

nd

欢迎阅读

fork=0:99

fori=l:512

forj=l:512

if(m(i,j,k+l)==l)

,1

plot3(i,j,k+l,r-.);

holdon

nd

nd