单目相机测距

2023-10-29

提示:文章写完后,目录可以自动生成,如何生成可参考右边的帮助文档

一、理论

先说缺点:单目相机测距需要庞大的样本量(需要事先知道物体的大致高度)。
理论依据:小孔成像,相似三角形。
是:像素坐标系==>相机坐标系==>大地坐标系,三个坐标系之间的转换。

1.像素坐标系

通过图像,找到图像的感兴区域所在的真实位置:
请添加图片描述
p为感兴区域,(u,v)为p在像素坐标系下的像素点,(u0,v0)为摄像机光轴和图像的交点,(x,y)为p在图像坐标系下的位置。

1.相机坐标系

请添加图片描述

O2为相机所在位置,(X2,Y2)是相机坐标系,f为焦距(相机到图像的距离)。

1.转换到大地坐标系

请添加图片描述
H是相机离地面的高度,α为相机的俯仰角,(X2,Y2)为大地坐标系,是P在大地的真实位置。

二、代码实现

1.使用matlab找到相机参数

使用matlab里面的图像处理和视觉的工具包:
在这里插入图片描述
导入标定板图片(Add Images):(从网上下载,最好导入15张以上)
在这里插入图片描述
通过Camera Calibrator获取相机的参数(Calibrate):
在这里插入图片描述

通过以上步骤,matlab的工作区里会有我们所需要的相机参数:
在这里插入图片描述
得到内参矩阵(IntrinsicMatrix):
在这里插入图片描述
4378.63641114583 0 0
5.21127842219180 4376.52506346791 0
2862.91763839858 2221.66776871603 1

和畸变系数k1,k2(RadialDistortion)、p1,p2(TangentialDistortion):
k1, k2 = 0.119834359859580, -0.518385077316403
P1, p2 = -8.93400883990887e-05, -0.00140607066202787

内参矩阵和畸变矩阵用来矫正图片。

通过内参矩阵可以得到焦距fx和fy:(测距用)
4378.63641114583 4376.52506346791

2.使用opencv进行测距

1.对标定板进行矫正:

import numpy as np

def show(name):
    cv2.imshow('img', name)
    cv2.waitKey(0)
    cv2.destroyAllWindows()

#内参数矩阵
k = np.array([[4378.63641114583 , 5.21127842219180	, 2862.91763839858	],
                         [0, 4376.52506346791,  2221.66776871603],
                         [0,0,1]])
#畸变系数矩阵,k1 k2 p1 p2 k3
d = np.array([0.119834359859580,  -0.518385077316403,  -8.93400883990887e-05, -0.00140607066202787, 0])
img = cv2.imread(r'./img.jpg')
h,w = img.shape[:2]

mapx, mapy = cv2.initUndistortRectifyMap(k,d,None,k,(w,h),5)

dst = cv2.remap(img, mapx,mapy,cv2.INTER_LINEAR)

h,w = dst.shape[:2]
h_1,w_1 = h/10,w/10
dst = cv2.resize(dst,(int(w_1),int(h_1)))
show(dst)

得到图片:
在这里插入图片描述

2.对目标进行测距(预先知道物体在图像中的感兴区域和物体的大致高度):
代码如下(示例):

import cv2

def show(name):
    cv2.imshow('img', name)
    cv2.waitKey(0)
    cv2.destroyAllWindows()
cap = cv2.imread(r'./cap.jpg')
h,w = cap.shape[:2]
h_1,w_1 = h/10, w/10
cap = cv2.resize(cap,(int(w_1),int(h_1)))
# cap_cutout = cap[95:310,260:400]
cap = cv2.rectangle(cap,(260,95),(400,310),(0,255,0),5)
fx = 4378.63641114583
fy = 4376.52506346791
y = 310-95
x = 400-260
distance = (fy*0.1)/(y*10)  #只求了y坐标的值。
distance = round(distance*100,2)
distance = str(distance)+'cm'
cap = cv2.putText(cap, distance,(260,85),cv2.FONT_HERSHEY_SIMPLEX,1,(0,255,0),2)
print(distance)
show(cap)

得到的结果为:
在这里插入图片描述


总结

单目相机测距需要大量的数据库作为测距的支撑,在项目中也需要通过yolo获取想要测量物体在图像中的感兴区域,通过对水杯的测距得到误差在0.36cm左右。

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

单目相机测距 的相关文章

随机推荐

  • 【OpenCV4】绘制直方图,调用calcHist()函数(图片美化版本)

    文章目录 calcHist 函数详解 绘制直方图的代码 测试图片 生成的直方图 calcHist 函数详解 具体请参考 https blog csdn net qq 42067550 article details 122533459 绘制
  • 论文阅读:Image Partial Blur Detection and Classification

    1 论文主要解决的问题 1 检测图片的模糊区域 2 将模糊区域正确分类 主要分为等向性模糊 near isotropoc blur 和定向运动模糊 directional motion blur 下图 a 为定向运动模糊 b 为等向模糊 2
  • opencv4.x +vs2017 图像学习

    学习来源 哔哩哔哩universityhttps www bilibili com video BV1i54y1m7tw p 2 spm id from pageDriver 学习环境 opencv4 2 vs2017 第一讲总结 imre
  • React中使用react-websocket

    react websocket 安装 npm install save react websocket 用法 import Websocket from react websocket 这里引入的是根目录下index jsx 如果需要引入j
  • springboot链接oracle,配置数据库连接池

    快速搭建一个springboot web项目 SpringBoot中的数据库连接池 SpringBoot链接oracle oracle由于权限问题需要自己下载jar包 然后构建到maven仓库中 下载oracle 11g并安装 安装后拷贝j
  • java在类找不到主方法_错误:在类MyClass中找不到主方法,请将主方法定义为…...

    当您使用java命令从命令行运行Java应用程序 例如 java some AppName arg1 arg2 该命令加载您指定的类 然后查找调用的入口点方法 main 更具体地说 它正在寻找一种声明如下的方法 package some p
  • Fifteenth

    h4 该页面传递一个参数QQ 直线下是接受参数页面的内容 h4 hr
  • linux wget返回值_在Linux命令行上下载文件的5个wget示例

    wget是Linux命令行实用程序 广泛用于从Linux命令行下载文件 有许多选项也可用于从远程服务器下载文件 wget与浏览器窗口中的open url相同 1 使用Wget下载文件 下面的示例将从服务器下载文件到当前本地目录 wget h
  • 【单片机毕业设计】【mcuclub-dz-062】基于单片机的蓝牙语音寻物器

    最近设计了一个项目基于单片机的蓝牙语音寻物器控制系统 与大家分享一下 一 基本介绍 项目名 蓝牙语音寻物器 项目编号 mcuclub dz 062 单片机 STC89C52 功能简介 1 通过蓝牙连接主从机 当主从机配对成功后 主从机的蜂鸣
  • Java求矩形三角形圆形梯形的面积和周长小程序页面

    用java实现图形面积周长计算器 1 一个接口 public interface Shape default void GetArea 面积 default void GetPerimetter 周长 2 圆形类 import java u
  • 孟岩:什么是通证(token),通证可以干什么?

    孟岩是最早将token翻译成为通证的区块链大咖 这个翻译已经得到到了越来越人的认可 原来它叫代币 孟岩建议把它翻译成通证 以下是孟岩关于通证的注解 孟岩 柏链道捷CEO CSDN副总裁 区块链通证派代表人物 通证的本质 通证 token 可
  • 矩阵求逆操作的复杂度分析(逆矩阵的复杂度分析)

    矩阵求逆操作的复杂度分析 逆矩阵的复杂度分析 1 背景 之前写过一篇关于矩阵复杂度分析的文章 没有想到阅读人数那么多 对于IT相关人士来说 从代码层次再结合基本数学知识 就能够很好地理解矩阵的复杂度如何计算得到和分析 其中一位读者提出 矩阵
  • 现代循环神经网络-1.门控循环单元(GRU)【动手学深度学习v2】

    文章目录 1 门控循环单元 GRU 1 1 门控隐状态 A 重置门与更新门 B 候选隐状态 C 隐状态 1 2 GRU的实现 A 从零实现 B 简洁实现 1 门控循环单元 GRU GRU是一个比较新的提出来的 在LSTM之后提出 但是相比L
  • Java接口分页调用的设计

    设计一 Test void t1 int pageNo 0 int pageSize 20 int total do pageNo 调用接口 返回数据 获取总页数 total 102 System out println 第 pageNo
  • 教务管理系统遇到的问题(回答面试官提问用)

    面试官 你做这个项目遇到了什么问题 怎么解决的 我 当时写了一个接口 是关于多个用户查询的 用户有老师 学生和管理员 他们查询的有奖项和文件 奖项和文件是两张表 这造成了很多问题 第一个问题是请求参数过多 有关于文件表和奖项表判断的参数 有
  • 华为OD机试真题-过滤组合字符串【2023Q1】

    题目内容 数字0 1 2 3 4 5 6 7 8 9分别关联 a z 26个英文字母 0 关联 a b c 1 关联 d e f 2 关联 g h i 3 关联 j k l 4 关联 m n o 5 关联 p q r 6 关联 s t 7
  • G2的一些常用配置项

    gt 本文介绍下G2的一些常用配置项 Chart图表的属性 图表的属性直接在代码上介绍 可直接复制到项目上根据实际需求可以选择配置 const chart new G2 Chart container c1 必选 指定图表容器 可以传入该
  • M600Pro操作流程

    ipad平板地面站操作 1 DJI GO 检查项 飞机状态是否为 GPS 模式 且显示起飞准备完毕 2 飞行模式为 P 档 遥控器模式是否正确 遥控器电量是否充足 其他箭头指向处 为正常 电池电量 99 各电池之间差异不超过 3 4 电池电
  • winCE中采用DMA传输数据的方法

    转载请标明是引用于 http blog csdn net chenyujing1234 对于文章中有什么意见或是需要代码的可以留言联系我 1 DMA入口 DDM控制器能服务的外围设备配置信息的结构体 typedef struct BOOL
  • 单目相机测距

    提示 文章写完后 目录可以自动生成 如何生成可参考右边的帮助文档 单目相机测距 一 理论 1 像素坐标系 1 相机坐标系 1 转换到大地坐标系 二 代码实现 1 使用matlab找到相机参数 2 使用opencv进行测距 总结 一 理论 先