本文实例为大家分享了Opencv提取连通区域轮廓的具体代码,供大家参考,具体内容如下
在进行图像分割后,可能需要对感兴趣的目标区域进行提取,比较常用的方法是计算轮廓。
通过轮廓可以获得目标的一些信息:
(1)目标位置
(2)目标大小(即面积)
(3)目标形状(轮廓矩)
当然,轮廓不一定代表希望目标区域,阈值分割时可能造成一部分信息丢失,因此可以计算轮廓的质心坐标,再进行漫水填充。
程序中有寻找质心+填充,但效果不好,因此就不放填充后的图了。
实验结果:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
|
#include "opencv2/imgproc/imgproc.hpp" #include "opencv2/highgui/highgui.hpp" #include <stdio.h> using namespace cv; using namespace std; vector<vector<Point> > contours; //轮廓数组 vector<Point2d> centers; //轮廓质心坐标 vector<vector<Point> >::iterator itr; //轮廓迭代器 vector<Point2d>::iterator itrc; //质心坐标迭代器 vector<vector<Point> > con; //当前轮廓 int main() { double area; double minarea = 100; double maxarea = 0; Moments mom; // 轮廓矩 Mat image,gray,edge,dst; namedWindow( "origin" ); namedWindow( "connected_region" ); image = imread( "view.jpg" ); cvtColor(image, gray, COLOR_BGR2GRAY); blur(gray, edge, Size(3,3)); //模糊去噪 threshold(edge,edge,200,255,THRESH_BINARY); //二值化处理 /*寻找轮廓*/ findContours( edge, contours, CV_RETR_EXTERNAL, CV_CHAIN_APPROX_SIMPLE ); itr = contours.begin(); //使用迭代器去除噪声轮廓 while (itr!=contours.end()) { area = contourArea(*itr); if (area<minarea) { itr = contours.erase(itr); //itr一旦erase,需要重新赋值 } else { itr++; } if (area>maxarea) { maxarea = area; } } dst = Mat::zeros(image.rows,image.cols,CV_8UC3); /*绘制连通区域轮廓,计算质心坐标*/ Point2d center; itr = contours.begin(); while (itr!=contours.end()) { area = contourArea(*itr); con.push_back(*itr); if (area==maxarea) drawContours(dst,con,-1,Scalar(0,0,255),2); //最大面积红色绘制 else drawContours(dst,con,-1,Scalar(255,0,0),2); //其它面积蓝色绘制 con.pop_back(); //计算质心 mom = moments(*itr); center.x = ( int )(mom.m10/mom.m00); center.y = ( int )(mom.m01/mom.m00); centers.push_back(center); itr++; } imshow( "origin" ,image); imshow( "connected_region" ,dst); waitKey(0); /*漫水填充连通区域*/ Point2d seed; int new_scalar = 0; int loDiff = 8, upDiff = 8; int connectivity = 4; itrc = centers.begin(); while (itrc!=centers.end()) { seed = *itrc; floodFill(image,seed,Scalar::all(new_scalar),NULL, Scalar::all(loDiff),Scalar::all(upDiff),connectivity); itrc++; } waitKey(0); return 0 ; } |
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持服务器之家。
原文链接:https://blog.csdn.net/mingtian715/article/details/53389604