一 Contour Finding
Contours使用 STL-style vector<> 表示,如 vector<cv::Point>, vector<cv::Point2f>。opencv中,使用函数 cv::findContours() 寻找contours, 具体函数定义如下:
void cv::findContours(cv::InputOutputArray image,
cv::OutputArrayofArrays contours,
cv::OutputArray hierarchy,
int mode,
int method,
cv::Point offset = cv::Point());
void cv::findContours(cv::InputOutputArray image,
cv::OutputArrayofArrays contours,
int mode,
int method,
cv::Point offset = cv::Point());
参数 image 为8位单通道输入图像,一般情况下,该图像可能由 cv::threshold(), cv::adaptiveThreshold 生成。当 image 由 cv::Canny() 生成时,cv::findContours() 仅当边缘图像为宽度为1的细区域图像,对于闭合边缘,可以使用内边缘或者外边缘代替边缘图像进行后续分析;但对于非闭合边缘,个人认为 cv::Canny() 生成的边缘图像不适合使用 cv::findContours() 进行查找;替代方案是使用边缘跟踪算法,将 cv::Canny() 生成边缘保存在 vector<cv::Point> 中,然后使用 Contours 相关分析进行更多分析。
参数 contours 为 vector<vector<cv::Point>>, vector<vector<cv::Point2f>>, 使用 array of arrays 结构可以同时保存多条 contours。
参数 hierarchy 为 vector<cv::Vec4i>, 每个元素对应一个 contour, 元素中4个整数表示该 contour 与其他 contour 之间的关系,具体关系如下:
其中, c 表示 contour, h 表示 hole.
opencv 如何生成 contour 的 hierarchy? 个人认为(暂时没有研究 opencv 具体实现),一种可行的思路是使用区域增长,大概流程如下:
1)将种子点设置在图像边角处,使用背景颜色进行区域增长,区域增长在前景处停止;
2)使用形态学算子对当前增长后区域进行 dilate 操作,使 c0, c1,... 层级上的 contours 包含在当前增长区域中;
3)继续使用前景对图像增长,将 h00, h01,..., h10, h11,... 等层级上的 contours 包含在当前增长区域中;
4)如此循环,即可生成 contour 的 hierarchy 结构。
要使用 vector<cv::Vec4i> 表示 以上树形结构,每个 cv::Vec4i 元素值含义如下:
1) index 0: next contour(same level);
2) index 1: privous contour(same level);
3) index 2: first child(next level down);
4) index 3: parent(next level up);
参数 mode 表示如何提取 contours, 具体如下:
1) cv::PETR_EXTERNAL: 仅提取最外层 contour, 在以上图例中,仅提取 c0;
2) cv::PETR_LIST: 提取所有 contours, 组成链表结构,在图例中为: c010->c001->c000->h01->h00->c0;
3) cv::PETR_CCOMP: 提取所有 contours, 根据类型组成 hole list 和 contour list, 在图例中为:c010->c001->c000->c0, h01->h00, c0->h01(此处将前两个链表连接起来);
4) cv::PETR_TREE: 提取所有 contours, 建立树形结构。
参数 method 对 contour 压缩,以减少数据量,包括:
1) cv::CHAIN_APPROX_NONE: 不做任何压缩;
2) cv::CHAIN_APPROX_SIMPLE: 压缩水平,垂直,对角方向上的线段;
4) cv::CHAIN_APPROX_TC89_L1/cv::CHAIN_APPROX_TC89_KCOS: 使用曲线曲率信息做更复杂的压缩。
参数 offset 将 contour 平移,在使用ROI进行 contour 提取时, 通过设置不同 offset 可以很方便的将 contour 绘制在原图上。
二 Contour Drawing
使用 cv::drawContours() 在图像上绘制 contours, 具体定义如下:
void cv::drawContours(cv::InputOutputArray image,
cv::InputArrayofArrays contours,
int contourIdx,
const cv::Scalar& color,
int thickness = 1,
int lineType = 8,
cv::InputArray hierarchy = noArray(),
int maxLevel = INT_MAX(),
cv::Point offset = cv::Point());
参数 image 为需要绘制 contours 图像。
参数 contours 为 cv::findContours() 生成。
参数 contourIdx 选择需要绘制的 contour, 如果其值为 -1, 则绘制所有 contours。
参数 color thickness lineType 分别定义线条颜色,线条宽度,线条类型(4连接 cv::LINE_4, 8连接 cv::LINE_8, 反走样cv::LINE_AA)。
参数 hierarchy maxLevel 共同控制 contours 层级。
参数 offset 表示绘制时平移值。
三 Contour Operation
1 cv::approxPolyDP() 对 contour 进行多边形近似,具体定义如下:
void cv::approxPolyDP(cv::InputArray curve, cv::OutputArray approxCurve, double epsilon, bool closed);
参数 curve 可以是 vector<cv::Point>, vector<cv::Point2f> 或者 arrays of size N*2, arrays of size N*1 with 2 channels。
参数 epsilon 表示近似精度,一般通过 contour 长度的百分比计算得出。
参数 closed 表示 contour 是否闭合。
cv::approxPolyDP() 函数使用 Douglas-Peucker approximation,基本思路如下:
1)在 contour 上寻找距离最大的两个点,将 contour 一分为二, 并连接两点构成线段S;
2) 在两个半边缘上分别寻找到到线段S上的最远点,将半边缘一分为二,连接以上四个点形成四边形;
3)继续寻找到四边形各条边上的最远点,构成多边形,知道最远点到对应边上距离小于epsilon停止。
2 double cv::arcLength(cv::InputArrray points, bool closed) 求 contour 长度。
3 double cv::contourArea(cv::InputArray points, bool oriented = false) 使用格林公式求 contour 所围成的面积,oriented = true 时返回带符号面积值,该函数对复杂区域(如自交区域)将返回错误结果。
4 cv::Rect cv::boundingRect(cv::InputArray points) 求 contour 所围成的矩形(无旋转矩形)。
5 cv::RotateRect cv::minAreaRect(cv::InputArray points) 求 contour 所围成的最小矩形(旋转矩形)。
6 void cv::minEnclosingCircle(cv::InputArray points, cv::Point2f& center, float& radius) 求 contour 所围成的最小圆形。
7 cv::RotateRect cv::fitEllipse(cv::InputArray points) 使用 contour 拟合椭圆(使用最小化代价函数)。
8 cv::fitLine() 拟合直线,具体定义如下:
void cv::fitLine(cv::InputArray points, cv::OutputArray line, int distType, double param, double reps, double aeps);
参数 points 接受二维或者三维点,拟合成二维平面直线或者三维空间直线。
参数 line 为 vec4f 或者 vec6f, 前两个(或三个)元素给出直线方向,后两个(或三个)给出直线上一个点。
参数 distType 表示距离度量范数,一般使用cv::DIST_L2, cv::DIST_L1。
参数 reps aeps 表示直线方向与直线点的计算精度,一般赋值为 .01 即可。
9 void cv::convexHull(cv::InputArray points, cv::OutputArray hull, bool clockwise = false, bool returnPoints = true) 计算 contour 凸包,凸包点集计算比较简单。
10 bool cv::isContourConvex(cv::InputArray contour) 判断 contour 凸性。
11 double cv::pointPolygonTest(cv::InputArray contour, cv::Point2f pt, bool measureDist) 判断点是否在 contour 中,如果 measureDist = true, 返回点到最近边缘距离。
参考资料 Learning OpenCV 3 Adrian Kaehler & Gary Bradski