图像细化

 在我们进行图像处理的时候，有可能需要对图像进行细化，提取出图像的骨架信息，进行更加有效的分析。

     图像细化（Image Thinning），一般指二值图像的骨架化（Image Skeletonization） 的一种操作运算。

     所谓的细化就是经过一层层的剥离，从原来的图中去掉一些点，但仍要保持原来的形状，直到得到图像的骨架。骨架，可以理解为图象的中轴。

     好的细化算法一定要满足：

• 收敛性；
• 保证细化后细线的连通性；
• 保持原图的基本形状；
• 减少笔画相交处的畸变；
• 细化结果是原图像的中心线；
• 细化的快速性和迭代次数少；

    这里，我们对“Zhang并行快速细化算法”进行了实现（注意，该算法为并行算法，而我们在实现过程中并没有并行化处理，所以，效率并没有达到最好）。

    参考资料

• 细化算法
  
• 论文 A fast parallel algorithm for thinning digital patterns
  [cpp]  view plain  copy

   

  1. #include <opencv2/opencv.hpp>  
  2. #include <opencv2/core/core.hpp>  
  3. #include <iostream>  
  4. #include <vector>  
  5.   
  6.   
  7. /** 
  8.  * @brief 对输入图像进行细化 
  9.  * @param src为输入图像,用cvThreshold函数处理过的8位灰度图像格式，元素中只有0与1,1代表有元素，0代表为空白 
  10.  * @param maxIterations限制迭代次数，如果不进行限制，默认为-1，代表不限制迭代次数，直到获得最终结果 
  11.  * @return 为对src细化后的输出图像,格式与src格式相同，元素中只有0与1,1代表有元素，0代表为空白 
  12.  */  
  13. cv::Mat thinImage(const cv::Mat & src, const int maxIterations = -1)  
  14. {  
  15.     assert(src.type() == CV_8UC1);  
  16.     cv::Mat dst;  
  17.     int width  = src.cols;  
  18.     int height = src.rows;  
  19.     src.copyTo(dst);  
  20.     int count = 0;  //记录迭代次数  
  21.     while (true)  
  22.     {  
  23.         count++;  
  24.         if (maxIterations != -1 && count > maxIterations) //限制次数并且迭代次数到达  
  25.             break;  
  26.         std::vector<uchar *> mFlag; //用于标记需要删除的点  
  27.         //对点标记  
  28.         for (int i = 0; i < height ;++i)  
  29.         {  
  30.             uchar * p = dst.ptr<uchar>(i);  
  31.             for (int j = 0; j < width; ++j)  
  32.             {  
  33.                 //如果满足四个条件，进行标记  
  34.                 //  p9 p2 p3  
  35.                 //  p8 p1 p4  
  36.                 //  p7 p6 p5  
  37.                 uchar p1 = p[j];  
  38.                 if (p1 != 1) continue;  
  39.                 uchar p4 = (j == width - 1) ? 0 : *(p + j + 1);  
  40.                 uchar p8 = (j == 0) ? 0 : *(p + j - 1);  
  41.                 uchar p2 = (i == 0) ? 0 : *(p - dst.step + j);  
  42.                 uchar p3 = (i == 0 || j == width - 1) ? 0 : *(p - dst.step + j + 1);  
  43.                 uchar p9 = (i == 0 || j == 0) ? 0 : *(p - dst.step + j - 1);  
  44.                 uchar p6 = (i == height - 1) ? 0 : *(p + dst.step + j);  
  45.                 uchar p5 = (i == height - 1 || j == width - 1) ? 0 : *(p + dst.step + j + 1);  
  46.                 uchar p7 = (i == height - 1 || j == 0) ? 0 : *(p + dst.step + j - 1);  
  47.                 if ((p2 + p3 + p4 + p5 + p6 + p7 + p8 + p9) >= 2 && (p2 + p3 + p4 + p5 + p6 + p7 + p8 + p9) <= 6)  
  48.                 {  
  49.                     int ap = 0;  
  50.                     if (p2 == 0 && p3 == 1) ++ap;  
  51.                     if (p3 == 0 && p4 == 1) ++ap;  
  52.                     if (p4 == 0 && p5 == 1) ++ap;  
  53.                     if (p5 == 0 && p6 == 1) ++ap;  
  54.                     if (p6 == 0 && p7 == 1) ++ap;  
  55.                     if (p7 == 0 && p8 == 1) ++ap;  
  56.                     if (p8 == 0 && p9 == 1) ++ap;  
  57.                     if (p9 == 0 && p2 == 1) ++ap;  
  58.   
  59.                     if (ap == 1 && p2 * p4 * p6 == 0 && p4 * p6 * p8 == 0)  
  60.                     {  
  61.                         //标记  
  62.                         mFlag.push_back(p+j);  
  63.                     }  
  64.                 }  
  65.             }  
  66.         }  
  67.   
  68.         //将标记的点删除  
  69.         for (std::vector<uchar *>::iterator i = mFlag.begin(); i != mFlag.end(); ++i)  
  70.         {  
  71.             **i = 0;  
  72.         }  
  73.   
  74.         //直到没有点满足，算法结束  
  75.         if (mFlag.empty())  
  76.         {  
  77.             break;  
  78.         }  
  79.         else  
  80.         {  
  81.             mFlag.clear();//将mFlag清空  
  82.         }  
  83.   
  84.         //对点标记  
  85.         for (int i = 0; i < height; ++i)  
  86.         {  
  87.             uchar * p = dst.ptr<uchar>(i);  
  88.             for (int j = 0; j < width; ++j)  
  89.             {  
  90.                 //如果满足四个条件，进行标记  
  91.                 //  p9 p2 p3  
  92.                 //  p8 p1 p4  
  93.                 //  p7 p6 p5  
  94.                 uchar p1 = p[j];  
  95.                 if (p1 != 1) continue;  
  96.                 uchar p4 = (j == width - 1) ? 0 : *(p + j + 1);  
  97.                 uchar p8 = (j == 0) ? 0 : *(p + j - 1);  
  98.                 uchar p2 = (i == 0) ? 0 : *(p - dst.step + j);  
  99.                 uchar p3 = (i == 0 || j == width - 1) ? 0 : *(p - dst.step + j + 1);  
  100.                 uchar p9 = (i == 0 || j == 0) ? 0 : *(p - dst.step + j - 1);  
  101.                 uchar p6 = (i == height - 1) ? 0 : *(p + dst.step + j);  
  102.                 uchar p5 = (i == height - 1 || j == width - 1) ? 0 : *(p + dst.step + j + 1);  
  103.                 uchar p7 = (i == height - 1 || j == 0) ? 0 : *(p + dst.step + j - 1);  
  104.   
  105.                 if ((p2 + p3 + p4 + p5 + p6 + p7 + p8 + p9) >= 2 && (p2 + p3 + p4 + p5 + p6 + p7 + p8 + p9) <= 6)  
  106.                 {  
  107.                     int ap = 0;  
  108.                     if (p2 == 0 && p3 == 1) ++ap;  
  109.                     if (p3 == 0 && p4 == 1) ++ap;  
  110.                     if (p4 == 0 && p5 == 1) ++ap;  
  111.                     if (p5 == 0 && p6 == 1) ++ap;  
  112.                     if (p6 == 0 && p7 == 1) ++ap;  
  113.                     if (p7 == 0 && p8 == 1) ++ap;  
  114.                     if (p8 == 0 && p9 == 1) ++ap;  
  115.                     if (p9 == 0 && p2 == 1) ++ap;  
  116.   
  117.                     if (ap == 1 && p2 * p4 * p8 == 0 && p2 * p6 * p8 == 0)  
  118.                     {  
  119.                         //标记  
  120.                         mFlag.push_back(p+j);  
  121.                     }  
  122.                 }  
  123.             }  
  124.         }  
  125.   
  126.         //将标记的点删除  
  127.         for (std::vector<uchar *>::iterator i = mFlag.begin(); i != mFlag.end(); ++i)  
  128.         {  
  129.             **i = 0;  
  130.         }  
  131.   
  132.         //直到没有点满足，算法结束  
  133.         if (mFlag.empty())  
  134.         {  
  135.             break;  
  136.         }  
  137.         else  
  138.         {  
  139.             mFlag.clear();//将mFlag清空  
  140.         }  
  141.     }  
  142.     return dst;  
  143. }  
  144.   
  145.   
  146. int main(int argc, char*argv[])  
  147. {  
  148.     //获取图像  
  149.     if (argc != 2)  
  150.     {  
  151.         std::cout << "参数个数错误！" << std::endl;  
  152.         return -1;  
  153.     }  
  154.     cv::Mat src = cv::imread(argv[1], cv::IMREAD_GRAYSCALE);  
  155.     if (src.empty())  
  156.     {  
  157.         std::cout << "读取文件失败！" << std::endl;  
  158.         return -1;  
  159.     }  
  160.   
  161.     //将原图像转换为二值图像  
  162.     cv::threshold(src, src, 128, 1, cv::THRESH_BINARY);  
  163.     //图像细化  
  164.     cv::Mat dst = thinImage(src);  
  165.     //显示图像  
  166.     dst = dst * 255;  
  167.     cv::namedWindow("src1", CV_WINDOW_AUTOSIZE);  
  168.     cv::namedWindow("dst1", CV_WINDOW_AUTOSIZE);  
  169.     cv::imshow("src1", src);  
  170.     cv::imshow("dst1", dst);  
  171.     cv::waitKey(0);  
  172. }  

  运行效果

  1原图像

  
  

  2.运行效果

  
  

•