#include "opencv2/imgproc/imgproc.hpp"#include "opencv2/highgui/highgui.hpp"#include <stdio.h> using namespace cv;using namespace std; vector<vector<Point> > contours; //輪廓數(shù)組vector<Point2d> centers; //輪廓質(zhì)心坐標(biāo) vector<vector<Point> >::iterator itr; //輪廓迭代器vector<Point2d>::iterator itrc; //質(zhì)心坐標(biāo)迭代器vector<vector<Point> > con; //當(dāng)前輪廓  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);  /*繪制連通區(qū)域輪廓，計算質(zhì)心坐標(biāo)*/ 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); //其它面積藍(lán)色繪制 con.pop_back();  //計算質(zhì)心 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);  /*漫水填充連通區(qū)域*/ 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 ;}