首页 > 代码库 > 图像填充算法
图像填充算法
封闭连通域的图像填充是个常见的算法,最近有机会接触到大图像的例子,做一下总结。
这类问题最基本的算法是种子填充。即先给出封闭区域内的一点,从这点出发搜索邻域,只要不到边界,就把相邻点纳入连通域,赋予填充色。边界的判断比较灵活,可以使用固定颜色,也可以用一定阈值的色彩容差,类似photoshop中的魔棒。其他更复杂的计算自然也可以。
邻域的搜索是填充的重点。最简单的算法就是递归,写出来也就几行代码,像下面的样子:
FillShape(int x, int y) { if( !IsBrim(x,y) ) { SetPixel(x,y); FillShape(x+1, y); FillShape(x-1, y); FillShape(x, y+1); FillShape(x, y-1); } }
只要当前点不是边缘,就填充,然后递归相邻的点。这个算法是4连域,改成8连域也是很简单的事。
上述算法简单直观,但存在一些问题。首先它需要逐点搜索,效率不高。其次,在实际的编程中,递归运算是很消耗资源的事情。函数的递归需要压栈,即把当前函数的地址和状态量放到系统的堆栈中,进入下一个状态。实际操作系统的资源总是有限的,如果递归的层数太多,系统的堆栈会溢出,这个是用户控制不了的。这样当我们处理稍大一些的图像时,往往会出现堆栈溢出的错误,使程序无法运行。
针对这些问题,我们先做修改,把填充目标由点改为线。也即使用线扫描的方式搜索连通域。进入一个种子点后,在X方向从左向右逐点搜索连通点,填充,直到遇到中断点停止,然后从这些连通点开始,取上下点递归,伪代码见下:
FillShape(int x, int y) { if( !IsBrim(x,y) ) { // 向左填充,FillToLeft(x,y); nleft = x; while(!IsBrim(nleft,y)) { SetPixel(nleft,y); nleft--; } // 向右填充,FillToRight(x,y); nright = x; while(!IsBrim(nright,y)) { SetPixel(nright,y); nright++; } // 上下层的点递归 for(i=nleft+1;i<right;i++) { FillShape(x,y+1); FillShape(x,y-1); } } }
这样做减少了递归次数,提高了效率。但是还没有解决递归的根本缺陷,如果遇到大区域填充,仍然可能出现堆栈溢出。根本的解决之道是放弃递归。研究表明,任何递归算法都是可以修改为使用循环的非递归算法。修改的关键是两步:一、设计并实现一个自己的栈,保存原来递归出现的中间状态量,这样资源的利用率大大提高。二、在循环处理代码中,至少实现起始层和下一层递归的功能,这样才能把中间状态压入自己的栈中以备处理。
经过修改的非递归算法如下:(为了简洁起见,我把具体的代码用函数名代替):
// 构建堆栈代码 PushStack();// 压栈 PopStack();// 出栈 SetStackEmpty();// 清空栈 int IsStackEmpty();// 判断栈是否为空 FillShape(int x, int y) { FillToLeft(x,y); FillToRight(x,y); SetStackEmpty(); PushStack(); while(!IsStackEmpty()) { PopStack(); xLeft = getstack_left(); xRight = getstack_right(); // 处理上边 y=y-1; FillToLeft(xLeft,y); i=xLeft; while( i <= xRight) { FillToRight(i,y); PushStack(); } // 处理下边 y=y+2; FillToLeft(xLeft,y); i=xLeft; while( i <= xRight) { FillToRight(i,y); PushStack(); } } }
这段代码的思路是,仍然使用扫描线进行邻域填充,使用自己的堆栈来记录中间状态。完成一条扫描线的填充后,把前一个扫描线状态压入栈,弹出时,向上下搜索相邻的扫描线段,填充,压栈。直到栈中状态都被弹出处理为止。
至此,填充算法提高了效率,也避免了系统堆栈溢出。而递归算法,更适合用于原理说明和较少层次的运算,对图像的处理应当慎用并修改之。
http://blog.sina.com.cn/s/blog_6fa6b8fe0100r4gn.html
图像填充算法