【blog算法原理】Opencv中直线的表示方法

#include "stdafx.h"
#include <iostream>
#include "opencv2/core/core.hpp"
#include "opencv2/highgui/highgui.hpp"
#include "opencv2/imgproc/imgproc.hpp"
//【blog算法原理】Opencv中直线的表示方法
// jsxyhelu 2016年1月20日
using namespace std;
using namespace cv;
void main()
{
    Mat    src;
    Mat board;//用于将识别出来直线绘制出来
    Mat board2;
    vector<Vec4i> lines;
    vector<Vec2f> linesf;
    ///////////////////////主要流程///////////////////////////////////////
    src = imread("PureLine.jpg",0);
    board = Mat::zeros(src.size(),src.type());
    board2 = Mat::zeros(src.size(),src.type());
    ////HoughLineP测试
    HoughLinesP(src, lines, 1, CV_PI/180, 50, 50, 10 );
    ////HoughLine测试
    HoughLines(src, linesf, 1, CV_PI/180, 100, 0, 0 );
    ///////////////////////显示结果///////////////////////////////////////
    for( size_t i = 0; i < lines.size(); i++ ){
        Vec4i l = lines[i];
        line( board, Point(l[0], l[1]), Point(l[2], l[3]), Scalar(255), 1, CV_AA);
    }
    for( size_t i = 0; i < linesf.size(); i++ ){
        float rho = linesf[i][0], theta = lines[i][1];
        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*(a));
        pt2.x = cvRound(x0 - 1000*(-b));
        pt2.y = cvRound(y0 - 1000*(a));
        line( board2, pt1, pt2, Scalar(255),1, CV_AA);
    }
    imshow("src",src);
    waitKey();
}