背景说明
自身项目需要通过C语言实现二值图像轮廓及连通区域(外区域)搜索,后期整理网上资料,发现大部分是通过OpenCV,matlab等接口函数实现,但也有通过堆栈方法实现,不过代码量较多,后面综合评估,决定动手重新设计该算法;经过短期的调试,现整理出实现算法和代码如下;算法思想
(1)依次遍历图像的每个像素(本算法以列优先搜索),若
BinaryDiff[row][col]
为1,则以当前像素点为起点,开始连通区域搜索;(2)定义搜索模板,本文采用逆时针顺序搜索方式,起点为下方的点;若当前像素点为Cur,则准确的搜索顺序如下图所示,从0->9:
(3)设Cur[0]
代表Cur下方的像素值;取0~7八个方向的最小的非零像素值,设为Cur[i]
(若存在多个并列最小,则按搜索顺序,取第一个值); 若存在,Cur[i]
即表示下次需要搜索的点,并遍历当前节点;若i
不存在,则只遍历当前节点,并直接跳到步骤(7);否则,若Cur
0~7方向的像素值都非零,则只遍历当前节点,不更新搜索节点(降低搜索权重,重新以改点进行搜索);(4)遍历当前节点:将
Cur
的像素值+1;更新搜索节点:将Cur=Cur[1]
往下递归搜索;(5)若出现
Cur -> Cur[j](j取值为:1,3,5,7)
的节点更新,则将Cur[j-1]
的像素值也做遍历,降低其搜索权重;(6)若更新后的搜索节点非起点,则执行步骤(3),否则执行步骤(7);
(7)当前区域搜索完成,判断是否还存在为遍历的搜索节点,若存在,则跳到步骤(2),继续执行;否则,执行步骤(8);
(8)所有连通区域搜索完成,算法结束;
代码实现
因代码中带了文件log打印,因此整体代码看起来较多,如果需要去掉,可以将DEBUG_MODE设置为0;另外本代码未保存搜索路径,只打印在Trace.txt文件中
#include "stdafx.h"
#include "stdio.h"
#define g_DriverTotal 13
#define g_SensorTotal 27
#define g_CellTotal (g_DriverTotal*g_SensorTotal)
#define TraceFileName "D:\\trace.txt" // 遍历路径
#define ResultFileName "D:\\result.txt" // 遍历结果
#define DEBUG_MODE 1
typedef unsigned short u8;
typedef signed short s8;
typedef int u16;
u8 BinaryDiff[g_SensorTotal][g_DriverTotal] =
{
0,0,0,0,0,0,0,0,0,0,0,0,0,
0,0,0,0,0,0,0,0,0,0,1,1,1,
0,0,0,0,0,0,0,0,0,1,1,1,1,
0,0,0,0,0,0,0,0,0,1,1,1,0,
0,0,