我和opencv 8 最大类间方差法

最近的一个月一直在根据项目的具体问题来学习一些新的算法，但是在写最大类间方差的时候由于调用了很多的mat类逻辑运算和求和运算来求最大方差，导致程序的运行时间直接跑到了4秒开外，所以最近准备学习我一开始就回避的向量方法。
借用同学的一句话：没有用到指针和内存分配，学C++干嘛？
先上一个简单的向量遍历的方法（灰度图像），并累计每个像素值的个数：
有点类似与calcHist()函数

bool fiter(Mat src){    int gray;//灰度值    int histogram[256]={0};    for(int i=0;i<src.rows;i++)    {        uchar* p=src.ptr<uchar>(i);//给p指向行的首地址        //printf("A[%d][j]=%d\n",i,(*p));        for(int j=0;j<src.cols;j++)        {            gray=*(p+j);            histogram[gray]++;        }    }    for(int k=0;k<256;k++)    printf("histogram[%d]= %d\n",k,histogram[k]);    return 1;}

最大类间方差

//这段代码耗时太长，不建议采纳    //for(int T=0;T<256;T++)    //{    //  U=(g_filter<T)/255;    //  //cout<<sum(U)[0]/N<<endl;    //  w0=sum(U)[0];//背景的个数    //  w1=N-w0;//N为像素总个数    //  u0=sum(g_filter.mul(U))[0]/w0;    //  u1=sum(((U==0)/255).mul(g_filter))[0]/w1;    //  w0=w0/N;    //  w1=1-w0;    //  S=w1*w0*(u0-u1)*(u0-u1);    //  if(S>Smax)    //      Smax=S;    //}

下面用指针的方法可大大节约程序运行时间

int histogram[256]={0};    //int gray;    double S,Smax=0;    int t_OTSU;    for (int i=0;i<g_filter.rows;i++)        {            uchar* p=g_filter.ptr<uchar>(i);            //uchar* p=src.ptr<uchar>(i);            for(int j=0;j<g_filter.cols;j++)                //gray=*(p+j);                histogram[*(p+j)]++;        }    float w0=0,w1=0,u0=0,u1=0;//背景和前景的比例及灰度均值    int N=g_filter.rows*g_filter.cols;    for(int T=50;T<150;T++)//遍历阈值        {        for(int k=0;k<256;k++)            {                //printf("histogram[%d]= %d\n",k,histogram[k]);                if(k<T)                 {                    w0+=histogram[k];//背景的个数                    u0+=histogram[k]*k;                }                else                 {                    w1+=histogram[k];                    u1+=histogram[k]*k;                }            }            u0=u0/w0;//灰度均值            w0=w0/N;//比例            u1=u1/w1;            w1=w1/N;        //  double S,Smax=0;            S=w1*w0*(u0-u1)*(u0-u1);            if(S>Smax)            {Smax=S;t_OTSU=T;}        }    cout<<"最大类间方差阈值为 "<<t_OTSU<<endl;    cout<<"最大类间方差为 "<<Smax<<endl;