MATLAB二值图像连通域快速标记demo.rar

% 二值图像连通域标记（三维） % % by HPC_ZY 191018 function label = label3d(bw) [M,N,D] = size(bw); % 待分类图 label = zeros(M,N,D); % 分类结果 class = 0; % 类别 stopflag = 0; % 停止标记 cp = 1; % 遍历起始点 stack = zeros(M*N*D/2,3,'uint32'); % 缓存栈 stackidx = 0; neib = [-1,0,0;0,-1,0;1,0,0;0,1,0;0,0,-1;0,0,1]; while 1 % 寻找种子点 for k = cp:M*N*D if bw(k) z = ceil(k/M/N); y = ceil((k-(z-1)*M*N)/M); x = k-(z-1)*M*N-(y-1)*M; stackidx = stackidx+1; stack(stackidx,:) = [x,y,z]; class = class+1; label(x,y,z) = class; bw(x,y,z) = 0; cp = k+1; break end if k == M*N*D stopflag = 1; end end % 结束 if stopflag break end % 连通域 while stackidx x = stack(stackidx,1); y = stack(stackidx,2); z = stack(stackidx,3); stackidx = stackidx-1; for n = 1:6 dx = x+neib(n,1); dy = y+neib(n,2); dz = z+neib(n,3); if dx<1||dx>M||dy<1||dy>N||dz<1||dz>D break end if bw(dx,dy,dz) stackidx = stackidx+1; stack(stackidx,:) = [dx,dy,dz]; bw(dx,dy,dz) = 0; label(dx,dy,dz) = class; end end end end end