opencv contours 遇到边界

问题：目标物和图像边界相连，怎么填充内部孔洞？

上一篇  opencv contours 的问题 提到 如果先通过 findContours() 找到轮廓，再通过 drawContours() 画出该轮廓，两者的形状不会发生改变。

就像这样：

Mat src = imread("test2.bmp", 1);
	Mat gray;
	cvtColor(src, gray, CV_BGR2GRAY);
	Mat thre = gray > 1;
	imshow("thre", thre);
	int area1 = countNonZero(thre);
	cout<<"area1="<<area1<<endl;

	vector< vector<Point> > contours;
	vector< Vec4i> hierarchy;	
	findContours(thre, contours, hierarchy, CV_RETR_LIST, CV_CHAIN_APPROX_SIMPLE);
	Mat mask1 = Mat::zeros(src.size(), CV_8UC1);
	for(int i=0; i<contours.size(); ++i) {
		drawContours(mask1, contours, i, 255, CV_FILLED, 8, hierarchy);
	}
	int area2 = countNonZero(mask1);
	cout<<"area2="<<area2<<endl;
	imshow("mask1", mask1); imwrite("mask1.bmp", mask1);

原图和结果图：

处理前后计算面积：

两者是一样的。

但是今天发现一个问题，当目标物和边界相连的时候，以上规律就不对了。结果变成这样：

计算结果：

处理后的图片左边界向右移了一个像素，计算出的面积小了 1550-1500=50.  也就是左边界的边长：76-27+1=50（通过四个顶点计算出）。

把method 由CV_CHAIN_APPROX_SIMPLE  改成 CV_CHAIN_APPROX_NONE 结果是一样的。

从contours里面每个像素点的坐标看，是从左上角（x=1,y=27）开始，往下沿着左边界到左下角的点（x=1,y=76），接着下边界，逆时针绕一圈，最后是点（x=2,y=27)。看来是忽略掉边界上的点了。

恰好我需要把孔洞填充，结果变成这样：

左边界上的目标物如图：

经过findContours() 函数提取过后，这个被保存为一个轮廓，一个完整的外轮廓，且没有内轮廓。

而我想要填充黑色部分，怎么办呢？

我曲折的进行了以下处理：

	Mat src = imread("mask4.bmp", 1);
	Mat gray;
	cvtColor(src, gray, CV_BGR2GRAY);
	Mat mask4 = gray > 1;
	imshow("mask4", mask4);

	//--提取轮廓填充，得到磨粒区域mask-------
	vector< vector<Point> > contours1;
	vector< Vec4i> hierarchy1;	
	Mat mask5 = Mat::zeros(mask4.size(), mask4.type());
	findContours(mask4, contours1, hierarchy1, CV_RETR_LIST, CV_CHAIN_APPROX_NONE);
	for(int i=0; i<contours1.size(); ++i) {
		drawContours(mask5, contours1, i, 255, CV_FILLED, 8, hierarchy1);
	}
	imshow("mask5", mask5); imwrite("mask5.bmp", mask5);

	//----以上处理后，边界处的轮廓无法填充-------
	//----反色，提取靠近边界的目标物-------
	Mat mask5_inv;
	bitwise_not(mask5, mask5_inv);  imwrite("mask5_inv.bmp", mask5_inv);
	vector< vector<Point> > contours2;
	findContours(mask5_inv, contours2, CV_RETR_EXTERNAL, CV_CHAIN_APPROX_NONE);
	Mat mask5_border = Mat::zeros(mask4.size(), CV_8UC1);
	for(int i=0; i<contours2.size(); ++i){
		int area = contourArea(contours2[i]);
		if(area<2000)
			drawContours(mask5_border, contours2, i, 255, CV_FILLED, 8);
	}
	imshow("mask5_border", mask5_border);

	//----把靠近边界的磨粒区域合并---
	Mat mask6;
	addWeighted(mask5, 1.0, mask5_border, 1.0, 0.0, mask6);
	imshow("mask6", mask6); imwrite("mask6.bmp", mask6);

把 mask5 反色后提取所有外轮廓，通过面积大小排除大片背景区域，取得靠近边界线的目标物区域mask5_border，再叠加到mask5上面。

中间图及结果图：

问题暂时解决。