常规情况与opencv图像中,计算直线与矩形框的交点

news2025/1/2 3:37:10

文章目录

  • 1、普通方式
    • 1.1、普通计算过程
    • 1.2、优化方式
  • 2、图像中的情况
    • 2.1、常规处理
    • 2.2、opencv中的处理
      • 2.2.1、cv::clipLine函数
      • 2.2.2、测试代码
      • 2.2.3、测试结果

1、普通方式

已知矩形框左上(x1,y1)、右下(x2,y2)点,直线方程 y = kx+b,求交点。

这个问题转换为:计算直线与矩形框的四条线段的交点,因此最多执行4次。
首先计算直线与线段所在直线交点,若交点的在线段坐标的范围内,那么认为是有交点。

1.1、普通计算过程

这里不考虑直线平行的常规情况,如下
{ a 1 x + b 1 y + c 1 = 0 a 2 x + b 2 y + c 2 = 0 \left\{\begin{matrix} a_1x + b_1y+c_1 = 0 \\ a_2x + b_2y+c_2 = 0 \end{matrix}\right. {a1x+b1y+c1=0a2x+b2y+c2=0
计算得到横坐标
x = b 2 c 1 − b 1 c 2 b 2 a 1 − b 1 a 2 x=\frac{b_2c_1 - b_1c_2}{b_2a_1-b_1a_2} x=b2a1b1a2b2c1b1c2

当其位于某一条线段的横坐标坐标取值范围 x ∈ [ x 1 , x 2 ] x \in [x_1, x_2] x[x1,x2] 时,认为存在交点
{ x = b 2 c 1 − b 1 c 2 b 2 a 1 − b 1 a 2 y = − a 2 x − c 2 b 2 \left\{\begin{matrix} x=\frac{b_2c_1 - b_1c_2}{b_2a_1-b_1a_2} \\ y = \frac{-a_2x-c_2}{b_2} \end{matrix}\right. {x=b2a1b1a2b2c1b1c2y=b2a2xc2

1.2、优化方式

我们分两种情况。直线垂直和水平,以及普通情况。

  • (1)垂直和水平情况
    在这里插入图片描述
    这种方式简单,当垂直时直线 x = m x = m x=m,若 m ∈ [ x 1 , x 2 ] m \in [x_1, x_2] m[x1,x2],那么交点为 ( m , y 1 ) (m,y_1) (m,y1) ( m , y 2 ) (m,y_2) (m,y2);同理,当水平式直线 y = n y = n y=n,若 n ∈ [ y 1 , y 2 ] n \in [y_1, y_2] n[y1,y2],那么交点为 ( x 1 , n ) (x_1,n) (x1,n) ( x 2 , n ) (x_2,n) (x2,n)

  • (2)不相交时
    在下图中我们绘制斜率分别为正、负的情况,且与矩形不相交的示意图,大概有8种情况下,矩形框的4个顶点均在某条直线的一侧。
    在这里插入图片描述
    因此,我们判断直线是否与矩形相交,可以直接判定矩形框四个顶点 ( x 1 , y 1 ) , ( x 2 , y 1 ) , ( x 2 , y 2 ) , ( x 1 , y 2 ) (x_1,y_1),(x_2,y_1),(x_2,y_2),(x_1,y_2) (x1,y1),(x2,y1),(x2,y2),(x1,y2) 以顺时针带入直线方程函数 f ( x , y ) = a x + b y + c f(x,y) = ax+by+c f(x,y)=ax+by+c 的结果是否都同号即可

  • (3)相交时
    在下图中绘制了所有相交的情况,可以看到 在(2)基础上,直线与矩形框线段有交点时,矩形框2个顶点一定在直线的两侧,或者某个端点在直线上
    在这里插入图片描述
    我们顺时针将矩形框框四个顶点 ( x 1 , y 1 ) , ( x 2 , y 1 ) , ( x 2 , y 2 ) , ( x 1 , y 2 ) (x_1,y_1),(x_2,y_1),(x_2,y_2),(x_1,y_2) (x1,y1),(x2,y1),(x2,y2),(x1,y2) 带入直线方程函数 f ( x , y ) = a x + b y + c f(x,y) = ax+by+c f(x,y)=ax+by+c ,若相邻两个点符号不同(不同为正,不同为负;或 +、-、0的组合),认为是有交点的。

图像中矩形看过顶点已知,且边界线为水平线和垂直线,因此计算简单。
直线方程为 y = k x + b y = kx + b y=kx+b,矩形框边界的线按顺时针排列:水平线 y = y 1 y = y_1 y=y1, 垂直线 x = x 2 x = x_2 x=x2,水平线 y = y 2 y=y_2 y=y2, 垂直线 x = x 1 x = x_1 x=x1。 后面计算过程,不做赘述。

2、图像中的情况

这里直线与矩形框的交点,我们考虑并非线段(若直线延长后与矩形框有交点,如果用两个点描述直线必须要求延长后与矩形框有交点),而是延长后的线段作为实际直线。

通常矩形框在画面中,例如下图中黑框灰色区域(宽度w)为画面区域,蓝色矩形框位于其中。绿色为线段,延长后与矩形框相交为棕色点。
在这里插入图片描述

2.1、常规处理

通常我们计算的逻辑,会考虑直线在画面中的多种情况,例如 垂直、水平,以及实际与画面无交集的无效直线等。
可以参考前述计算方式,这里不在赘叙。

2.2、opencv中的处理

我们直接以图像x方向左、右边界作为线段的延长线,得到两点(0, y0)(w, y1),之后使用cv::clipLine()计算交点。

2.2.1、cv::clipLine函数

先给出函数定义

bool clipLine(  // true if any part of line in imgRect
	cv::Rect imgRect,  // rectangle to clip to
	cv::Point& pt1,  // first endpoint of line overwritten
	cv::Point& pt2  // second endpoint of line, overwritten
);
bool clipLine(  // true if any part of line in image size
	cv::Size imgSize,  // size of image,implies rectangle at 0,0
	cv::Point& pt1,  // first endpoint of line,overwritten
	cv::Point& pt2  // second endpoint of line,overwritten
);

第一种函数的形式使用了cv::Rect,直线和这个矩形比较;第二个函数只有cv::Size,该形式表示矩形的范围是从(0,0)开始的。

只有当喜爱那段完全在指定的矩形范围之外时,函数cv::clipLine()才会返回false

注意,端点是引用,当线段与矩形框相交时,函数返回true,并改写pt1pt2.

2.2.2、测试代码

首先Mat创建一个800x800的画面区域, 之后创建一个 cv::Rect(100, 100, 400, 400) 区域。
我们流程如下:

  1. 直线方程:通过2个点,使用cv::fitLine 拟合出直线的方程
  2. 线段端点:之后以图像左、右边界,确定直线的两个端点
  3. 求解交点: 使用 cv::clipLine 计算矩形区域与线段的交点

完整测试代码如下

   {
        cv::Mat img(800, 800, CV_8UC3, {255,255,255});

        cv::Rect rect(100, 100, 400, 400);

        cv::rectangle(img, rect, {0,0,0}, 3);

        // ---- 直线拟合 
        cv::Point temp[2] = {{100, 100}, {50,-100}};

        cv::Vec4f line_para;
        cv::fitLine(std::vector<cv::Point>{ temp[0], temp[1]/*{-40,-20},{810,830}*/}, line_para, cv::DIST_L2, 0, 1e-2, 1e-2);
        
        // ---- 直线端点
        cv::Point point0(line_para[2], line_para[3]);

        //计算直线的端点(y = k(x - x0) + y0)
        cv::Point point1, point2; 
        if(line_para[1] == 1 || line_para[1] == -1 ) {
            point1 = cv::Point(line_para[2],0);
            point2 = cv::Point(line_para[2], img.rows);
        }
        else {
            float k = line_para[1] / line_para[0];
            point1.x = 0;
            point1.y = k * (0 - point0.x) + point0.y;
            point2.x = img.cols;
            point2.y = k * (point2.x - point0.x) + point0.y;
        }

        cv::Point line1[] = {point1,point2};   
        cv::line(img, line1[0], line1[1], {255,0,0},2);

        // ----- 直线与矩形框交点
        if(cv::clipLine(rect, line1[0], line1[1])) {
            cv::circle(img, line1[0], 3, {0,0,255}, -1);
            cv::circle(img, line1[1], 3, {0,0,255}, -1);
        }


        // -- 绘制
        cv::circle(img, temp[0], 4, {0,255,0}, 2);
        cv::circle(img, temp[1], 4, {0,255,0}, 2);

        cv::imshow("1",img);
        cv::waitKey(0);

        return 0;
    }

2.2.3、测试结果

测试结果如下图。实测任何对与任何直线,以上代码都能适用。

在这里插入图片描述

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1898322.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

Ubuntu / Debian安装FTP服务

本章教程,记录在Ubuntu中安装FTP服务的具体步骤。FTP默认端口:21 1、安装 pure-ftpd sudo apt-get install pure-ftpd2、修改默认配置 # 与 centos 不同,这里需要在 /etc/pure-ftpd/conf 文件夹下执行下列命令,增加对应配置文件: # 创建 /etc/pure-ftpd/conf/PureDB 文件…

4.2 投影

一、投影和投影矩阵 我们以下面两个问题开始&#xff0c;问题一是为了展示投影是很容易视觉化的&#xff0c;问题二是关于 “投影矩阵”&#xff08;projection matrices&#xff09;—— 对称矩阵且 P 2 P P^2P P2P。 b \boldsymbol b b 的投影是 P b P\boldsymbol b Pb。…

金属3D打印如何精准选材

随着3D打印技术的飞跃发展&#xff0c;模具制造领域迎来了前所未有的创新机遇。在众多3D打印技术中&#xff0c;SLM金属3D打印以其精度高、复杂结构成型能力&#xff0c;成为众多行业的优选。然而&#xff0c;金属打印材料&#xff0c;如何精准选择&#xff0c;以最大化满足项目…

ASP.NET MVC-razor编写-2-svg中使用js+添加事件监听

环境&#xff1a;win10 效果 初始状态&#xff1a; 鼠标移入某个text&#xff08;比如KS primer&#xff09;时&#xff0c;text和连接的线条与箭头都变色&#xff1a; 鼠标移出时回复正常。 如果是移入另一种红色的text&#xff08;比如Cell Sceening Tag&#xff09;&…

Python学习笔记29:进阶篇(十八)常见标准库使用之质量控制中的数据清洗

前言 本文是根据python官方教程中标准库模块的介绍&#xff0c;自己查询资料并整理&#xff0c;编写代码示例做出的学习笔记。 根据模块知识&#xff0c;一次讲解单个或者多个模块的内容。 教程链接&#xff1a;https://docs.python.org/zh-cn/3/tutorial/index.html 质量控制…

【数据结构】(6.2)堆的应用——Top-K问题(C语言)

系列文章目录 文章目录 系列文章目录问题引入一、TopK 问题 是什么&#xff1f;二、TopK 问题解决思路2.1 TopK 思路2.2 随机产生数字2.2 完整代码2.3 验证结果 问题引入 TopK 问题 (在一堆数据里面找到前 K 个最大 / 最小的数)。 一、TopK 问题 是什么&#xff1f; 生活中也…

太速科技-FMC209-基于FMC的4路125MAD输入、2路1GDA输出子卡

FMC209-基于FMC的4路125MAD输入、2路1GDA输出子卡 一、板卡概述 本子卡基于FMC连接器实现4路125M采样率AD输出&#xff0c;两路1G采样率DA输出子卡&#xff0c;板卡默认由FMC连接器12V供电&#xff0c;支持外参考时钟&#xff0c;外输入时钟&#xff0c;外触发。 …

全端面试题15(canvas)

在前端开发领域&#xff0c;<canvas> 元素和相关的 API 是面试中经常被提及的主题。下面是一些常见的关于 HTML5 Canvas 的面试问题及解答示例&#xff1a; 1. 什么是 <canvas> 元素&#xff1f; <canvas> 是 HTML5 引入的一个用于图形渲染的标签。它本身并…

使用ChatGPT写论文,只需四步突破论文写作瓶颈!

欢迎关注&#xff0c;为大家带来最酷最有效的智能AI学术科研写作攻略。关于使用ChatGPT等AI学术科研的相关问题可以和作者七哥&#xff08;yida985&#xff09;交流 地表最强大的高级学术AI专业版已经开放&#xff0c;拥有全球领先的GPT学术科研应用&#xff0c;有兴趣的朋友可…

RT-Thread和freeRTOS启动流程

一. freeRTOS启动流程 二. RT-Thread启动流程 因为RT-Thread中我们定义了补丁函数也叫做钩子函数--$Sub$$main()--作为一个新功能函数&#xff0c;可以将原有函数劫持下来&#xff0c;并在之后的程序运行中加上$Super $ $前缀来重新调用原始函数。 所以启动流程是$Sub$$main(…

3033.力扣每日一题7/5 Java

博客主页&#xff1a;音符犹如代码系列专栏&#xff1a;算法练习关注博主&#xff0c;后期持续更新系列文章如果有错误感谢请大家批评指出&#xff0c;及时修改感谢大家点赞&#x1f44d;收藏⭐评论✍ 目录 思路 解题方法 时间复杂度 空间复杂度 Code 思路 首先创建一个与…

CV01_相机成像原理与坐标系之间的转换

目录 0.引言&#xff1a;小孔成像->映射表达式 1. 相机自身的运动如何表征&#xff1f;->外参矩阵E 1.1 旋转 1.2 平移 2. 如何投影到“像平面”&#xff1f;->内参矩阵K 2.1 图像平面坐标转换为像素坐标系 3. 三维到二维的维度是如何丢失的&#xff1f;…

重载一元运算符

自增运算符 #include<iostream> using namespace std; class CGirl { public:string name;int ranking;CGirl() { name "zhongge"; ranking 5; }void show() const{ cout << "name : "<<name << " , ranking : " <…

matplotlib下载安装

matplotlib下载安装过程同之前写的pygame很类似。 Pygame下载安装 python官网 1.搜索matplotlib 直接点进去 查看历史版本&#xff0c;因为新版本可能出现与python不匹配问题。 我选择3.6.3版本&#xff0c;因为我安装的python是3.8&#xff0c;可以匹配版本。同时window操…

C++:拷贝构造函数

拷贝构造函数的引入 用对象来初始化对象 (1)简单变量定义时&#xff0c;可以直接初始化&#xff0c;也可以用另一个同类型变量来初始化。举例说明 (2)用class来定义对象时&#xff0c;可以直接初始化&#xff0c;也可以用另一个对象来初始化。举例说明 testperson xiaohong(na…

搜索+动态规划

刷题刷题刷题刷题 ​​​​​​​​​​​​​​Forgery - 洛谷 | 计算机科学教育新生态 (luogu.com.cn) 思路&#xff1a; 需要两个数组&#xff0c;一个数组全部初始化为".",另一个数组输入数据&#xff0c;每碰到一个“.”就进行染色操作&#xff0c;将其周围的…

鸿蒙:1.入门

概述 简介 鸿蒙操作系统&#xff08;HarmonyOS&#xff09;是华为公司发布的一款智能终端系统&#xff0c;是基于微内核的面向全场景的分布式操作系统。它致力于提供更加安全、高效、低延迟、低功耗的操作体验&#xff0c;可通过技术手段对应用程序和设备进行智能协同&#xf…

「植物大战僵尸杂交版」保姆级攻略大全以及下载指南

植物大战僵尸杂交版自推出以来&#xff0c;以其独特的植物组合和策略玩法&#xff0c;迅速赢得了玩家们的喜爱。如果你正准备加入这场植物与僵尸的战斗&#xff0c;或者已经在战斗中寻求突破&#xff0c;那么这份保姆级的攻略大全将是你的得力助手。同时&#xff0c;我们也提供…

【收藏级神丹】Liae384_刘亦菲_直播可用,平衡度最高的原创神丹,独家珍稀资源

Liae384_刘亦菲_DFL神丹&#xff1a;点击下载 此丹较重&#xff0c;小卡可以使用但不能训练&#xff0c;实测复训适合24G卡8G、12G、16G卡下载练好的专丹直接使用即可384的Liae对各类杂论视频兼容比较好&#xff0c;高参也能容忍高分辨率的DST复用方式: 非必要不用删除AB&…

解决前后端同一个端口跨域问题

前端起了一个代理 如果url是api开头的自动代理访问8080端口&#xff08;解决前后端端口不一致要么是前端代理&#xff0c;要么是后端加过滤器&#xff09; proxy:{/api:{target:http://localhost:8080,changeOrigin : true,// 替换去掉路径上的api// rewrite:(path)>path.r…