boundaries(BW, conn, dir)冈萨雷斯图像处理源代码.doc

function B = boundaries(BW, conn, dir) %BOUNDARIES Trace object boundaries. % B = BOUNDARIES(BW) traces the exterior boundaries of objects in % the binary image BW. B is a P-by-1 cell array, where P is the % number of objects in the image. Each cell contains a Q-by-2 % matrix, each row of which contains the row and column coordinates % of a boundary pixel. Q is the number of boundary pixels for the % corresponding object. Object boundaries are traced in the % clockwise direction. % % B = BOUNDARIES(BW, CONN) specifies the connectivity to use when % tracing boundaries. CONN may be either 8 or 4. The default % value for CONN is 8. % % B = BOUNDARIES(BW, CONN, DIR) specifies the direction used for % tracing boundaries. DIR should be either cw (trace boundaries % clockwise) or ccw (trace boundaries counterclockwise). If DIR % is omitted BOUNDARIES traces in the clockwise direction. % Copyright 2002-2004 R. C. Gonzalez, R. E. Woods, S. L. Eddins % Digital Image Processing Using MATLAB, Prentice-Hall, 2004 % $Revision: 1.6 $ $Date: 2003/11/21 14:22:07 $ if nargin 3 dir = cw; end if nargin 2 conn = 8; end L = bwlabel(BW, conn); % The number of objects is the maximum value of L. Initialize the % cell array B so that each cell initially contains a 0-by-2 matrix. numObjects = max(L(:)); if numObjects 0 B = {zeros(0, 2)}; B = repmat(B, numObjects, 1); else B = {}; end % Pad label matrix with zeros. This lets us write the % boundary-following loop without worrying about going off the edge % of the image. Lp = padarray(L, [1 1], 0, both); % Compute the linear indexing offsets to take us from a pixel to its % neighbors. M = size(Lp, 1); if conn == 8 % Order is N NE E SE S SW W NW. offsets = [-1, M - 1, M, M + 1, 1, -M + 1, -M, -M-1]; else % Order is N E S W. offsets = [-1, M, 1, -M]; end % next_search_direction_lut is a lookup table. Given the direction % from pixel k to pixel k+1, wh