OpenCV学习笔记9 （Hough线变换）

最新推荐文章于 2020-02-17 22:43:31 发布
最新推荐文章于 2020-02-17 22:43:31 发布
阅读量1.1k 收藏 1
点赞数
分类专栏： OpenCV
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u014488388/article/details/51524101
版权


1.hough变换调试程序


#include "stdafx.h"
#include "cv.h"
#include "highgui.h"
#include "cxcore.h"



int _tmain(int argc, _TCHAR* argv[])
{
	IplImage* src = cvLoadImage("D:\\PERSONAL\\VC++\\2016\\毕业论文\\论文图表文件he文字资料\\模式识别\\picture\\模式1和2处理\\canny\\J11.jpg",CV_LOAD_IMAGE_GRAYSCALE);
	IplImage* src3 = cvCreateImage(cvGetSize(src),IPL_DEPTH_8U,3);//三通道显示
	CvMemStorage *storage = cvCreateMemStorage();  
    CvSeq *lines = 0;  

	cvMerge(src,src,src,NULL,src3);
	cvNamedWindow("1");
	cvShowImage("1",src);
	lines = cvHoughLines2(src,storage,CV_HOUGH_STANDARD,1,CV_PI/180,100,0,0); 


	int i = 0; 
H: 
	if(i<lines->total)
	{
		float* line = (float*)cvGetSeqElem(lines,i);
		float rho = line[0];
		float theta = line[1];
		CvPoint pt1, pt2;
		double a = cos(theta), b = sin(theta);
		double x0 = a*rho, y0 = b*rho;
		pt1.x = cvRound(x0 + 1000*(-b));
		pt1.y = cvRound(y0 + 1000*(a));
		pt2.x = cvRound(x0 - 1000*(-b));
		pt2.y = cvRound(y0 - 1000*(a));
		cvLine( src3, pt1, pt2, CV_RGB(255,0,0), 1, CV_AA, 0 );
		cvCircle(src3,pt1,5,CV_RGB(0,255,0),2,8,0);
		cvCircle(src3,pt2,5,CV_RGB(0,255,0),2,8,0);
		cvShowImage("3",src3);
		printf("i = %d\n",i);
		printf("a = %f,b = %f\n",a,b);
		printf("roh = %f,theta = %f\n",rho,theta);
		printf("x0 = %f,y0 = %f\n",x0,y0);
		printf("x1 = %d,y1 = %d\n",pt1.x,pt1.y);
		printf("x2 = %d,y2 = %d\n",pt2.x,pt2.y);
		i++;
	}
        

	char c(0);
	while (1)
	{
		cvShowImage("3",src3);
		c = cvWaitKey(1);
		if(c == 13)
			goto H;
		if(c == 27)
			break;
	}

	cvDestroyWindow("1");
	cvReleaseImage(&src);
	
	return 0;
}

// if(theta>2.9 || theta<0.2) //vertical


2.hough提取水平线程序

// deal.cpp : 定义控制台应用程序的入口点。
//

#include "stdafx.h"
#include "cv.h"
#include "highgui.h"
#include "cxcore.h"



int _tmain(int argc, _TCHAR* argv[])
{
	IplImage* src = cvLoadImage("D:\\...\\H11.jpg",CV_LOAD_IMAGE_GRAYSCALE);
	IplImage* src3 = cvCreateImage(cvGetSize(src),IPL_DEPTH_8U,3);//三通道显示
	CvMemStorage *storage = cvCreateMemStorage();  
    CvSeq *lines = 0;  

	cvMerge(src,src,src,NULL,src3);
	cvNamedWindow("1");
	cvShowImage("1",src);
	lines = cvHoughLines2(src,storage,CV_HOUGH_STANDARD,1,CV_PI/180,100,0,0);
	int y_median(0);
	int temp(0);
	int n_horizon(0);

	for(int i = 0; i < MIN(lines->total,100); i++ )
	{
		float* line = (float*)cvGetSeqElem(lines,i);
		float rho = line[0];
		float theta = line[1];
		CvPoint pt1, pt2;
		if(theta>1.3 && theta<1.7)
		{
			double a = cos(theta), b = sin(theta);
			double x0 = a*rho, y0 = b*rho;
			pt1.x = cvRound(x0 + 1000*(-b));
			pt1.y = cvRound(y0 + 1000*(a));
			pt2.x = cvRound(x0 - 1000*(-b));
			pt2.y = cvRound(y0 - 1000*(a));
		//	cvLine( src3, pt1, pt2, CV_RGB(255,0,0), 1, CV_AA, 0 );
			n_horizon++;
			y_median += ( (pt1.y + pt2.y)/2 );
		/*	if(y_median!=0)
				y_median = (temp + y_median)/2;
			else
				y_median = temp;*/

		}
		else
		{
			double a = cos(theta), b = sin(theta);
			double x0 = a*rho, y0 = b*rho;
			pt1.x = cvRound(x0 + 1000*(-b));
			pt1.y = cvRound(y0 + 1000*(a));
			pt2.x = cvRound(x0 - 1000*(-b));
			pt2.y = cvRound(y0 - 1000*(a));
//			cvLine( src3, pt1, pt2, CV_RGB(255,0,0), 1, CV_AA, 0 );	
		}
		printf("i = %d\n",i);
		printf("x1 = %d,y1 = %d\n",pt1.x,pt1.y);
		printf("x2 = %d,y2 = %d\n",pt2.x,pt2.y);

	}
	y_median = y_median/n_horizon;
	cvRectangle(src3,cvPoint(0,y_median+25),cvPoint(640,y_median-25),CV_RGB(0,0,0),CV_FILLED,8,0);

    cvShowImage("3",src3);
	printf("y_median = %d\n",y_median); 
	cvSaveImage("C:\\xiu.jpg",src3);

	cvWaitKey(0);


	cvDestroyWindow("1");
	cvReleaseImage(&src);
	
	return 0;
}
//if(theta>2.9 || theta<0.2) //vertical

/**

for(int i = 0; i < MIN(lines->total,100); i++ )
    {
        float* line = (float*)cvGetSeqElem(lines,i);
        float rho = line[0];
        float theta = line[1];
        CvPoint pt1, pt2;
        double a = cos(theta), b = sin(theta);
        double x0 = a*rho, y0 = b*rho;
        pt1.x = cvRound(x0 + 1000*(-b));
        pt1.y = cvRound(y0 + 1000*(a));
        pt2.x = cvRound(x0 - 1000*(-b));
        pt2.y = cvRound(y0 - 1000*(a));
        cvLine( src3, pt1, pt2, CV_RGB(255,0,0), 1, CV_AA, 0 );
		cvShowImage("3",src3);
		printf("i = %d\n",i);
		printf("x1 = %d,y1 = %d\n",pt1.x,pt1.y);
		printf("x2 = %d,y2 = %d\n",pt2.x,pt2.y);
		
	
	}
**/



3.直线方程计算

通过测试可以知道，以上代码计算出来的直线，其实是由两个坐标值超出图像显示区域的点连线而成的。

笔者在程序中需要对直线区域进行填充，使用的是cvrectangle（fill），发现在水平和垂直线上填充效果很好，但是斜线往往有蛮大偏差，不能完全覆盖直线区域。最后发现原因就是没有搞清楚以上代码中的pt1、pt2的含义。使用pt3、pt4才是我想要的效果。



所以需要利用pt1、pt2计算出直线的方程，然后推知pt3、pt4的坐标，代码修改如下：

			if(theta>1.3 && theta<1.7)   //horizon
			{  
				double a = cos(theta), b = sin(theta);  
				double x0 = a*rho, y0 = b*rho;  
				pt_1.x = cvRound(x0 + 1000*(-b));  //88
				pt_1.y = cvRound(y0 + 1000*(a));  
				pt_2.x = cvRound(x0 - 1000*(-b));  //88
				pt_2.y = cvRound(y0 - 1000*(a));  	
				double KK;
				KK = (double)(pt_1.y - pt_2.y)/(double)(pt_1.x - pt_2.x);
				int yL = cvRound( KK*(0 - pt_2.x) + pt_2.y );
				int yR = cvRound( KK*(640 - pt_2.x) + pt_2.y );
				CvPoint pt1;CvPoint pt2;CvPoint pt3;CvPoint pt4;
				pt1.x = 0;	pt1.y = yL - 25;
				pt2.x = 0;	pt2.y = yL + 25;
				pt3.x = 640;	pt3.y = yR + 25;
				pt4.x = 640;	pt4.y = yR - 25;

				fillArbitaryRectangle(temp,pt1,pt2,pt3,pt4);
				cvLine( src_original_3, pt_1, pt_2, CV_RGB(160,32,240), 1, CV_AA, 0 ); 
			}
			if(theta>2.9 || theta<0.2) //vertical
			{  
				double a = cos(theta), b = sin(theta);  
				double x0 = a*rho, y0 = b*rho;  
				pt_1.x = cvRound(x0 + 1000*(-b));  
				pt_1.y = cvRound(y0 + 1000*(a));  //888
				pt_2.x = cvRound(x0 - 1000*(-b));  
				pt_2.y = cvRound(y0 - 1000*(a));  //888
				double KK;
				KK = (double)(pt_1.x - pt_2.x)/(double)(pt_1.y - pt_2.y);
				int xT = cvRound( KK*(0 - pt_2.y) + pt_2.x );
				int xB = cvRound( KK*(480 - pt_2.y) + pt_2.x );
				CvPoint pt1;CvPoint pt2;CvPoint pt3;CvPoint pt4;
				pt1.x = xT - 25;	pt1.y = 0;
				pt2.x = xT + 25;	pt2.y = 0;
				pt3.x = xB + 25;	pt3.y = 480;
				pt4.x = xB - 25;	pt4.y = 480;
				fillArbitaryRectangle(temp,cvPoint(0,0),cvPoint(0,480),pt4,pt1);
				fillArbitaryRectangle(temp,pt2,pt3,cvPoint(640,480),cvPoint(640,0));
				cvLine( src_original_3, pt_1, pt_2, CV_RGB(160,32,240), 1, CV_AA, 0 ); 
			}
		}

其中的填充函数fillArbitaryRectangle如下： 

//src 3 channels and will be filled
void CALEXDlg::fillArbitaryRectangle(IplImage* src,CvPoint pt1,CvPoint pt2,CvPoint pt3,CvPoint pt4)
{
	int arr[1];  
	arr[0] = 4;  
	CvPoint ** pt = new CvPoint*[1];  
	pt[0] = new CvPoint[4];  
	pt[0][0] = pt1;  
	pt[0][1] = pt2;  
	pt[0][2] = pt3;  
	pt[0][3] = pt4;  
	cvPolyLine( src_original_hough, pt, arr, 1, 1, CV_RGB(255,97,0));  
	cvFillPoly(src,pt,arr,1,CV_RGB(0,0,0)); 

}







丐帮二帮主
关注
专栏目录
OpenCV学习笔记-Hough变换怎么用
Home_Wood的博客
03-22 283
在本教程中，您将学习如何： - 使用OpenCV函数HoughLines（）和HoughLinesP（）来检测图像中的线条。 - 使用OpenCV函数HoughCircles（）来检测图像中的圆圈。
OpenCV学习笔记（十三）：霍夫变换：HoughLines(),HoughLinesP(),HoughCircles( )
Leon_Chen的博客
09-09 6705
OpenCV学习笔记（十三）：霍夫变换：HoughLines(),HoughLinesP(),HoughCircles( ) 1、霍夫线变换 OpenCV支持三种不同的霍夫线变换，它们分别是： 1）标准霍夫变换(Standard Hough Transform，SHT) 2）多尺度霍夫变换（Multi-Scale Hough Transform，MSHT） 3）累计概率霍夫变换(Progressi...
(一)OpenCV图像处理基础_19_Hough线变换
weixin_45553439的博客
02-17 134
二进制图像中的任何点都可能属于某些可能的线。 前提条件：边缘检测已经完成; 平面空间到极坐标空间转换(ρ,θ),从[0~360]空间，可以得到r的大小; 属于同一条直线上点在极坐标空(r, theta)必然在一个点上有最强的信号出现，根据此反算到平面坐标中就可以得到直线上各点的像素坐标,从而得到直线。 标准的霍夫变换 cv::HoughLines从平面坐标转换到霍夫空间，最终输出是(ρ,θ)表示...
Opencv图像处理---Hough线变换
LYKymy的博客
10-19 344
理论 霍夫线变换是用于检测直线的变换。 为了应用变换，首先需要边缘检测预处理。 它是如何工作的？ 图像空间中的一条线可以用两个变量表示。 例如： 在笛卡尔坐标系中：参数：（m，b） 在极坐标系中：参数：（r，θ） 对于霍夫变换，我们将在极坐标系统中表达线条。 因此，线方程可以写成： 通常对于每个点（x0，y0），我们可以定义经过该点的行： 如果对于给定的（x0，y0）我们绘制经过它的线...
【杂谈opencv】OpenCV中的cvRound()、cvFloor()、 cvCeil()函数讲解
热门推荐
拜小白的成长之路，告别小白
12-20 3万+
cvRound()：返回跟参数最接近的整数值，即四舍五入；cvFloor()：返回不大于参数的最大整数值，即向下取整；cvCeil()：返回不小于参数的最小整数值，即向上取整；
opencv 关于霍夫变换中pt1,pt2点设定问题
kalp_yp
07-02 993
    首先把毛星云的霍夫变换检测直线中确定pt1,pt2两点的代码片段贴出来。 Point pt1, pt2; double a = cos(theta), b = sin(theta); double x0 = a*rho, y0 = b*rho; pt1.x = cvRound(x0 + 1000 * (-b)); pt1.y = cvRound(y0 + 1000 ...
cvRound, cvFloor, cvCeil
RAMBO14541516的博客
04-07 256
函数 cvRound, cvFloor, cvCeil 用一种舍入方法将输入浮点数转换成整数。 cvRound 返回和参数最接近的整数值。 cvFloor 返回不大于参数的最大整数值。cvCeil 返回不小于参数的最小整数值。在某些体系结构中该函数 工作起来比标准 C 操作起来还要快。 注意下面的编程技巧： typedefunionCv32suf{ int i;...
【opencv学习笔记 21 直线检测】
12-21
本篇学习笔记将探讨如何使用霍夫变换（Hough Transform）来检测图像中的直线。霍夫变换是一种参数空间搜索算法，它可以将图像中的直线转换为极坐标系中的点，从而找出图像中的所有直线。 首先，霍夫直线变换的前提...
opencv教程个人学习笔记总结及C++示例.zip
最新发布
11-23
这个压缩包文件“opencv教程个人学习笔记总结及C++示例.zip”显然是一个学习OpenCV的资料集合，其中包含了一些用C++编写的示例代码，帮助初学者理解和实践OpenCV的功能。 首先，我们来看看这些示例文件名，它们揭示...
OpenCV学习笔记大集锦.rar
09-17
本资源"OpenCV学习笔记大集锦"是专为初学者设计的一份资料合集，旨在帮助新手快速入门计算机视觉，并掌握OpenCV的基本使用。 在C语言背景下学习OpenCV，首先需要理解基础的编程概念，因为OpenCV库主要用C++编写，但...
OpenCV中Hough检测直线中pt1、pt2点的确定
一步一步走
04-10 2804
学习hough变换，一般会看到pt1.x = cvRound(x0 + 1000*(-b));这段代码。一开始可能不是很好理解。查阅了资料和原理，现在写下来分享。 这里的pt1和pt2是一条直线的两个端点，这里已知下图中的rh0 和θ，现在只需要求图中“任意”两点， 使用OpenCV的cvLine函数画出pt1点到pt2的直线 。 看看下图 就明白这里1000什么的是为什么了。
opencv中的标准霍夫线变换（HoughLines）和统计霍夫变换（HoughLinesP）
longwinyang的博客
10-11 1万+
本文主要讲解标准霍夫线变换（HoughLines）和统计霍夫变换（HoughLinesP）函数，并对函数在调用中的一些细节进行些讲解。  针对标准霍夫线变换（HoughLines）和统计霍夫变换（HoughLinesP）的一些细节和原理可以参考：  http://www.opencv.org.cn/opencvdoc/2.3.2/html/doc/tutorials/imgproc/imgtr
OpenCV霍夫变换系列（前篇）-经典霍夫线变换
qingyun_wudaoletu的博客
09-08 8740
最近新来的的我的大学室友（现在也是我的学弟）在研究霍夫线变换，我之前只是知道这玩意可以拿来做直线检测，并没有深入研究，那既然提到了，还是按照我们的老规矩，原理，示例以及OpenCV这一套流程走下来
边沿检测与提取，轮廓跟踪与Hough变换(转)
萧萧雨轩－－云计算/网格计算技术/智能仿生算法&算法优化/
07-19 4633
边沿检测与提取，轮廓跟踪首先,引入了模板操作。边沿检测我们给出一个模板 和一幅图象 。不难发现原图中左边暗，右边亮，中间存在着一条明显的边界。进行模板操作后的结果如下： 。可以看出，第3、4列比其他列的灰度值高很多，人眼观察时，就能发现一条很明显的亮边，其它区域都很暗，这样就起到了边沿检测的作用。为什么会这样呢？仔细看看那个模板就明白了，它的意思是将右邻点的灰度值减左
OpenCV的霍夫变换(Hough Transform)直线检测
不眠旅行
10-20 7624
霍夫变换(Hough Transform)的主要思想：一条直线在平面直角坐标系(x-y)中可以用y=ax+b式表示，对于直线上一个确定的点(x0,y0)，总符合y0-ax0=b，而它可以表示为参数平面坐标系(a-b)中的一条直线。因此，图像中的一个点对应参数平面的一条直线，同样，图像中的一条直线对应参数平面上的一个点。 笔记包括两个测试，基本Hough变换和概率Hough变换的直线检测。
OpenCV 学习（Hough 变换提取直线)
Ivan 的专栏
12-03 2万+
OpenCV 学习（Hough 变换提取直线)在机器视觉应用中，我们经常要提取图像中的各种特征，最基本的特征就是图像中的线条、拐角等。这篇笔记就来讲讲如何提取图像中的直线。这里使用的方法叫做 Hough 变换。Hough 变换这个名称最早是在 Richard Duda 和 Peter Hart 两人于 1972 年合写的发表于 Comm. ACM 文章 《Use of the Hough Trans
opencv标准霍夫变换检测直线以及描绘每条线段中pt1、pt2点的确定
liumoude6的博客
08-29 981
#include "opencv2/core/core.hpp" #include "opencv2/imgproc/imgproc.hpp" #include "opencv2/highgui/highgui.hpp" #include "opencv2/opencv.hpp" #include using namespace cv; using namespace std; using s
OpenCV中的Hough变换：检测线与圆
"利用Hough变换原理，OpenCV实现线与圆的检测。示例代码展示了如何使用cvHoughLines2和cvHoughCircles函数进行图像处理，适用于线段和圆形的检测。" 在计算机视觉领域，图像中的几何形状识别是常见的任务之一，...
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值