### 霍夫变换检测直线并计算两线交点 在 C++ 和 OpenCV 中，可以利用 `cv::HoughLines` 或 `cv::HoughLinesP` 函数来检测图像中的直线。这两种方法分别返回极坐标形式的直线参数 (ρ, θ) 和端点表示的线段[^1]。 #### 使用 HoughLines 进行直线检测 `cv::HoughLines` 返回的是每条直线的 ρ 和 θ 参数。这些参数可以通过以下方式转换为笛卡尔坐标系下的斜率和截距： \[ m = -\frac{\cos(\theta)}{\sin(\theta)}, \quad b = \frac{\rho}{\sin(\theta)} \] 通过解析几何的方法，可以求解两条直线的交点 \( P(x, y) \)，其公式如下： 对于两条直线方程分别为： \[ L_1: Ax + By + C = 0 \] \[ L_2: Dx + Ey + F = 0 \] 则它们的交点可由下述矩阵运算得出： \[ x = \frac{CE-BF}{AE-BD}, \quad y = \frac{AF-CD}{BD-AE} \][^1] 以下是完整的代码实现： ```cpp #include <opencv2/opencv.hpp> #include <iostream> using namespace cv; using namespace std; Point calculateIntersection(const vector<Vec2f>& lines) { float rho1 = lines[0][0]; float theta1 = lines[0][1]; float cosTheta1 = cos(theta1); float sinTheta1 = sin(theta1); float rho2 = lines[1][0]; float theta2 = lines[1][1]; float cosTheta2 = cos(theta2); float sinTheta2 = sin(theta2); double det = cosTheta1 * sinTheta2 - cosTheta2 * sinTheta1; if (fabs(det) < 1e-6) { // Lines are parallel or coincident return Point(-1, -1); // Invalid point } double x = (sinTheta2*rho1 - sinTheta1*rho2) / det; double y = (-cosTheta2*rho1 + cosTheta1*rho2) / det; return Point((int)x, (int)y); } void drawLines(Mat& img, const vector<Vec2f>& lines) { for (size_t i = 0; i < lines.size(); ++i) { float rho = lines[i][0], theta = lines[i][1]; double a = cos(theta), b = sin(theta); double x0 = a * rho, y0 = b * rho; double length = 1000; Point pt1(cvRound(x0 + length*(-b)), cvRound(y0 + length*(a))); Point pt2(cvRound(x0 - length*(-b)), cvRound(y0 - length*(a))); line(img, pt1, pt2, Scalar(0, 0, 255), 2, LINE_AA); } } int main() { Mat src = imread("image.jpg", IMREAD_GRAYSCALE); if (src.empty()) { cout << "Could not open or find the image!" << endl; return -1; } Mat edges; Canny(src, edges, 50, 150, 3); vector<Vec2f> lines; HoughLines(edges, lines, 1, CV_PI / 180, 150); Mat dst; cvtColor(edges, dst, COLOR_GRAY2BGR); if (lines.size() >= 2) { Point intersection = calculateIntersection(vector<Vec2f>{lines[0], lines[1]}); circle(dst, intersection, 5, Scalar(0, 255, 0), -1, LINE_AA); } drawLines(dst, lines); imshow("Detected Lines and Intersection", dst); waitKey(); return 0; } ``` 此程序实现了霍夫变换检测直线的功能，并绘制了前两条直线及其交点的位置。 #### 注意事项 如果输入图片质量较差或者噪声较多，则可能会影响霍夫变换的结果精度。因此，在实际应用中通常会对原始灰度图进行高斯模糊处理以减少噪声影响[^1]。