您现在的位置是：首页 > 文章详情

最大轮廓和投影

日期：2015-07-15点击：419收藏

最大轮廓和投影

最近非常多的用到了最大轮廓和投影运算。回想起来，这两种算法的确是属于非常常见的基础算法。这里加以总结和提取。

最大轮廓：

前提是图像经过了灰度和阈值处理，也可以直接处理canny的结果，有些时候需要预先经过色彩域的转换。最后得到的结果，应该是一个contour，当然可以采用一定的方法处理得到外接矩形。

 //寻找最大的轮廓 
vector <cv : :Point > FindBigestContour(Mat src) 
{     
     int imax  =  0;  //代表最大轮廓的序号 
     int imaxcontour  =  - 1;  //代表最大轮廓的大小 
    std : :vector <std : :vector <cv : :Point >>contours;     
    findContours(src,contours,CV_RETR_LIST,CV_CHAIN_APPROX_SIMPLE); 
     for ( int i = 0;i <contours.size();i ++) 
    { 
         int itmp  =  contourArea(contours[i]); //这里采用的是轮廓大小 
         if (imaxcontour  < itmp ) 
        { 
            imax  = i; 
            imaxcontour  = itmp; 
        } 
    } 
     return contours[imax]; 
} 

调用方法,得到外接矩阵

Rect boundRect = boundingRect(Mat(FindBigestContour(canny)));

投影处理 ：

这里的前提和最大轮廓是非常相似的。投影主要关心的是通过投影图像，获得原始图像中的ROI，或者获得有多少个上波形多少个下波形这些定量的结果。

 
    Vector < int > vectorV;  //横向循环 
    Vector < int > vectorH;  //纵向循环 
    Vector < int > VUpper; 
    Vector < int > VDown; 
    vector < int > HUpper; 
    vector < int > HDower; 
 
    // 做横向循环 
     for ( int i = 0;i <ostu.cols;i ++) 
    { 
        Mat data  = ostu.col(i); 
         int itmp  = countNonZero(data); 
        vectorV.push_back(itmp); 
    } 
     //上波形为VUpper,下波形为VDown 
     for ( int i = 1;i <vectorV.size();i ++) 
    { 
         if (vectorV[i - 1]  ==  0  && vectorV[i] > 0) 
        { 
            VUpper.push_back(i); 
        } 
         if (vectorV[i - 1] > 0  && vectorV[i]  ==  0) 
        { 
            VDown.push_back(i); 
        } 
    } 
    //做纵向循环，这个往往处理的是横向循环的结果图片 
       for ( int j = 0;j <ostu.rows;j ++) 
        { 
            Mat data  = roitmp.row(j); 
             int itmp  = countNonZero(data); 
            vectorH.push_back(itmp); 
        } 
         for ( int j = 0;j <vectorH.size() - 1;j ++) 
        {    
             if (vectorH[j] > 0  && vectorH[j + 1]  ==  0) 
            { 
                HDower.push_back(j); 
            } 
             if (vectorH[j]  ==  0  && vectorH[j + 1] > 0) 
            { 
                HUpper.push_back(j); 
            } 
        } 
 
      //由于处理的上波形和下波形可能会有问题，需要进行一定的处理 
      //这里的一个处理就是提出哪些短暂的空白区域 
     for ( int j = 0;j <HDower.size() - 1;j ++) 
            { 
                 //得出之间空白的区域 
                 int iwidth  = HUpper[j + 1]  - HDower[j]; 
                 if (iwidth  >  10) 
                { 
                    iresult  = iresult  + 1; 
                } 
            } 