MATLAB求图形的最大内切圆代码2页.doc

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; %纵坐标的对应 endend %图像在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=le