OpenCV aprendiendo operaciones básicas de imágenes (12): transformación de Hough

Transformación de Hough

Coordenadas polares

x=\rho cos\Theta \ \ ,\ \ y=\rho sin\Theta \\ \rho^2 = x^2+y^2, tan\Theta = \frac {y}{x}(x\neq 0)

Detección de línea 

Los pasos del proceso:

  • Transforme puntos en coordenadas espaciales en el espacio de coordenadas polares (espacio Hough) a través de coordenadas polares
  • Las líneas representadas por los puntos de intersección en el espacio de Hough son las líneas donde se ubican estos puntos
  • Luego use coordenadas polares para transformar inversamente los puntos seleccionados por Hough en el espacio de Hough en líneas rectas en coordenadas rectangulares.

y = ax + b \\ y = (- \ frac {cos \ Theta} {sin \ Theta}) x + (\ frac {r} {sim \ Theta}) \\ r_ \ Theta = x_0 \ cdot cos \ Theta + y_0 \ cdot sin \ Theta

Detección de círculos

Los pasos del proceso

  •  Realice primero un filtrado medio en la imagen, porque la detección del círculo de Hough es más sensible al ruido de sal y pimienta en la imagen
  • Realice una detección de bordes astuta, primero filtre los puntos que pueden ser bordes
  • Utilice ecuaciones paramétricas para convertir el sistema de coordenadas cartesianas en un sistema de coordenadas polares

              \ frac {(xa) ^ 2} {R} + \ frac {(yb) ^ 2} {R} = 1 \\ \ left \ {\ begin {matrix} x = a + Rcos \ Theta \\ y = b + Rsin \ Theta \ end {matriz} \ right.

  • Suponiendo que cada píxel es el centro (a, b), después de determinar (a, b), cambie θ y R para dibujar varios círculos, y cuando estos círculos se cruzan en un punto, es el centro de un círculo y el radio es Puede haber varios (determinado por la distancia desde esos puntos hasta este punto)

Introducción a la API

Detección de línea

HoughLines(
InputArray src,//输入图像,必须是8-bit的灰度图
OutputArray lines,//输出图像,是极坐标形式的直线
double rho,//生成极坐标时候的像素扫描步长,以像素为单位的距离步长。
double theta,//生成极坐标时候的角度扫描步长,以弧度为单位的角度步长。
int threshold,//阈值,只有获得足够交点的极坐标才会被视为是一条有效的直线
double srn=0,//默认值为0,用于在多尺度霍夫变换中作为参数rho的除数,rho=rho/srn。
double stn=0,//默认值为0,用于在多尺度霍夫变换中作为参数theta的除数,theta=theta/stn。
double min_theta,//表示角度扫描范围0~180之间,默认即可
double max_theta = CV_PI//
)
//需要自己将他反变换回直角坐标系中

HoughLinesP
(
InputArray src,//输入图像,必须是8-bit的灰度图
OutputArray lines,//输出直线,是直角坐标形式的直线
double rho,//生成极坐标时候的像素扫描步长
double theta,//生成极坐标时候的角度扫描步长,通常取CV_PI/180
double minLineLength = 0,//最小直线长
double maxLineGap = 0,//最大间隔,即直线与直线之间的间隔 
)
//直接返回直角坐标系中的直线

Detección de círculos

HoughCircles(
    InputArray src,//输入图像,必须是8-bit的灰度图
    OutputArray lines,//输出结果,发现圆的信息
    Int method,//方法 HOUGH_GRADENT 基于梯度的
    Double dp,//dp = 1 ,DIP,即尺度
    Double mindist,//两个圆弧相隔的最短距离,比这个距离小就认为是同一个圆,大于则认为是同心圆
    Double param1,//canny的弱边缘阈值
    Double param2,//中心点累加器的阈值,即有多数个像素在同一弧度上,才判断为圆,设定此值需要参考半径
    Int minradiux,//最小半径
    Int maxradiux,//最大半径
)

Código y práctica

Detección de línea


#include <iostream>
#include <math.h>
#include <opencv2/opencv.hpp>

using namespace cv;
using namespace std;

int main(int argc, char* argv[])
{
	//src = imread("src.jpg");
	Mat src = imread("1.PNG");
	Mat edge_img;
	if (!src.data)
	{
		cout << "cannot open image" << endl;
		return -1;
	}
	Mat gray_src;
	cvtColor(src, gray_src, COLOR_BGR2GRAY);
	Canny(src, edge_img, 150, 250);
	
	
	imshow("edge image", edge_img);

	vector<Vec4f> plines;
	//创建向量存储检测结果
	HoughLinesP(edge_img, plines, 1, CV_PI/180.0, 10, 0, 5);
	//检测直线
	Scalar color = Scalar(0, 0, 255);

	for (size_t i = 0; i < plines.size(); i++)
	{
		Vec4f hline = plines[i];
		line(gray_src, Point(hline[0], hline[1]), Point(hline[2], hline[3]), color, 3, LINE_AA);
	}
	imshow("OUTPUT", gray_src);
	
	waitKey(0);
	
	return 0;
}

#include <iostream>
#include <math.h>
#include <opencv2/opencv.hpp>

using namespace cv;
using namespace std;


int main(int argc, char* argv[])
{
	Mat src = imread("5.jpg");
	//Mat src = imread("1.PNG");
	Mat edge_img;
	if (!src.data)
	{
		cout << "cannot open image" << endl;
		return -1;
	}

	Mat moutput;
	medianBlur(src, moutput,3);
	cvtColor(moutput, moutput, COLOR_BGR2GRAY);


	vector<Vec3f> pcircles;
	//创建向量存储检测结果
	HoughCircles(moutput, pcircles, HOUGH_GRADIENT, 1, 15, 30, 35, 15, 50);
	//检测圆
	Scalar color = Scalar(0, 255, 0);
	for (size_t i = 0; i < pcircles.size(); i++)
	{
		Vec3f cc = pcircles[i];
		circle(src, Point(cc[0], cc[1]), cc[2], color, 2);
		circle(src, Point(cc[0], cc[1]), 2, color, 2);
	}
	imshow("OUTPUT", src);
	
	waitKey(0);
	
	return 0;
}

Supongo que te gusta

Origin blog.csdn.net/fan1102958151/article/details/107448800
Recomendado
Clasificación