工业现场相机坐标系和机械手坐标系的标定_labview 9点标定-程序员宅基地

技术标签: opencv  手眼标定  

原文:https://blog.csdn.net/kaychangeek/article/details/73878994

参考:https://blog.csdn.net/qq_16481211/article/details/79764730

工业现场使用视觉时一般需要相机坐标系和机械手臂坐标系的转化,这里介绍一种比较简单的标定方案。没有使用到标定板。经过几个项目的测试,精度还算可以,如果要求高精度的场合,就用标定板标定吧!【可以购买专用的标定板,或者自己制作(像我这种穷逼),哈哈,可以阅读这个文章这个软件自己做一个标定板嘿!】

 

        如上图所示:OXY为机械手坐标系,O'X'Y'为相机坐标系。theta为两个坐标系之间的夹角。假设P点在图像上的位置如图,则P在机械手坐标系有一个坐标,在图像坐标系也有一个坐标。我们要做的工作就是图像上的任意一点都可以转化为机械手坐标系上的坐标点:P(Machine) = f(P(Image))。接下来介绍如何找到这个关系。

OXY为机械手坐标系 O'X'Y'为相机坐标系从上图可以看出坐标转化关系:

x = x' * r * cos(theta) - y' * r * sin(theta) + x0;

y = x' * r * sin(theta) + y' * r * cos(theta) + y0;

其中r是毫米像素比、(mm/pixel)就是一个毫米有几个像素,theta为两个坐标系之间的夹角,(x0,y0)为图像坐标原点到机械坐标原点的距离。

简化抽象公式,假设:

a = r * cos(theta);

b = r * sin(theta);

c = x0;

d = y0;

得到:

 

x = x' * a - y' * b + c;

y = x' * b + y' * a + d;

很显然,要解出这个方程,需要两组对应关系,就是两组对应的坐标点。设两组坐标点,如下:

第一组:图像坐标点:(xImage1,yImage1) 对应的机械坐标点:(xMachine1,yMachine1)

第二组:图像坐标点:(xImage2,yImage2) 对应的机械坐标点:(xMachine2,yMachine2)

则可以解出a ,b, c, d。如下:

 

 
  1. a = ((xMachine1 - xMachine2)*(xImage1- xImage2) + (yMachine1 - yMachine2)*(yImage1 - yImage2))   

  2.     / ((xImage1 - xImage2)*(xImage1 - xImage2) + (yImage1 - yImage2)*(yImage1 - yImage2));  

  3. b = ((yMachine1 - yMachine2)*(xImage1 - xImage2) - (xMachine1 - xMachine2)*(yImage1 - yImage2))   

  4.     / ((xImage1 - xImage2)*(xImage1 - xImage2) + (yImage1 - yImage2)*(yImage1 - yImage2));  

  5. c = xMachine1 - a*xImage1 + b*yImage1;  

  6. d = yMachine1 - b*xImage1 - a*yImage1;

 

所以,就得出了图像上任意一点的像素坐标转成机械手坐标的关系。

 

以下是我写的一个求解a ,b, c, d软件:需要的点这里下载。下面是软件截图:

 

 

        下面举个栗子说明一下操作吧!这是我实际项目中的机械手和相机布局情况,画图真T``M`累!,如下图:

 

 

        首先将产品上的目标点搞到机械手的Z轴的中心,然后移动到相机事业范围内,让识别产品上的目标点,这时,你会读到一组机械手坐标(也就是上面所说的(xMachine1,yMachine1)),和一组相机坐标(也就是上面所说的(xImage1,yImage1)),再在相机视野范围内移动机械手,就可以得到第二组数据(xMachine2,yMachine2)和(xImage2,yImage2),写入标定软件就可以求出a, b, c, d。就求出了图像上任意一点对应的机械坐标,然后你下次移动第二个产品到相机视野时,首先识别目标点的像素坐标,经过对应的转换就可以得出机械坐标,然后进行相应的偏移即可实现你想要的操作!

        当然,有些时候,实际现场并没有如你所愿,现在出现一种情况就是以机械手自带的Z轴无法移动到视野中去,这样就无法进行上面的操作了,系不系!这也是我遇到的情况,SO,我的解决方案如下:先上示意图:

如图,就是在机械手Z轴的横杆上加上一个横条,使得产品可以移动到视野范围内。然后进行上面的操作,虽然现在可以将产品移动到视野范围内了,但是你会发现此时机械手的坐标并不是产品的坐标,因为他们之间隔着一个横条,那我们要怎么换算过来呢!

        首先,我们可以这么假定的认为,他们是机械手的坐标就是产品的坐标,然后就可以用上面取两组对应点的方法求出图像上每个点的机械坐标了,但是,此时的机械坐标并不是真正的机械坐标,但是,我们可以利用这些坐标找出机械手U轴的旋转中心就是Z轴(就是实际机械坐标所在的位置)在我们所建立的不是真正的机械坐标系中的坐标。是不是有点不好理解,我也不大清楚怎么描述!求U轴的旋转中心的过程如下:上图:

        如图:将产品的目标点绕机械手的U轴的旋转中心旋转得到3的目标点的不是真正的机械坐标系的坐标P1,P2,P3,必须让每次旋转的目标点在视野范围内,然后通过圆弧上的3点就可以求得圆心。此时的圆心坐标为不是真正的机械坐标系的坐标。然后我们就可以求出产品中心和机械手U轴旋转中心的deltaX,deltaY,然后就可以结合真正的机械坐标系建立工件坐标系,要对机械手有些了解才比较好理解这些拗口的话!然后对于每一个新的产品都可以建立合适的工件坐标系,然后进行补正和一些操作!(上面那个软件也有三点求圆心的工具)

        或者,最近我制作了一个更加强大一点的软件(并且修复上面的哪个软件不能输入负数的缺陷),也可以直接算出这边第二种情况的对应关系(详细细节见该软件的帮助信息),看下这篇文章的介绍,里面有提供下载地址。略贵!哈哈,不过毕竟是辛辛苦苦做出来哒!

        相机绑定在机械手上的标定方法参考这里工业现场相机坐标系和机械手坐标系的标定(2)-相机和机械手绑定的情况,以上是我最近总结的心得体会,希望对大家有所帮助!如果帮助到您,点个赞!^-^

-------------------------------------------------------------------------------------------------

 

工业现场相机坐标系和机械手坐标系的标定(2)-相机和机械手绑定的情况

 之前,我写过一篇标定文章,但是只说明了相机和机械手分离的情况的标定方法,为了提供完整的标定方案,这边介绍一下怎么利用我做的标定助手完成相机和机械手绑定的情况的标定方法。

        由文章“工业现场相机坐标系和机械手坐标系的标定(1)-相机和机械手分离的情况”我们知道只要找到两组对应点即可完成相机坐标系和机械手坐标系的关系转换。那么,现在相机跟随着机械手运动(它们绑定在一起了),如果还将标定靶标也固定在机械手上的话,三个对象都相对运动,就无法标定了!那要怎么利用标定助手完成标定呢?经过思考,我给出了以下标定思路:

整个标定系统的初始状态

        相机和机械手绑定在一起并且机械手旋转中心处于机械原点,这时候相机中心和机械手中心有一个相对的偏移并且是固定的(相机中心和机械手旋转中心的相对距离设为(detX,detY)),将标定靶标放置在工作区域(实际机械手的工作区域)中的合适位置。

第一步:求取第一组图像点坐标C1

        移动机械手(相机会跟着移动)到相机能够清晰拍摄到标定靶标的位置(设该位置的机械坐标为(x0,y0))。然后调整标定靶标的识别参数,使得能够实时的识别标定靶标的中心位置,手动移动靶标最好让它远离图像中心,越远越好(精度越高)但是不能超出相机的视野范围,我就把靶标放在相机视野的左上角的区域吧,然后我们识别当前靶标的图像位置C1,就是第一组数据的图像点。

第二步:求取第一组机械点坐标M1

        进入“标定”选项卡,勾选“显示十字标”,这时候会在实时图像中间显示一个蓝色的十字光标,然后,慢速移动机械手,直到C1点和十字光标的交点重合(也就是说移动机械手让实时识别到的靶标中心坐标处于图像的中心,比如采集的图像为2048*1536大小的,你需要缓慢移动机械手(靶标保持不动)直到靶标的中心处于(1024,768)的这个位置),然后记下当前的机械手坐标也就是第一组数据的机械点M1(注意:此时实际是让相机中心和靶标中心重合,并非机械手旋转中心和靶标中心重合,它们之间差了一个detX,detY的距离,具体怎么换算,等下再说)。

第三步:求取第二组图像点坐标C2

        将机械手移动回到初始的工作坐标点(x0,y0),将靶标放在相机视野的右下角的区域,识别当前的靶标图像坐标C2,就是第二组数据的图像点。

第四步:求取第二组机械点坐标M2

        同理,慢速移动机械手,直到C2点和十字光标的交点重合(就是移动机械手让实时识别到的靶标中心坐标处于图像的中心),然后记下当前的机械手坐标也就是第二组数据的机械点M2。

第五步:求取整个机械手工作区域坐标系和相机坐标系的对应关系

        不知道你有没有发现,上面求得工作区域只对工作区域(x0,y0)开始到相机视野结束的位置有效(比如你在机械手移动到工作区域(x0,y0)为起点的位置识别靶标的图像位置坐标为(100,100),但是你将机械手移动到工作区域(x1,y1)为起点的位置识别靶标的图像位置坐标也可能为(100,100))。这样换算过来的机械坐标就是错误的,其实我们上面求得是小区域的转换关系,那我们要怎么扩展到整个机械手的工作区域呢?我们只需要这样做:在机械手工作的时候是可以知道自己在哪里的,比如:现在机械手移动到工作区域(x1,y1)为起点的位置识别靶标的图像位置坐标为(100,100),我们将(100,100)代入转换关系trans得出来的机械坐标是对应(x0,y0)的,要转换到对应(x1,y1)的只需要将换算处理的坐标加上(x1 - x0,y1 - y0)即可!抽象的公式如下:

 

 
  1. current_machine.x = trans(current_image.x,current_image.y) + (x1 - x0);

  2. current_machine.y = trans(current_image.x,current_image.y) + (y1 - y0);

其中trans函数就是将图像坐标转换为相对于(x0,y0)工作区域的机械坐标,再加上坐标(x1 - x0,y1 - y0)就是当前正确的机械坐标。到此,我们得到了图像坐标和整个机械手工作区域的坐标的关系,但是此时的机械坐标指的是相机的中心,并不是机械手的旋转中心。

 

第六步:将相机中心转换为机械手的旋转中心

        首先,你可以利用目测法,游标卡尺法等等粗略的算出detX,detY的值,然后,进行实时识别和移动测试进行细调,因为粗略的估算有误差,你可以经过多次微调直到精确到机械手的旋转中心移动到产品的上方就成功了。记下此时的detX,detY,代入以下公式就可以算出了最终的机械手应该去的地方了。最终的换算公式如下:

 

 

 
  1. current_machine.x = trans(current_image.x,current_image.y) + (x1 - x0) + detX;

  2. current_machine.y = trans(current_image.x,current_image.y) + (y1 - y0) + detY;

        整个过程的伪代码运算过程如下:

 

 

 
  1. [a,b,c,d] = f([C1,M1],[C2,M2])

  2. trans(double xImage,double yImage){

  3. double xMachine = xImage * a - yImage * b + c;

  4. double yMachine = xImage * b + yImage * a + d;

  5. return xMachine,yMachine;

  6. }

  7. current_machine.x = trans(current_image.x,current_image.y) + (x1 - x0) + detX;

  8. current_machine.y = trans(current_image.x,current_image.y) + (y1 - y0) + detY;

        利用标定助手求取[a,b,c,d]外参矩阵的输入参数[C1,M1]和[C2,M2]的操作方法如下图所示:

 

 

        可能说的比较抽象,请见谅,后期有时间画几张过度图解释!

-------------------------------------------------------------------------------------------------------------------------------------------------------------------------

六轴机器人的标定

所谓手眼系统,就是人眼镜看到一个东西的时候要让手去抓取,就需要大脑知道眼镜和手的坐标关系。如果把大脑比作B,把眼睛比作A,把手比作C,如果A和B的关系知道,B和C的关系知道,那么C和A的关系就知道了,也就是手和眼的坐标关系也就知道了。 
这里写图片描述
相机知道的是像素坐标,机械手是空间坐标系,所以手眼标定就是得到像素坐标系和空间机械手坐标系的坐标转化关系。 
在实际控制中,相机检测到目标在图像中的像素位置后,通过标定好的坐标转换矩阵将相机的像素坐标变换到机械手的空间坐标系中,然后根据机械手坐标系计算出各个电机该如何运动,从而控制机械手到达指定位置。这个过程中涉及到了图像标定,图像处理,运动学正逆解,手眼标定等。

常用的标定方法有:九点标定

九点标定:

这里写图片描述
九点标定直接建立相机和机械手之间的坐标变换关系。 
让机械手的末端去走这就9个点得到在机器人坐标系中的坐标,同时还要用相机识别9个点得到像素坐标。这样就得到了9组对应的坐标。 
由下面的式子可知至少需要3个点才能求出标定的矩阵。 
这里写图片描述

(1)、标定,Halcon中进行9点标定的算子

(2)、求解

affine_trans_point_2d(HomMat2D,Row2,Column2,Qx,Qy)
%由像素坐标和标定矩阵求出机器人基础坐标系中的坐标

一些特殊情况的解释:

有些情况中我们看到相机固定在一个地方,然后拍照找到目标,控制机械手去抓取,这种就很好理解。我们也叫做eye-to-hand 
还有一种情况是相机固定在机械手上面,这种情况的标定过程实际上和相机和机械手分离的标定方法是一样的,因为相机拍照时,机械手会运动到相机标定的时候的位置,然后相机拍照,得到目标的坐标,再控制机械手,所以简单的相机固定在末端的手眼系统很多都是采用这种方法,标定的过程和手眼分离系统的标定是可以相同对待的。我们也叫做eye-in-hand
这里写图片描述 
这里写图片描述

 

 

 

 

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/liangchunjiang/article/details/81031027

智能推荐

c# 调用c++ lib静态库_c#调用lib-程序员宅基地

文章浏览阅读2w次,点赞7次,收藏51次。四个步骤1.创建C++ Win32项目动态库dll 2.在Win32项目动态库中添加 外部依赖项 lib头文件和lib库3.导出C接口4.c#调用c++动态库开始你的表演...①创建一个空白的解决方案,在解决方案中添加 Visual C++ , Win32 项目空白解决方案的创建:添加Visual C++ , Win32 项目这......_c#调用lib

deepin/ubuntu安装苹方字体-程序员宅基地

文章浏览阅读4.6k次。苹方字体是苹果系统上的黑体,挺好看的。注重颜值的网站都会使用,例如知乎:font-family: -apple-system, BlinkMacSystemFont, Helvetica Neue, PingFang SC, Microsoft YaHei, Source Han Sans SC, Noto Sans CJK SC, W..._ubuntu pingfang

html表单常见操作汇总_html表单的处理程序有那些-程序员宅基地

文章浏览阅读159次。表单表单概述表单标签表单域按钮控件demo表单标签表单标签基本语法结构<form action="处理数据程序的url地址“ method=”get|post“ name="表单名称”></form><!--action,当提交表单时,向何处发送表单中的数据,地址可以是相对地址也可以是绝对地址--><!--method将表单中的数据传送给服务器处理,get方式直接显示在url地址中,数据可以被缓存,且长度有限制;而post方式数据隐藏传输,_html表单的处理程序有那些

PHP设置谷歌验证器(Google Authenticator)实现操作二步验证_php otp 验证器-程序员宅基地

文章浏览阅读1.2k次。使用说明:开启Google的登陆二步验证(即Google Authenticator服务)后用户登陆时需要输入额外由手机客户端生成的一次性密码。实现Google Authenticator功能需要服务器端和客户端的支持。服务器端负责密钥的生成、验证一次性密码是否正确。客户端记录密钥后生成一次性密码。下载谷歌验证类库文件放到项目合适位置(我这边放在项目Vender下面)https://github.com/PHPGangsta/GoogleAuthenticatorPHP代码示例://引入谷_php otp 验证器

【Python】matplotlib.plot画图横坐标混乱及间隔处理_matplotlib更改横轴间距-程序员宅基地

文章浏览阅读4.3k次,点赞5次,收藏11次。matplotlib.plot画图横坐标混乱及间隔处理_matplotlib更改横轴间距

docker — 容器存储_docker 保存容器-程序员宅基地

文章浏览阅读2.2k次。①Storage driver 处理各镜像层及容器层的处理细节,实现了多层数据的堆叠,为用户 提供了多层数据合并后的统一视图②所有 Storage driver 都使用可堆叠图像层和写时复制(CoW)策略③docker info 命令可查看当系统上的 storage driver主要用于测试目的,不建议用于生成环境。_docker 保存容器

随便推点

网络拓扑结构_网络拓扑csdn-程序员宅基地

文章浏览阅读834次,点赞27次,收藏13次。网络拓扑结构是指计算机网络中各组件(如计算机、服务器、打印机、路由器、交换机等设备)及其连接线路在物理布局或逻辑构型上的排列形式。这种布局不仅描述了设备间的实际物理连接方式,也决定了数据在网络中流动的路径和方式。不同的网络拓扑结构影响着网络的性能、可靠性、可扩展性及管理维护的难易程度。_网络拓扑csdn

JS重写Date函数,兼容IOS系统_date.prototype 将所有 ios-程序员宅基地

文章浏览阅读1.8k次,点赞5次,收藏8次。IOS系统Date的坑要创建一个指定时间的new Date对象时,通常的做法是:new Date("2020-09-21 11:11:00")这行代码在 PC 端和安卓端都是正常的,而在 iOS 端则会提示 Invalid Date 无效日期。在IOS年月日中间的横岗许换成斜杠,也就是new Date("2020/09/21 11:11:00")通常为了兼容IOS的这个坑,需要做一些额外的特殊处理,笔者在开发的时候经常会忘了兼容IOS系统。所以就想试着重写Date函数,一劳永逸,避免每次ne_date.prototype 将所有 ios

如何将EXCEL表导入plsql数据库中-程序员宅基地

文章浏览阅读5.3k次。方法一:用PLSQL Developer工具。 1 在PLSQL Developer的sql window里输入select * from test for update; 2 按F8执行 3 打开锁, 再按一下加号. 鼠标点到第一列的列头,使全列成选中状态,然后粘贴,最后commit提交即可。(前提..._excel导入pl/sql

Git常用命令速查手册-程序员宅基地

文章浏览阅读83次。Git常用命令速查手册1、初始化仓库git init2、将文件添加到仓库git add 文件名 # 将工作区的某个文件添加到暂存区 git add -u # 添加所有被tracked文件中被修改或删除的文件信息到暂存区,不处理untracked的文件git add -A # 添加所有被tracked文件中被修改或删除的文件信息到暂存区,包括untracked的文件...

分享119个ASP.NET源码总有一个是你想要的_千博二手车源码v2023 build 1120-程序员宅基地

文章浏览阅读202次。分享119个ASP.NET源码总有一个是你想要的_千博二手车源码v2023 build 1120

【C++缺省函数】 空类默认产生的6个类成员函数_空类默认产生哪些类成员函数-程序员宅基地

文章浏览阅读1.8k次。版权声明:转载请注明出处 http://blog.csdn.net/irean_lau。目录(?)[+]1、缺省构造函数。2、缺省拷贝构造函数。3、 缺省析构函数。4、缺省赋值运算符。5、缺省取址运算符。6、 缺省取址运算符 const。[cpp] view plain copy_空类默认产生哪些类成员函数

推荐文章

热门文章

相关标签