julibh1.m

来自「包括MATLAB7.0实用指南（下）中所有例子的源代码」· M 代码 · 共 13 行

center1=-10;
center2=-center1;
dist=sqrt(2*(2*center1)^2);
radius=dist/2*1.4;
%lims=[floor(center1-1.2*radius) cell(center2+1.2*radius)];
lims=[floor(center1-1.2*radius) ceil(center2+1.2*radius)];
[x,y]=meshgrid(lims(1):lims(2));
bw1=sqrt((x-center1).^2+(y-center1).^2)<=radius;
bw2=sqrt((x-center2).^2+(y-center2).^2)<=radius;
bw=bw1|bw2;
figure, imshow(bw)
d=bwdist(~bw);
figure, imshow(d,[])