Matlab图像练习程序：bwboundaries功能实现

基本思路：
1. 预处理后遍历图像，判断L(i,j)为前景 且 未被标记 且 该像素点所在的区域未被标记。
2. 保存当前区域标签、像素点坐标。
3. 对该像素点进行八邻域搜索，顺时针进行，初始方向设定为前进方向的反方向的下一方向（即左上）。
4. 当遇到下一坐标点时进行标记并更新当前像素点坐标（步骤2、3）。
5. 若八邻域搜索结束后找不到符合目标的像素点，则再次进行遍历，寻找新的区域，即回到步骤1。
6. 每个区域的结束标记为当前点的坐标等于保存坐标数组的第一个且数组长度大于2且下一像素点的坐标等于保存坐标数组的第二个像素点。
7. 反复进行上述操作，直到整幅图像遍历完成。

源码：

function coor = edge_search(I)
%边缘跟踪
%输入：灰度图矩阵（可用Matlab自带图像 rice.png 进行测试
%输出：米粒边缘坐标（顺时针

%%
%转化为标记矩阵
I = I >130;
L = bwlabel(I);

%标记连通图像
[row,col] = size(L);
connected = zeros(row,col);
labeled = [];
coor = {};
curr_coor = [];

%8邻域坐标
offsets = [0 -1; -1 -1; -1