ch9_4_1.m

来自「常用图像处理源码合集」· M 代码 · 共 16 行

%距离变换
center1 = -10; 
center2 = -center1; 
dist = sqrt(2*(2*center1)^2); 
radius = dist/2 * 1.4; 
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,[])
D1 = bwdist(~bw); %生成灰度取反的距离变换图像
figure, imshow(D1,[])