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;
                }
            }