c++opencv中线条细化算法

要达到的效果就是将线条尽量细化成单像素，按照论文上的Hilditch算法试了一下，发现效果不好，于是自己尝试着写了一下细化的算法，基本原理就是从上下左右四个方向向内收缩。

1.先是根据图片中的原则确定了以下16种情况

2.调试过后发现，迭代次数多了之后，原来连接着的线条会断开，分析原因如下图

3.修改了一下判断条件

4.调试过后发现还是会出现断的地方，再次分析原因如下图

5.又加了判断条件，如下图

最终实现的效果如下

 

对比图

对规则曲线的效果比较好

但是圆的效果不太好，有待改进

附上代码，测试了一天，终于弄完了，啊哈哈哈！然而后面还有更艰苦的路要走。加油！！！

 //四周细化算法
 void Refine(Mat& image)
 {
     int p[];
     int top=, down=, right=, left=;
     vector<Point> del;
     int grayvalue = ;
     int height = image.rows;   //获取图像高度
     int width = image.cols;       //获取图像宽度
     Mat *im = reinterpret_cast<Mat*>((void*)&image);    //获取像素点信息
     //上下收缩
     for (int i = ; i < height-; i++)
     {
         for (int j = ; j < width-; j++)
         {
             grayvalue = Get_gray(im, j, i);  //获取指定点灰度值
             if (grayvalue != )   //判断中心点是否为前景
             {
                 p[] = (Get_gray(im, j + , i) == ) ?  : ;
                 p[] = (Get_gray(im, j + , i - ) == ) ?  : ;
                 p[] = (Get_gray(im, j, i - ) == ) ?  : ;
                 p[] = (Get_gray(im, j - , i - ) == ) ?  : ;
                 p[] = (Get_gray(im, j - , i) == ) ?  : ;
                 p[] = (Get_gray(im, j - , i + ) == ) ?  : ;
                 p[] = (Get_gray(im, j, i + ) == ) ?  : ;
                 p[] = (Get_gray(im, j + , i + ) == ) ?  : ;
                 if (i < height - )
                     down = (Get_gray(im, j, i + ) == ) ?  : ;
                 else
                     down = ;
                 //  横向直线
                 if (p[] && (p[] || p[] || p[] || p[]) && !(p[] || p[]) && p[] ==  && down)
                 {
                     del.push_back(Point(j, i));
                 }
                 if (p[] && (p[] || p[] || p[] || p[]) && !( p[] || p[]) && p[] == )
                 {
                     del.push_back(Point(j, i));
                 }
             }
         }
     }

     for (int i = ; i < height - ; i++)
     {
         grayvalue = Get_gray(im, , i);
         if (grayvalue != )
         {
             if ( Get_gray(im, , i - ) && Get_gray(im, , i - ) && Get_gray(im, , i + )== && Get_gray(im, , i)==) //上2，上1，右上1，下1=0，右1=0
             {
                 del.push_back(Point(, i));
             }
             if (Get_gray(im, , i - ) ==  && Get_gray(im, , i + ) && Get_gray(im, , i) ==  && Get_gray(im, , i+))//上1=0，下1，右下1，右1=0，下2
             {
                 del.push_back(Point(, i));
             }
         }
         if (grayvalue != )
         {
             if (Get_gray(im, width - , i - ) && Get_gray(im, width - , i - ) && Get_gray(im, width - , i + ) ==  && Get_gray(im, width - , i) == ) //上2，上1，左上1，下1=0，左1=0
             {
                 del.push_back(Point(width - , i));
             }
             if (Get_gray(im, width - , i - ) ==  && Get_gray(im, width - , i + ) && Get_gray(im, width - , i) ==  && Get_gray(im, width - , i + ))//上1=0，下1，左下1，左1=0，下2
             {
                 del.push_back(Point(width - , i));
             }
         }
     }
     for (int i = ; i < del.size();i++)
     {
         uchar* data = image.ptr<uchar>(del[i].y);
         data[del[i].x]=;
     }

     //左右收缩
     for (int i = ; i < height - ; i++)
     {
         for (int j = ; j < width - ; j++)
         {
             grayvalue = Get_gray(im, j, i);  //获取指定点灰度值
             if (grayvalue != )   //判断中心点是否为前景
             {
                 p[] = (Get_gray(im, j + , i) == ) ?  : ;
                 p[] = (Get_gray(im, j + , i - ) == ) ?  : ;
                 p[] = (Get_gray(im, j, i - ) == ) ?  : ;
                 p[] = (Get_gray(im, j - , i - ) == ) ?  : ;
                 p[] = (Get_gray(im, j - , i) == ) ?  : ;
                 p[] = (Get_gray(im, j - , i + ) == ) ?  : ;
                 p[] = (Get_gray(im, j, i + ) == ) ?  : ;
                 p[] = (Get_gray(im, j + , i + ) == ) ?  : ;
                 if (j < width - )
                     right = (Get_gray(im, j + , i) == ) ?  : ;
                 else
                     right = ;

                 //竖直线
                 if (p[] && (p[] || p[] || p[] || p[]) && !(p[] || p[]) && p[] ==  && right)
                 {
                     del.push_back(Point(j, i));
                 }
                 if (p[] && (p[] || p[] || p[] || p[]) && !(p[] || p[]) && p[] == )
                 {
                     del.push_back(Point(j, i));
                 }

             }
         }
     }

     for (int j = ; j < width - ; j++)
     {
         grayvalue = Get_gray(im, j, );
         if (grayvalue != )
         {
             if (Get_gray(im, j - , ) ==  && Get_gray(im, j + , ) && Get_gray(im, j + , ) && Get_gray(im, j, ) ==  && Get_gray(im, j+, )) //左1=0，右1，右2，下1=0，右下1
             {
                 del.push_back(Point(j, ));
             }
             if (Get_gray(im, j - , ) && Get_gray(im, j+, )== && Get_gray(im, j, ) ==  && Get_gray(im, j-, ))//左1，右1=0，下1=0，左下1
             {
                 del.push_back(Point(j, ));
             }
         }
     }
     for (int j = ; j < width - ; j++)
     {
         grayvalue = Get_gray(im, j, height-);
         if (grayvalue != )
         {
             if (Get_gray(im, j - , height - ) ==  && Get_gray(im, j + , height - ) && Get_gray(im, j + , height - ) && Get_gray(im, j, height - ) ==  && Get_gray(im, j + , height - )) //左1=0，右1，右2，下1=0，右下1
             {
                 del.push_back(Point(j, height - ));
             }
             if (Get_gray(im, j - , height - ) && Get_gray(im, j + , height - ) ==  && Get_gray(im, j, height - ) ==  && Get_gray(im, j - , height - ))//左1，右1=0，下1=0，左下1
             {
                 del.push_back(Point(j, height - ));
             }
         }
     }

     for (int i = ; i < del.size(); i++)
     {
         uchar* data = image.ptr<uchar>(del[i].y);
         data[del[i].x] = ;
     }
 }