图像处理基本算法 链码 边界跟踪

链码在图像提取的后期即模式识别是一个很重要的特征，比如进行数字识别或者文字识别都会用到链码的特征，而链码的提取则可以借助于边界跟踪算法获取边界序列，注意是边界序列而不是边界，边界很容易获取，但是要想把边界的点按照一定的顺序输出则要费些功夫。下面采用边界跟踪算法获取边界，并存储在堆栈中，（这里的堆栈实际是C++容器类，是虚拟堆栈）。


利用点的八邻域信息，选择下一个点作为边界点，这个算法需要选择一个开始点，可以选择图像上是目标点，在最上，最左的点。然后查看它的八邻域的点，从右下方45°的位置开始寻找，如果是目标点，将沿顺时针90°作为下一次寻找的方向，如果不是，则逆时针45°继续寻找，一旦找到重复上面的过程。

具体的步骤在算法中有讲解。

/************************************************************************/
/* 查找物体的边界，输出已排序的边界序列 适应于单一区域        */
/************************************************************************/
//若能够输出边界点的序列则是比较有用的
#include
#include 
#include 
#include 
using namespace std;



int main(){ 
	IplImage * image,*image2;
	image = cvLoadImage("E:\\image\\mapleleaf.tif",0);
	cvNamedWindow("image",1);
	cvShowImage("image",image);

	image2 = cvCreateImage(cvSize(image->width, image->height),image->depth,1);
	cvZero(image2);//image2 赋值为0
	//寻找区域的左上角点
	CvPoint startPoint = cvPoint(0,0);
	bool bFindStartpoint = false;
	int i ,j;
	unsigned char * ptr,*dst;
	stack board;//奇数位存储x坐标，偶数位存储y坐标

	//当前扫描点
	CvPoint currentPoint = cvPoint(0,0);
	//邻域的8个点的方向
	int directions[8][2] = {{0,1},{1,1},{1,0},{1,-1},{0,-1},{-1,-1},{-1,0},{-1,1}}; 
	int beginDirection = 0;
	bool bFindBoardpoint = false;//寻找到邻域的边界点的判定
	for (i = 0 ; i< image->height && bFindStartpoint == false; i++)
	{
		for (j = 0 ; j< image->width && bFindStartpoint == false; j++)
		{
			ptr = (unsigned char *)(image->imageData + i*image->widthStep + j);
			if (*ptr == 255)
			{
				startPoint.x = j;
				startPoint.y = i;
				bFindStartpoint = true;
				//cout<<"x:  " << j <<"y :  " <imageData + (currentPoint.y + directions[beginDirection][1])* image->widthStep + currentPoint.x + directions[beginDirection][0]);
			if (*ptr == 255)
			{
				bFindBoardpoint = true;
				currentPoint.x +=  directions[beginDirection][0];
				currentPoint.y  += directions[beginDirection][1];
				/************************************************************************/
				/*  此处添加序列存储的代码                    */
				/************************************************************************/
				//一、将边界存储到图片中
				dst  = (unsigned char *)image2->imageData + currentPoint.y * image2->widthStep + currentPoint.x;
				*dst = 255;

				//二、将边界点的序列存储到一个堆栈中
				board.push(currentPoint.x);
				board.push(currentPoint.y);

				if (currentPoint.x == startPoint.x  && currentPoint.y == startPoint.y )
				{
					bFindStartpoint = true;
				}
				//改变下次首先开始扫描的方向
				beginDirection -= 2;
				if (beginDirection < 0)
				{
					beginDirection += 8;
				}
				
				
				
			}
			else
			{
				beginDirection ++;
				beginDirection = beginDirection%8;
			}
		}
		//cout<<"currentPoint    "<

/************************************************************************/  
/* 轮廓跟踪算法获取物体的轮廓序列 生成边界链码   */  
/************************************************************************/  
#include  
#include   
#include   
#include   
using namespace std;  
  
  
  
int main(){   
    IplImage * image,*image2,*image3;  
    image = cvLoadImage("E:\\image\\bottle2.tif",0);  
    cvNamedWindow("image",1);  
    cvShowImage("image",image);  
  
    image2 = cvCreateImage(cvSize(image->width, image->height),image->depth,1);  
    image3 = cvCreateImage(cvSize(image->width, image->height),image->depth,1);  
    cvZero(image2);//image2 赋值为0  
    cvZero(image3);  
    //寻找区域的左上角点  
    CvPoint startPoint = cvPoint(0,0);  
    bool bFindStartpoint = false;  
    int i ,j;  
    unsigned char * ptr,*dst;  
    stack board;//奇数位存储x坐标，偶数位存储y坐标  
  
    //当前扫描点  
    CvPoint currentPoint = cvPoint(0,0);  
    //邻域的8个点的方向  
    int directions[8][2] = {{0,1},{1,1},{1,0},{1,-1},{0,-1},{-1,-1},{-1,0},{-1,1}};   
    int beginDirection = 0;  
    bool bFindBoardpoint = false;//寻找到邻域的边界点的判定  
    for (i = 0 ; i< image->height && bFindStartpoint == false; i++)  
    {  
        for (j = 0 ; j< image->width && bFindStartpoint == false; j++)  
        {  
            ptr = (unsigned char *)(image->imageData + i*image->widthStep + j);  
            if (*ptr == 255)  
            {  
                startPoint.x = j;  
                startPoint.y = i;  
                bFindStartpoint = true;  
                //cout<<"x:  " << j <<"y :  " <imageData + (currentPoint.y + directions[beginDirection][1])* image->widthStep + currentPoint.x + directions[beginDirection][0]);  
            if (*ptr == 255)  
            {  
                bFindBoardpoint = true;  
                currentPoint.x +=  directions[beginDirection][0];  
                currentPoint.y  += directions[beginDirection][1];  
                /************************************************************************/  
                /*  此处添加序列存储的代码                    */  
                /************************************************************************/  
                //一、将边界存储到图片中  
                dst  = (unsigned char *)image2->imageData + currentPoint.y * image2->widthStep + currentPoint.x;  
                *dst = 255;  
  
                //二、将边界点的序列存储到一个堆栈中  
                board.push(currentPoint.x);  
                board.push(currentPoint.y);  
  
                if (currentPoint.x == startPoint.x  && currentPoint.y == startPoint.y )  
                {  
                    bFindStartpoint = true;  
                }  
                //改变下次首先开始扫描的方向  
                beginDirection -= 2;  
                if (beginDirection < 0)  
                {  
                    beginDirection += 8;  
                }  
                  
                  
                  
            }  
            else  
            {  
                beginDirection ++;  
                beginDirection = beginDirection%8;  
            }  
        }  
        //cout<<"currentPoint    "<imageData + i*image->widthStep + j;  
         *ptr = 255;  
          
    }  
    cvNamedWindow("image3",1);  
    cvSaveImage("E:\\image\\bottle2lianma.bmp",image3);  
    cvShowImage("image3",image3);  
  
  
    cvWaitKey(0);  
    return 0;  
}  

转自http://blog.csdn.net/renshengrumenglibing/article/details/7251072