边缘检测 霍夫检测

#include <opencv2/highgui/highgui.hpp>
#include <opencv2/imgproc/imgproc.hpp>
#include <opencv2/core/core.hpp>
#include <opencv2/calib3d/calib3d.hpp>
#include <iostream>
using namespace cv;
using namespace std;
#define PATH "/Users/mac/ClionProjects/opencv_test/"

/**
 Cannya检测：
 1.高斯滤波消除噪声
 2.计算梯度和方向，利用sobel算子计算x,y方向梯度，合成后有幅值和方向，方向背量化为0，45，90，135
 3.非极大值抑制，消除非边缘像素
 4.滞后阈值（双阈值）
    高过高阈值则保留为边缘像素
    低于低阈值则被排除
    处于中间的若连接至高阈值像素时保留


 霍夫检测： 过一点的直线族利用极坐标形式表示，可得过定点的直线族的r和theta存在正弦函数关系，如果把不同直线族（过不同定点）
 的交点求出，及可说明这条直线经过这些定点
 */
int main()
{

    //
    Mat src = imread(PATH"home.jpg");
    Mat dstImg, grayImg,edge;
    Mat grad_x, grad_y, absGradX, absGradY,dst;
    // 梯度有正负，8位正负数用16位有符号数
    Sobel(src, grad_x, CV_16S, 1, 0, 3, 1, 1, BORDER_DEFAULT);
    convertScaleAbs(grad_x, absGradX);
    imshow("x", absGradX);
    Sobel(src, grad_y, CV_16S, 0, 1, 3, 1, 1, BORDER_DEFAULT);
    convertScaleAbs(grad_y, absGradY);
    imshow("Y",absGradY);

    addWeighted(absGradX, 0.5, absGradY, 0.5, 0, dst);
    imshow("合成结果", dst);
    

    waitKey(0);
    return 0;

}

 

#include <opencv2/highgui/highgui.hpp>
#include <opencv2/imgproc/imgproc.hpp>
#include <opencv2/core/core.hpp>
#include <opencv2/calib3d/calib3d.hpp>
#include <iostream>
using namespace cv;
using namespace std;
#define PATH "/Users/mac/ClionProjects/opencv_test/"

/**
 Cannya检测：
 1.高斯滤波消除噪声
 2.计算梯度和方向，利用sobel算子计算x,y方向梯度，合成后有幅值和方向，方向背量化为0，45，90，135
 3.非极大值抑制，消除非边缘像素
 4.滞后阈值（双阈值）
    高过高阈值则保留为边缘像素
    低于低阈值则被排除
    处于中间的若连接至高阈值像素时保留


 霍夫检测： 过一点的直线族利用极坐标形式表示，可得过定点的直线族的r和theta存在正弦函数关系，如果把不同直线族（过不同定点）
 的交点求出，及可说明这条直线经过这些定点
 */
int main()
{

    //
    Mat srcImg = imread(PATH"home.jpg");
    Mat midImg, dstImg, colCanny;
    cvtColor(srcImg, midImg, CV_BGR2GRAY);
    Canny(midImg, dstImg, 50, 200, 3);
    // 检测结果作为掩码从而得到彩色的边缘
    srcImg.copyTo(colCanny, dstImg);
    imshow("边缘图", colCanny);

    vector<Vec2f> lines; // 存放检测结果
    HoughLines(dstImg, lines, 1, CV_PI/180, 300);
    // 在图中绘制结果
    for(size_t i=0; i<lines.size();i++){
        float rho = lines[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(dstImg, pt1, pt2, Scalar(0, 0, 255), 1, LINE_AA);


    }
    imshow("ans", dstImg);
    waitKey(0);
    return 0;

}