MATLAB求图形的最大内切圆代码.docVIP

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