骨架提取

BOOL IsContourP(int x,int y, IplImage *Src_Img)
{
 BOOL p[10] ={0};
 int LineBytes =Src_Img->widthStep;
 BYTE *lpPtr= (BYTE*)(Src_Img->imageData+LineBytes*y)+x;

 p[2]=*(lpPtr-LineBytes) ? true:false;
 p[3]=*(lpPtr-LineBytes+1) ? true:false;
 p[4]=*(lpPtr+1) ? true:false;
 p[5]=*(lpPtr+LineBytes+1) ? true:false;
 p[6]=*(lpPtr+LineBytes) ? true:false;
 p[7]=*(lpPtr+LineBytes-1) ? true:false;
 p[8]=*(lpPtr-1) ? true:false;
 p[9]=*(lpPtr-LineBytes-1) ? true:false;

 int Np=0;//邻域不为零节点总数
 int Tp=0;//邻域节点由0变成1的次数
 for (int i=2; i<10; i++)
 {
  Np += p[i];
  int k= (i<9) ? (i+1) : 2;
  if ( p[k] -p[i]>0)
  {
   Tp++;
  }
 }
 int p246= p[2] && p[4] && p[6];
 int p468= p[4] && p[6] && p[8];

 int p24= p[2] && !p[3] && p[4] && !p[5] && !p[6] && !p[7] && !p[8] && !p[9];
 int p46= !p[2] && !p[3] && p[4] && !p[5] && p[6] && !p[7] && !p[8] && !p[9];
 int p68= !p[2] && !p[3] && !p[4] && !p[5] && p[6] && !p[7] && p[8] && !p[9];
 int p82= p[2] && !p[3] && !p[4] && !p[5] && !p[6] && !p[7] && p[8] && !p[9];

 int p782= p[2] && !p[3] && !p[4] && !p[5] && !p[6] && p[7] && p[8] && !p[9];
 int p924= p[2] && !p[3] && p[4] && !p[5] && !p[6] && !p[7] && !p[8] && p[9];
 int p346= !p[2] && p[3] && p[4] && !p[5] && p[6] && !p[7] && !p[8] && !p[9];
 int p568= !p[2] && !p[3] && !p[4] && p[5] && p[6] && !p[7] && p[8] && !p[9];

 int p689= !p[2] && !p[3] && !p[4] && !p[5] && p[6] && !p[7] && p[8] && p[9];
 int p823= p[2] && p[3] && !p[4] && !p[5] && !p[6] && !p[7] && p[8] && !p[9];
 int p245= p[2] && !p[3] && p[4] && p[5] && !p[6] && !p[7] && !p[8] && !p[9];
 int p467= !p[2] && !p[3] && p[4] && !p[5] && p[6] && p[7] && !p[8] && !p[9];

 int p2468= p24 || p46 || p68 || p82;
 int p3333= p782 || p924 || p346 || p568 || p689 || p823 || p245 || p467;

 //判定条件第一个由数字图像处理上得到，由于结果不够满意，又加上两个条件
 return ( !p246 && !p468 && (Np<7) && (Np>1) && (Tp==1) ) || p2468 || p3333;
}

void thin( IplImage *Src_Img)//细化轮廓，得到单像素轮廓
{
 int i, j;
 CvSize img_size = cvGetSize(Src_Img);
 int Remove_Num=0;

 do //循环调用，直至没有可以去掉的点
 {
  Remove_Num=0;
  for (j = 1; j < img_size.height-1; j++)
  {
   for(i = 1; i < img_size.width-1; i++)
   {
    BYTE gray_value = ((BYTE*)(Src_Img->imageData + Src_Img->widthStep*j))[i];
    if ( gray_value && IsContourP( i, j, Src_Img))//符合条件，去掉
    {
     ((BYTE*)(Src_Img->imageData + Src_Img->widthStep*j))[i]=(UCHAR)0;
     Remove_Num++;
    }  //if
   }  //for
  }  //for
 } while( Remove_Num);
}