findcontours获取轮廓后如何获得轮廓内平均灰度值

sckxxsjh

小白一枚，最近在做缺陷检测的东西，用到findcontours函数，找到的contours轮廓里已经有了很多信息但是不知道如何使用，现在想对轮廓内所有点的灰度值求平均，请问应该怎么做？

sckxxsjh

        for (int Size = 0; Size != size(twocontours); Size++)
        {
                int value = 0;
                int meanvalue = 0;
                for (int Size0 = 0; Size0 != size(twocontours[Size]); Size0++)
                {
                        value = value + crop.at<uchar>(twocontours[Size][Size0].x, twocontours[Size][Size0].y);
                        meanvalue = value / size(twocontours[Size]);
                }
                cout << "第" << Size+1 << "个轮廓内平均灰度值为：" << meanvalue << endl;
        }

自己写了一下现在能计算出平均灰度值但是要怎么样才能设置一个阈值把平均灰度值低于阈值的轮廓去掉呢求解

sckxxsjh

                if (!twocontours.empty())
        {
                vector<vector<Point> >::iterator it;
                for (it = twocontours.begin(); it != twocontours.end();)
                {
                        for (int Size = 0; Size != size(*it); Size++)
                        {
                                int value = 0;
                                int meanvalue = 0;
                                for (int Size0 = 0; Size0 != size(it[Size]); Size0++)
                                {
                                        value = value + crop.at<uchar>(it[Size][Size0].x, it[Size][Size0].y);
                                        meanvalue = value / size(it[Size]);
                                }
                                                                        if (meanvalue > 180)
                                        twocontours.erase(it);
                                        else it++;
                                cout << "第" << Size + 1 << "个轮廓内平均灰度值为：" << meanvalue << endl;
                        }
                }
        }
又写了一个我觉得没有什么问题的，但是运行出错，求教是哪里出了问题

sckxxsjh

value = value + crop.at<uchar>(it[Size][Size0].x, it[Size][Size0].y);
出错位置是这句，我调试了一下总是在数组的最后几个数的位置会出问题，不知道是哪里写错了，求大神解答啊