OpenCV利用霍夫變換進(jìn)行直線檢測
本文實(shí)例為大家分享了OpenCV利用霍夫變換進(jìn)行直線檢測的具體代碼,供大家參考,具體內(nèi)容如下
1.最簡單的霍夫變換是在圖像中識(shí)別直線。在平面直角坐標(biāo)系(x-y)中,一條直線可以用下式表示:y=kx+b。
這表示參數(shù)平面(k-b)中的一條直線。因此,圖像中的一個(gè)點(diǎn)對(duì)應(yīng)參數(shù)平面中的一條直線,圖像中的一條直線對(duì)應(yīng)參數(shù)平面中的一個(gè)點(diǎn)。對(duì)圖像上所有的點(diǎn)作霍夫變換,最終所要檢測的直線對(duì)應(yīng)的一定是參數(shù)平面中直線相交最多的那個(gè)點(diǎn)。這樣就在圖像中檢測出了直線。在實(shí)際應(yīng)用中,直線通常采用參數(shù)方程:p=x\cos\theta+y\sin\theta。
Opencv里有以下函數(shù)檢測直線(最基本的霍夫變換):
void HoughLines(InputArray image, OutputArray lines, double rho, double theta, int threshold, double srn=0, double stn=0 )
具體用法看代碼就知道了:(現(xiàn)在版本的OpenCV使用函數(shù)cvHoughLines2)
#include "opencv2/opencv.hpp" #define PI 3.1415926 int main(int argc, char *argv[]) { cv::Mat image = cv::imread ("road.jpg"); cv::Mat result; cv::cvtColor (image,result,CV_BGRA2GRAY); cv::Mat contours; cv::Canny (result,contours,125,350); //邊緣檢測 std::vector<cv::Vec2f> lines; /*霍夫變換,獲得一組極坐標(biāo)參數(shù)(rho,theta),每一對(duì)對(duì)應(yīng)一條直線,保存到lines 第3,4個(gè)參數(shù)表示在(rho,theta)坐標(biāo)系里橫縱坐標(biāo)的最小單位,即步長*/ cv::HoughLines (contours,lines,1,PI/180,80); std::vector<cv::Vec2f>::const_iterator iter = lines.begin (); std::cout<<lines.size ()<<std::endl; while(iter != lines.end()) { float rho = (*iter)[0]; float theta = (*iter)[1]; if(theta<PI/4.||theta>3.*PI/4) { //畫交點(diǎn)在上下兩邊的直線 cv::Point pt1(rho/cos(theta),0); cv::Point pt2((rho-result.rows*sin(theta))/cos(theta),result.rows); cv::line(image,pt1,pt2,cv::Scalar(255),1); } else { //畫交點(diǎn)在左右兩邊的直線 cv::Point pt1(0,rho/sin(theta)); cv::Point pt2(result.cols,(rho-result.cols*cos(theta)/sin(theta))); cv::line(image,pt1,pt2,cv::Scalar(255),1); } ++iter; } cv::namedWindow ("hough"); cv::imshow("hough",image); cv::waitKey (0); }
測試結(jié)果如下:
2.可以看出,上面的直線檢測存在以下問題:
1)只能檢測出線段所在的直線,而不知道具體線段位置,也不知道線段長度;
2)同一直線可能檢測出多條直線;
3)偶然地也可能誤判段直線。
針對(duì)這些問題,opencv有那么一個(gè)函數(shù):(現(xiàn)在版本的OpenCV使用同一個(gè)函數(shù)cvHoughLines2)
void HoughLinesP(InputArray image, OutputArray lines, double rho, double theta, int threshold, dou-
ble minLineLength=0, double maxLineGap=0)
這個(gè)方法是通過概率霍夫變換實(shí)現(xiàn)的:
1)隨機(jī)獲取邊緣圖片上的前景點(diǎn),映射到級(jí)坐標(biāo)系畫曲線;
2)當(dāng)極坐標(biāo)系里面有交點(diǎn)達(dá)到最小投票數(shù),將該點(diǎn)對(duì)應(yīng)x-y坐標(biāo)系的直線L找出來;
3)搜索邊緣圖片上前景點(diǎn),在直線L上的點(diǎn)(且點(diǎn)與點(diǎn)之間距離小于maxLineGap的)連成線段,然后這些點(diǎn)全部刪除,并且記錄該線段的參數(shù),就是起始點(diǎn)和終止點(diǎn)。(當(dāng)然這里線段長度要滿足最小長度的,否則就不用記錄了)
4)重復(fù)1),2),3)
其使用方法見代碼:
#include "opencv2/opencv.hpp" #define PI 3.1415926 class LineFinder { private: std::vector<cv::Vec4i> lines; // 直線對(duì)應(yīng)的點(diǎn)參數(shù)向量 double deltaRho; //步長 double deltaTheta; int minVote; // 判斷是直線的最小投票數(shù) double minLength; // 判斷是直線的最小長度 double maxGap; // 同一條直線上點(diǎn)之間的距離容忍度 public: LineFinder() : deltaRho(1), deltaTheta(PI/180), minVote(10), minLength(0.), maxGap(0.) {} //初始化 void setAccResolution(double dRho, double dTheta) // 設(shè)置步長 { deltaRho= dRho; deltaTheta= dTheta; } void setMinVote(int minv) // 設(shè)置最小投票數(shù) { minVote= minv; } void setLineLengthAndGap(double length, double gap) // 設(shè)置最小線段長度和線段間距容忍度 { minLength= length; maxGap= gap; } std::vector<cv::Vec4i> findLines(cv::Mat& binary) //尋找線段 { lines.clear(); cv::HoughLinesP(binary,lines, deltaRho, deltaTheta, minVote,minLength, maxGap); return lines; } void drawDetectedLines(cv::Mat &image, cv::Scalar color=cv::Scalar(255,255,255)) // 畫線段 { std::vector<cv::Vec4i>::const_iterator it2=lines.begin(); while (it2!=lines.end()) { cv::Point pt1((*it2)[0],(*it2)[1]); cv::Point pt2((*it2)[2],(*it2)[3]); cv::line( image, pt1, pt2, color); ++it2; } } }; int main(int argc, char *argv[]) { cv::Mat image = cv::imread ("road.jpg"); cv::Mat result; cv::cvtColor (image,result,CV_BGRA2GRAY); cv::Mat contours; cv::Canny (result,contours,125,350); //邊緣檢測 LineFinder finder; finder.setMinVote (80); finder.setLineLengthAndGap (100,20); finder.findLines (contours); finder.drawDetectedLines (image); std::vector<cv::Vec2f> lines; cv::HoughLines (contours,lines,1,PI/180,80); std::vector<cv::Vec2f>::const_iterator iter = lines.begin (); std::cout<<lines.size ()<<std::endl; while(iter != lines.end()) { float rho = (*iter)[0]; float theta = (*iter)[1]; if(theta<PI/4.||theta>3.*PI/4) { //畫交點(diǎn)在上下兩邊的直線 cv::Point pt1(rho/cos(theta),0); cv::Point pt2((rho-result.rows*sin(theta))/cos(theta),result.rows); cv::line(image,pt1,pt2,cv::Scalar(255),1); } else { //畫交點(diǎn)在左右兩邊的直線 cv::Point pt1(0,rho/sin(theta)); cv::Point pt2(result.cols,(rho-result.cols*cos(theta)/sin(theta))); cv::line(image,pt1,pt2,cv::Scalar(255),1); } ++iter; } cv::namedWindow ("hough"); cv::imshow("hough",image); cv::waitKey (0); }
測試結(jié)果如下:
以上就是本文的全部內(nèi)容,希望對(duì)大家的學(xué)習(xí)有所幫助,也希望大家多多支持我們。
上一篇:OpenCV使用鼠標(biāo)響應(yīng)裁剪圖像
欄 目:C語言
下一篇:OpenCV實(shí)現(xiàn)圖像角點(diǎn)檢測
本文標(biāo)題:OpenCV利用霍夫變換進(jìn)行直線檢測
本文地址:http://mengdiqiu.com.cn/a1/Cyuyan/549.html
您可能感興趣的文章
- 01-10利用C語言實(shí)踐OOP,以及new,delete的深入分析
- 01-10利用ace的ACE_Task等類實(shí)現(xiàn)線程池的方法詳解
- 01-10解析如何利用switch語句進(jìn)行字符統(tǒng)計(jì)
- 01-10利用C語言實(shí)現(xiàn)HashTable
- 01-10利用C++實(shí)現(xiàn)從std::string類型到bool型的轉(zhuǎn)換
- 01-10利用C++的基本算法實(shí)現(xiàn)十個(gè)數(shù)排序
- 01-10利用C++實(shí)現(xiàn)矩陣的相加/相稱/轉(zhuǎn)置/求鞍點(diǎn)
- 01-10c++回調(diào)之利用sink示例
- 01-10c++回調(diào)之利用函數(shù)指針示例
- 01-10利用c語言實(shí)現(xiàn)卷積碼編碼器示例


閱讀排行
本欄相關(guān)
- 04-02c語言函數(shù)調(diào)用后清空內(nèi)存 c語言調(diào)用
- 04-02func函數(shù)+在C語言 func函數(shù)在c語言中
- 04-02c語言的正則匹配函數(shù) c語言正則表達(dá)
- 04-02c語言用函數(shù)寫分段 用c語言表示分段
- 04-02c語言中對(duì)數(shù)函數(shù)的表達(dá)式 c語言中對(duì)
- 04-02c語言編寫函數(shù)冒泡排序 c語言冒泡排
- 04-02c語言沒有round函數(shù) round c語言
- 04-02c語言分段函數(shù)怎么求 用c語言求分段
- 04-02C語言中怎么打出三角函數(shù) c語言中怎
- 04-02c語言調(diào)用函數(shù)求fibo C語言調(diào)用函數(shù)求
隨機(jī)閱讀
- 01-10C#中split用法實(shí)例總結(jié)
- 01-11Mac OSX 打開原生自帶讀寫NTFS功能(圖文
- 01-10SublimeText編譯C開發(fā)環(huán)境設(shè)置
- 01-10delphi制作wav文件的方法
- 04-02jquery與jsp,用jquery
- 08-05織夢dedecms什么時(shí)候用欄目交叉功能?
- 08-05DEDE織夢data目錄下的sessions文件夾有什
- 01-10使用C語言求解撲克牌的順子及n個(gè)骰子
- 01-11ajax實(shí)現(xiàn)頁面的局部加載
- 08-05dedecms(織夢)副欄目數(shù)量限制代碼修改