机器视觉入门之路（三四，经典canny边缘检测（计算机视觉里程碑））

这一节讲traceedge函数，从字面翻译，显然和边缘痕迹相关，也有人叫他染色嵌套函数，借用别人的指点，这种编程方法对图像中具有相同性质（灰度值等）的联通区域进行染色操作是很有用的。

是的，正是这个traceedge函数的启示，我领悟了找斑函数（findblob），因为我爷爷辈是染布的，后来自己的找斑函数就命名为染色函数，以兹纪念。

以下都是我调试成功，久经考验的函数：

网上traceedge函数的第一种表达（c++），

void TraceEdge(int y, int x, int nThrLow, BYTE *pResult, int *pMag, SIZE sz)  
{   //TraceEdge(i(h),  j(w),        dThrLow,     N,               M,        sz);sz.cx=w;sz.cy=h
    //对8邻域像素进行查询   
    int xNum[8] = {1,1,0,-1,-1,-1,0,1};  
    int yNum[8] = {0,1,1,1,0,-1,-1,-1};  
    LONG yy,xx,k;  
    for(int k=0;k<8;k++)  
    {  
        yy = y+yNum[k];  
        xx = x+xNum[k];  
        if(pResult[(yy*sz.cx+xx)*4]==128 && pMag[yy*sz.cx+xx]>=nThrLow )  
        {  
            //该点设为边界点   
            pResult[(yy*sz.cx+xx)*4] = 255;
                        pResult[(yy*sz.cx+xx)*4+1] = 255;
            pResult[(yy*sz.cx+xx)*4+2] = 255;
            pResult[(yy*sz.cx+xx)*4+3] = 255;
            //以该点为中心再进行跟踪   
            TraceEdge(yy,xx,nThrLow,pResult,pMag,sz);  
        }  
    }  
}  

网上traceedge函数的第二种表达（c++），

void TraceEdge(int y, int x, int nThrLow, BYTE *pResult, int *pMag, SIZE sz)  
{   //TraceEdge(i(h),  j(w),        dThrLow,     N,               M,        sz);sz.cx=w;sz.cy=h
    //癸8綟办琩т  
    //int xNum[8] = {1,1,0,-1,-1,-1,0,1};  
    //int yNum[8] = {0,1,1,1,0,-1,-1,-1};  
    //LONG yy,xx/*,k*/;  
    //for(int k=0;k<8;k++)  
    //{  
    //    yy = y+yNum[k];  
    //    xx = x+xNum[k];  
    //    if(pResult[(yy*sz.cx+xx)]==128 && pMag[yy*sz.cx+xx]>=nThrLow )  
    //    {  
    //        pResult[(yy*sz.cx+xx)] = 255;
    //        TraceEdge(yy,xx,nThrLow,pResult,pMag,sz);  
    //    }  
    //} 
    int i,j;//2013.8.29
   for(i=-1;i<2;i++)
   for(j=-1;j<2;j++)//这里没考虑图像4边界情况，因为边界灰度都事先设为0。不可能是128
    if(pResult[(y+i)*sz.cx+x+j]==128&&pMag[(y+i)*sz.cx+x+j]>nThrLow)
    {
     pResult[(y+i)*sz.cx+x+j]=255;
     TraceEdge(y+i,x+j,nThrLow,pResult,pMag,sz);
    }
}

这是我c#的改进算法：以上使用嵌套递归，这里使用堆栈。

   private void TraceEdge(int y, int x, int nThrLow, ref byte[,] pResult, ref int[,] pMag)
         {
             Point start = new Point(x, y);
             Stack<Point> s = new Stack<Point>();
             int[] xNum = new int[8] { 1, 1, 0, -1, -1, -1, 0, 1 };
             int[] yNum = new int[8] { 0, 1, 1, 1, 0, -1, -1, -1 };
             long xx = 0, yy = 0;
             s.Push(start);
             do
             {
                 //Point pt1 = s.ElementAt(0); 
                 Point pt1 = s.Peek();
                 s.Pop();
                 for (int k = 0; k < 8; k++)
                 {
                     yy = y + yNum[k];
                     xx = x + xNum[k];
                     if (pResult[yy, xx] == 128 && pMag[yy, xx] >= nThrLow)
                     {
                         pResult[yy, xx] = 255;
                         Point pt = new Point((int)xx, (int)yy);
                         s.Push(pt);
                     }
                 }
             } while (s.Count != 0);

         }

以下是我的发展而来的findblob（找斑）函数（c++），C#博文中已经存在（染色函数），此处略去，但可以对比着看

void CBufferBase::FindBlob(BYTE *tempImage,int w,int h)//process binary image.8 bits one pixel.
 {
     for(int i=0;i<h;i++)
         for(int j=0;j<w;j++)
         {
             if(tempImage[i*w+j]>128)
             {
                 tempImage[i*w+j]=0;
             }
             else
             {
                 tempImage[i*w+j]=128;
             }
         }
         CPoint start;
         for(int m=0;m<h;m++)
             for(int n=0;n<w;n++)
             {
                 if(tempImage[m*w+n]==128)
                 {
                     start.x=n;
                     start.y=m;
                     break;
                 }
             }
      stack<CPoint> s;
            //癸8綟办琩т  
    int xNum[8] = {1,1,0,-1,-1,-1,0,1};  
    int yNum[8] = {0,1,1,1,0,-1,-1,-1};  
    LONG yy=0,xx=0;  

    tempImage[ start.y*w+start.x ]=255;
     

      s.push(start);

      do{
          CPoint pt1=s.top();
          s.pop();

          for(int k=0;k<8;k++)  
          {  
              yy = pt1.y+yNum[k];  
              xx =pt1. x+xNum[k];  
              if(tempImage[(yy*w+xx)]==128)  
              {  
                  tempImage[(yy*w+xx)] = 255;
                  CPoint pt(xx,yy);
                  s.push(pt);
              }  
          } 
      }while(!s.empty());
 }

没有谁能一次成功，只是跟着遛弯久了，自然就通透了。