01年数学建模A题程序.pdfVIP

程序 1：计算最大内切圆半径的粗略范围。 运行前请先将程序中 “ 图片 ”改为图片在你电脑中的实际储 存路径。 运行时用鼠标在切片图中点出内切圆的大致位置； 运行后显示的 r 是半径估计区间的下限， r1 是上限。 clear,clc 图片\); % 读入一张图片 figure imshow(ima) % 显示这张图像 [j,i]=ginput(1); % 在坐标上取圆心 i=round(i);j=round(j); r=1;s=0; [xx yy]=yuan(j,i,r); for k=1:length(xx) s=s+ima(yy(k),xx(k)); end while s==0 r=r+1; [xx yy]=yuan(j,i,r); s=0; for k=1:length(xx) s=s+ima(yy(k),xx(k)); end end r=r-1 % rmin r1=sqrt((sum(sum(ones(size(ima))))-sum(sum(ima)))/pi) 程序 2：计算最大内切圆半径、圆心坐标，作投影图与三维图。 运行前同样要修改路径。 运行时要调用函数，在命令窗口当 tt=100 时程序运行完毕，运行时 间一般不超过 2 分钟。 运行结束后，工作空间中， aax，aay 为利用计算出的数据拟合后的 球心横、纵坐标。 ar 为从每张图计算出的最大内切圆半径。 clear,clc r=28;tt=0; for ii=0:99 tt=tt+1 图片\ int2str(tt-1) .bmp]); ima1=ima; ima1([1 512],:)=0;ima1(:,[1 512])=0; for i=3:3:510 for j=3:3:510 if sum(sum(ima(i-1:i+1,[j-2 j+2])))+sum(sum(ima([i-2 i+2],j-1:j+1)))==12 ima1(i-1:i+1,j-1:j+1)=0; end end end %eg=edge(ima,sobel); 求出切片的边缘点的位置 %[yeg,xeg]=find(eg); 边缘点的横纵坐标 [yeg,xeg]=find(ima1); for k=1:length(xeg) [xx yy]=yuan(xeg(k),yeg(k),r); if rem(k,5)==1 fl=floor(r/sqrt(2)); ima(yeg(k)-fl:yeg(k)+fl,xeg(k)-fl:xeg(k)+fl)=1; end for k1=1:length(xx) ima(yy(k1),xx(k1))=1; % 剔除区域 C 以外的点 end end [yc,xc]=find(1-ima); len2=500000;len22=0; for k=1:length(xc) len2=500000; for k1=1:length(xeg) t=(xc(k)-xeg(k1))*(xc(k)-xeg(k1))+(yc(k)-yeg(k1))*(yc(k)-ye g(k1)); if tlen2 len2=t; end end if len2len22 len22=len2; z=k; end end len2=500000;len22=0; for k=:: for k0=:: len2=500000; for k1=1: