• OpenCV 闭合轮廓检测


     

     

    这个好像是骨头什么的,但是要求轮廓闭合,于是对图片进行一下膨胀操作,再次检测轮廓就好了。

     

    // A closed contour.cpp : 定义控制台应用程序的入口点。
    //
    
    #include "stdafx.h"
    
    
    // FindRotation-angle.cpp : 定义控制台应用程序的入口点。
    //
    
    // findContours.cpp : 定义控制台应用程序的入口点。
    //
    
    #include "stdafx.h"
    
    
    
    #include <iostream>
    #include <vector>
    #include <opencv2/opencv.hpp> 
    #include <opencv2/core/core.hpp>
    #include <opencv2/imgproc/imgproc.hpp>
    #include <opencv2/highgui/highgui.hpp>
    //#include "highlight"
    //#include "highgui.h"
    
    
    #pragma comment(lib,"opencv_core2410d.lib")        
    #pragma comment(lib,"opencv_highgui2410d.lib")        
    #pragma comment(lib,"opencv_imgproc2410d.lib")  
    
    #define PI 3.1415926
    
    using namespace std;
    using namespace cv;
    
    int main()
    {
    	// Read input binary image
    
    	char *image_name = "test.bmp";
    	cv::Mat image = cv::imread(image_name);
    	if (!image.data)
    		return 0; 
    
    	
    
    
    	
    	// 从文件中加载原图  
    	  // IplImage *pSrcImage = cvLoadImage(image_name, CV_LOAD_IMAGE_UNCHANGED);  
    	  Mat gray(image.size(),CV_8U);
    		  
    	  cvtColor(image,gray,CV_BGR2GRAY); 
    		 // 转为2值图
    	 threshold(gray,gray,145,255,cv::THRESH_BINARY_INV);
    	//cvThreshold(pSrcImage,pSrcImage,145,255,cv::THRESH_BINARY_INV);
    		   
    	
    	   //image = gray;
    
    	   cv::namedWindow("Binary Image");
    	   cv::imshow("Binary Image",gray);
    
    
    
    	   cv::Mat element(5,5,CV_8U,cv::Scalar(255));
    
    	   cv::dilate(gray,gray,element);
    	   //cv::erode(image,image,element);
    
    	   cv::namedWindow("dilate Image");
    	   cv::imshow("dilate Image",gray);
    
    
    	// Get the contours of the connected components
    	std::vector<std::vector<cv::Point>> contours;
    
    	cv::findContours(gray, 
    		contours, // a vector of contours 
    		CV_RETR_EXTERNAL , // retrieve the external contours
    		CV_CHAIN_APPROX_NONE); // retrieve all pixels of each contours
    
    	// Print contours' length
    	std::cout << "Contours: " << contours.size() << std::endl;
    	std::vector<std::vector<cv::Point>>::const_iterator itContours= contours.begin();
    	for ( ; itContours!=contours.end(); ++itContours) 
    	{
    
    		std::cout << "Size: " << itContours->size() << std::endl;
    	}
    
    	// draw black contours on white image
    	cv::Mat result(image.size(),CV_8U,cv::Scalar(255));
    	cv::drawContours(result,contours,
    		-1, // draw all contours
    		cv::Scalar(0), // in black
    		2); // with a thickness of 2
    
    	cv::namedWindow("Contours");
    	cv::imshow("Contours",result);
    
    
    
    
    
    
    
    
    
    	// Eliminate too short or too long contours
    
    	/*
    	int cmin= 100;  // minimum contour length
    	int cmax= 1000; // maximum contour length
    	std::vector<std::vector<cv::Point>>::const_iterator itc= contours.begin();
    	while (itc!=contours.end()) {
    
    		if (itc->size() < cmin || itc->size() > cmax)
    			itc= contours.erase(itc);
    		else 
    			++itc;
    	}
    	
    	*/
    
    	// draw contours on the original image
    	cv::Mat original= cv::imread(image_name);
    	cv::drawContours(original,contours,
    		-1, // draw all contours
    		cv::Scalar(255,255,0), // in white
    		2); // with a thickness of 2
    
    	cv::namedWindow("Contours on Animals");
    	cv::imshow("Contours on Animals",original);
    
    	
    
    	// Let's now draw black contours on white image
    	result.setTo(cv::Scalar(255));
    	cv::drawContours(result,contours,
    		-1, // draw all contours
    		cv::Scalar(0), // in black
    		1); // with a thickness of 1
    	image= cv::imread("binary.bmp",0);
    
    	// testing the bounding box 
    	
    
    
    	
    
    	std::vector<std::vector<cv::Point>>::const_iterator itc_rec= contours.begin();
    	while (itc_rec!=contours.end())
    	{
    		cv::Rect r0= cv::boundingRect(cv::Mat(*(itc_rec)));
    		cv::rectangle(result,r0,cv::Scalar(0),2);
    			++itc_rec;
    	}
    
    	/*
    	// testing the enclosing circle 
    	float radius;
    	cv::Point2f center;
    	cv::minEnclosingCircle(cv::Mat(contours[1]),center,radius);
    	cv::circle(result,cv::Point(center),static_cast<int>(radius),cv::Scalar(0),2);
    
    	//	cv::RotatedRect rrect= cv::fitEllipse(cv::Mat(contours[1]));
    	//	cv::ellipse(result,rrect,cv::Scalar(0),2);
    
    	// testing the approximate polygon
    	std::vector<cv::Point> poly;
    	cv::approxPolyDP(cv::Mat(contours[2]),poly,5,true);
    
    	std::cout << "Polygon size: " << poly.size() << std::endl;
    
    	// Iterate over each segment and draw it
    	std::vector<cv::Point>::const_iterator itp= poly.begin();
    	while (itp!=(poly.end()-1)) {
    		cv::line(result,*itp,*(itp+1),cv::Scalar(0),2);
    		++itp;
    	}
    	// last point linked to first point
    	cv::line(result,*(poly.begin()),*(poly.end()-1),cv::Scalar(20),2);
    
    	// testing the convex hull
    	std::vector<cv::Point> hull;
    	cv::convexHull(cv::Mat(contours[3]),hull);
    
    	// Iterate over each segment and draw it
    	std::vector<cv::Point>::const_iterator it= hull.begin();
    	while (it!=(hull.end()-1)) {
    		cv::line(result,*it,*(it+1),cv::Scalar(0),2);
    		++it;
    	}
    	// last point linked to first point
    	cv::line(result,*(hull.begin()),*(hull.end()-1),cv::Scalar(20),2);
    
    	// testing the moments
    
    	// iterate over all contours
    	itc= contours.begin();
    	while (itc!=contours.end()) {
    
    		// compute all moments
    		cv::Moments mom= cv::moments(cv::Mat(*itc++));
    
    		// draw mass center
    		cv::circle(result,
    			// position of mass center converted to integer
    			cv::Point(mom.m10/mom.m00,mom.m01/mom.m00),
    			2,cv::Scalar(0),2); // draw black dot
    	}
    
    	*/
    
    	cv::namedWindow("Some Shape descriptors");
    	cv::imshow("Some Shape descriptors",result);
    
    
    	cv::waitKey();
    	return 0;
    
    
    }


     

    实现效果:

     

  • 相关阅读:
    HDU4372 Count the Buildings
    Luogu4292 WC2010重建计划
    「学习笔记」二项式反演
    微软Power BI 每月功能更新系列——10月Power BI 新功能学习
    用Synoptic Panel自定义基于图形的可视化控件--制作一张剧场售票统计报表
    Power BI十大视觉效果,知多少?
    微软Power BI 每月功能更新系列——3月Power BI 新功能学习
    微软Power BI 每月功能更新系列——4月Power BI 新功能学习
    如何使用DAX函数解决动态图表标题
    如何用DAX实现查看每个月中不同类别排名前一位,以及一个简单的svg案例
  • 原文地址:https://www.cnblogs.com/wuyida/p/6301394.html
Copyright © 2020-2023  润新知