CN103438832B - 基于线结构光的三维影像测量方法 - Google Patents
基于线结构光的三维影像测量方法 Download PDFInfo
- Publication number
- CN103438832B CN103438832B CN201310390389.0A CN201310390389A CN103438832B CN 103438832 B CN103438832 B CN 103438832B CN 201310390389 A CN201310390389 A CN 201310390389A CN 103438832 B CN103438832 B CN 103438832B
- Authority
- CN
- China
- Prior art keywords
- line
- dimensional
- structured light
- coordinate
- dimensional image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 23
- 230000003287 optical effect Effects 0.000 claims description 28
- 230000005484 gravity Effects 0.000 claims description 13
- 238000005259 measurement Methods 0.000 claims description 9
- 239000000284 extract Substances 0.000 claims description 7
- 238000000605 extraction Methods 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
Landscapes
- Image Analysis (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Image Processing (AREA)
Abstract
本发明涉及一种,尤其是一种基于线结构光的三维影像测量方法。该方法不仅能利用线结构光的原理快速获取物体表面上的三维数据,同时还能获取空间物体的精确的边缘轮廓数据;同时,该方法步骤简单,在硬件上,只需在现有的三维影像测量仪上增加环形LED光源即可实现。
Description
技术领域
本发明涉及一种,尤其是一种基于线结构光的三维影像测量方法。
背景技术
现有技术中,线结构光测头已经在逆向工程和质量控制领域获得了广泛应用。该线结构光测头与传统的三坐标测量机接触式测量相比大大提高了测量速度。这类测头通常是由一个激光发射器和一个或两个电荷耦合器件(ChargeCoupledDevice,简称CCD)摄像机组成。其中激光发射器发出一条激光线,投射到物体表面;CCD摄像机可以接收被测物体上的漫反射光。通过对该系统进行标定就可将CCD摄像机像面上的像素坐标转化为光平面内的二维坐标。通过运动机构移动物体或移动测头就可以实现对物体的三维扫描测量,所获取的三维数据就是投射到物体上的光条的数据。这类测头的特点是一次可得到一条激光线上的数据,通常是几百个,具有较高的效率,光条上两个点之间的最小距离是由沿光条方向的像素大小决定的。这类测头的主要缺点是不能清晰准确地测量物体的边缘轮廓,物体的边缘通常对应光条断开的地方,这使得该测头无法获取物体边缘轮廓的准确数据。
另外,传统的二维影像测量是利用一个CCD摄像机检测二维平面上物体的边缘轮廓,具有很高的清晰度和精度,在精密检测领域已经得到了广泛的应用。但该方法只能检测同一个平面内的边缘特征,不能检测三维特征。
发明内容
本发明的目的是克服现有技术的上述缺陷,提出了一种将线结构光和影像测量相结合的三维影像测量方法,不仅能利用线结构光的原理快速获取物体表面上的三维数据,同时还能获取空间物体的精确的边缘轮廓数据。
本发明是采用以下的技术方案实现的:一种基于线结构光的三维影像测量方法,包括以下步骤:
(1)关闭环形LED光源,打开线结构光投射器,线结构光投射器发出的平面与被测物体相交产生一条很亮的激光光条,摄像机拍摄该光条的图像,对该光条的图像进行处理从而获得激光光条的中心位置,该中心位置通过O″uv坐标系的像素坐标(u,v)表示;
(2)关闭线结构光投射器,打开环形LED光源,摄像机获取被测物体的图像,对该图像进行亚像素边缘提取,得到物体特征的精确的图像信息;
(3)将步骤(1)的激光光条中心和步骤(2)中物体的亚像素边缘相交得到两个点P和Q,P点和Q点既在激光光条内,又在物体的边缘上,即P点和Q点是在光条上的边缘点,只有在光条上的点才能通过最终转换为空间三维坐标,重复上述三个步骤即可得到获取物体边缘轮廓的二维像面坐标;
(4)通过线结构光“内参数”标定可以根据二维像面坐标得到光平面中的二维坐标,OsXsYs是光平面中的二维坐标系;
(5)通过三维影像测量仪的“外参数”标定可以将光平面中的二维坐标转换到运动机构的三维坐标系中,从而实现三维扫描测量,OmXmYmZm是运动机构的三维坐标系。
本发明中,激光光条中心的提取采用重心法,即投射到物体上的激光光条的亮度近似高斯分布,中间最亮,向两边逐渐变暗,Ni-2到Ni+2是光条上像素位置,Vi-2到Vi+2是对应像素的灰度值,则重心法公式如下:
所述图像亚像素边缘提取采用基于灰度差重心的亚像素算法,设边缘处的灰度值为gi(i=-3,-2,…,3),相应的灰度差值表示为
点P(x,y)为像素级边缘点,点P沿扫描方向选取邻域像素点,对包括P点在内的像素点分别求灰度差fi(i=-2,-1,…,2),图像灰度变化最大的地方由fi关于灰度差坐标i的重心表示,对于左像素级边缘点来说,亚像素边缘提取公式如下:
本发明中内参数的标定过程中,物体边缘点的二维像面坐标确定后,通过f(Pr)=P2d公式得到光平面内的二维坐标,其中Pr是物体边缘点的二维像面坐标,f是该测头内参数所对应的函数,P2d是转换后在光平面中的二维坐标。
所述三维影像测量仪的外参数共有6个分量,包括3个平移分量和3个旋转分量,其中三个平移分量为三维影像测量仪工作原点在世界坐标系中的位置,3个旋转分量为三维影像测量仪二维坐标系相对于世界坐标系的三个旋转角度,设平移分量为M,旋转分量为Q,世界坐标系下的三维坐标系为P3d,则从光平面二维坐标到世界坐标系下的三维坐标为P3d=QP2d+M。
所述三维影像测量仪包括线结构光投射器和摄像机,其中线结构光投射器和摄像机固定在同一个水平基座上,线结构光投射器、摄像机与水平基座一起封装在机壳内,线结构光投射器的投射方向与摄像机的光轴方向的夹角为30°~35°,线结构光投射器的出光孔到摄像机像面中心的距离为90mm~100mm之间,摄像机镜头的前端安装由LED制成的环形LED光源。
本发明的有益效果是:该方法不仅能利用线结构光的原理快速获取物体表面上的三维数据,同时还能获取空间物体的精确的边缘轮廓数据;同时,该方法步骤简单,在硬件上,只需在现有的三维影像测量仪上增加环形LED光源即可实现。
附图说明
图1是基于线结构光的三维影像测量仪的结构示意图;
图2是线结构光内外参数转换关系示意图;
图3(a)是线结构光与物体的相交成像图;
图3(b)是线结构光光条中心的提取结果;
图4是激光光条亮度的分布示意图;
图5(a)是被测物体成像图;
图5(b)是被测物体特征亚像素边缘提取结果;
图6是灰度差重心法的示意图;
图7是结构光光条中心和物体边缘相交的结构示意图。
图1中:1摄像机;2镜头;3线结构光投射器;4环形LED光源。
具体实施方式
下面结合附图对本发明做进一步说明。
本发明所述方法由基于线结构光的三维影像测量仪来完成。所述的基于线结构光的三维影像测量仪包括线结构光投射器3和摄像机1,其中线结构光投射器3和摄像机1固定在同一个水平基座上,并连同该水平基座一起封装在机壳内。线结构光投射器3的投射方向与摄像机1的光轴方向的夹角为30°~35°之间,线结构光投射器3的出光孔到摄像机像面中心的距离为90mm~100mm之间。为了实现影像测量功能,在摄像机镜头的前端安装由红色LED做成的环形LED光源4。
本发明所述的基于线结构光的三维影像测量方法包括以下步骤:
第一步,利用线结构光测量直接由摄像机获取物体边缘轮廓的二维像面坐标,通过图像处理获得激光光条的在像面上的像素坐标(u,v),如图2,O″uv为像面坐标系。其中,获取物体边缘轮廓的二维像面坐标主要包括以下步骤:
(1)关闭环形LED光源4,打开线结构光投射器3,线结构光投射器3发出的平面与被测物体相交产生一条很亮的激光光条,摄像机1拍摄该光条的图像,如图3(a)所示。对该光条的图像进行处理从而获得激光光条的中心位置,如图3(b)所示,该中心位置通过O″uv坐标系的像素坐标(u,v)表示。本发明中,激光光条中心的提取采用重心法,即投射到物体上的激光光条的亮度近似高斯分布,中间最亮,向两边逐渐变暗,如图4所示,Ni-2到Ni+2是光条上像素位置,Vi-2到Vi+2是对应像素的灰度值,则重心法公式如下:
(2)关闭线结构光投射器3,打开环形LED光源4,摄像机1获取被测物体的图像,如图5(a)所示,对该图像进行亚像素边缘提取,得到物体特征的精确的图像信息,如图5(b)所示。本发明中,图像亚像素边缘提取采用基于灰度差重心的亚像素算法,设边缘处的灰度值为gi(i=-3,-2,…,3),相应的灰度差值可表示为
图像的边缘即为图像灰度变化最大的地方,对像素级边缘点进行亚像素提取实际上就是在求得的像素级边缘点的邻域内求出灰度值变化最大的位置,即为对应亚像素边缘点位置。基于上述分析,如图6所示,点P(x,y)为像素级边缘点,点A,B,C,D为点P沿扫描方向的邻域像素点。对上述五点分别求灰度差fi(i=-2,-1,…,2),图像灰度变化最大的地方可由fi关于灰度差坐标i的重心表示。对于左像素级边缘点来说,亚像素边缘提取公式如下:
(3)物体的边缘一般在激光光条断开处,因此将图3(b)中的激光光条中心和图5(b)中物体的亚像素边缘相交可以得到两个点P和Q,如图7所示。P点和Q点既在激光光条内,又在物体的边缘上。经过以上三个步骤可得到物体边缘点的二维像面坐标。
第二步,通过线结构光“内参数”标定可以根据二维像面坐标得到光平面中的二维坐标,OsXsYs是光平面中的二维坐标系。
本发明中才内参数的标定过程中,内参数的标定是建立从CCD像面到激光束或激光平面的映射关系,将CCD所接收的被测物体上的激光线转换为激光平面内的二维坐标。物体边缘点的二维像面坐标确定后,通过f(Pr)=P2d公式得到光平面内的二维坐标,其中Pr是物体边缘点的二维像面坐标,f是该测量仪内参数所对应的函数,P2d是转换后再光平面中的二维坐标。
第三步,通过三维影像测量仪的“外参数”标定可以将光平面中的二维坐标转换到运动机构的三维坐标系中,从而实现三维扫描测量,OmXmYmZm是运动机构的三维坐标系。
三维影像测量仪的外参数共有6个分量,3个平移分量,即三维影像测量仪工作原点在世界坐标系中的位置;3个旋转分量,即测头一维或二维坐标系相对于世界坐标系的三个旋转角度。设平移分量为M,旋转分量为Q,世界坐标系下的三维坐标系为P3d,则从三维影像测量仪二维坐标到世界坐标系下的三维坐标为P3d=QP2d+M。
通过上述步骤即可快速获取物体表面上的三维数据,同时还能获取空间物体的精确的边缘轮廓数据。
Claims (6)
1.一种基于线结构光的三维影像测量方法,其特征在于包括以下步骤:
1)关闭环形LED光源(4),打开线结构光投射器(3),线结构光投射器(3)发出的平面与被测物体相交产生一条激光光条,摄像机(1)拍摄该光条的图像,对该光条的图像进行处理从而获得激光光条的中心位置,该中心位置通过O″uv坐标系的像素坐标(u,v)表示;
2)关闭线结构光投射器(3),打开环形LED光源(4),摄像机(1)获取被测物体的图像,对该图像进行亚像素边缘提取,得到物体特征的精确的图像信息;
3)将步骤1)的激光光条中心和步骤2)中物体的亚像素边缘相交得到两个点P和Q,P点和Q点既在激光光条内,又在物体的边缘上,也就是说P点和Q点是在光条上的边缘点,重复上述两个步骤即可得到获取物体边缘轮廓的二维像面坐标;
4)通过线结构光内参数标定可以根据二维像面坐标得到光平面中的二维坐标,OsXsYs是光平面中的二维坐标系;
5)通过三维影像测量仪的外参数标定可以将光平面中的二维坐标转换到运动机构的三维坐标系中,从而实现三维扫描测量,OmXmYmZm是运动机构的三维坐标系。
2.根据权利要求1所述的基于线结构光的三维影像测量方法,其特征在于:所述激光光条中心的提取采用重心法,即投射到物体上的激光光条的亮度近似高斯分布,中间最亮,向两边逐渐变暗,Ni-2到Ni+2是光条上像素位置,Vi-2到Vi+2是对应像素的灰度值,则重心法公式如下:
3.根据权利要求1所述的基于线结构光的三维影像测量方法,其特征在于:所述图像亚像素边缘提取采用基于灰度差重心的亚像素算法,设边缘处的灰度值为gi(i=-3,-2,…,3),相应的灰度差值表示为:
点P(x,y)为像素级边缘点,点P沿扫描方向选取邻域像素点,对包括P点在内的像素点分别求灰度差fi(i=-2,-1,…,2),图像灰度变化最大的地方由fi关于灰度差坐标i的重心表示,对于左像素级边缘点来说,亚像素边缘提取公式如下:
4.根据权利要求1所述的基于线结构光的三维影像测量方法,其特征在于:所述内参数的标定过程中,物体边缘点的二维像面坐标确定后,通过f(Pr)=P2d公式得到光平面内的二维坐标,其中Pr是物体边缘点的二维像面坐标,f是该测量仪内参数所对应的函数,P2d是转换后在光平面中的二维坐标。
5.根据权利要求1所述的基于线结构光的三维影像测量方法,其特征在于:所述三维影像测量仪的外参数共有6个分量,包括3个平移分量和3个旋转分量,其中三个平移分量为三维影像测量仪工作原点在世界坐标系中的位置,3个旋转分量为三维影像测量仪一维或二维坐标系相对于世界坐标系的三个旋转角度,设平移分量为M,旋转分量为Q,世界坐标系下的三维坐标系为P3d,则从光平面二维坐标到世界坐标系下的三维坐标为P3d=QP2d+M。
6.根据权利要求1或5所述的基于线结构光的三维影像测量方法,其特征在于:所述三维影像测量仪包括线结构光投射器(3)和摄像机(1),线结构光投射器(3)和摄像机(1)固定在同一个水平基座上,线结构光投射器(3)、摄像机(1)与水平基座一起封装在机壳内,线结构光投射器(3)的投射方向与摄像机(1)的光轴方向的夹角为30°~35°,线结构光投射器(3)的出光孔到摄像机像面中心的距离为90mm~100mm之间,摄像机镜头的前端安装由LED制成的环形LED光源(4)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310390389.0A CN103438832B (zh) | 2013-08-30 | 2013-08-30 | 基于线结构光的三维影像测量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310390389.0A CN103438832B (zh) | 2013-08-30 | 2013-08-30 | 基于线结构光的三维影像测量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103438832A CN103438832A (zh) | 2013-12-11 |
| CN103438832B true CN103438832B (zh) | 2016-01-20 |
Family
ID=49692528
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310390389.0A Expired - Fee Related CN103438832B (zh) | 2013-08-30 | 2013-08-30 | 基于线结构光的三维影像测量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103438832B (zh) |
Families Citing this family (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6320051B2 (ja) * | 2014-01-17 | 2018-05-09 | キヤノン株式会社 | 三次元形状計測装置、三次元形状計測方法 |
| CN104677305B (zh) * | 2015-02-11 | 2017-09-05 | 浙江理工大学 | 一种基于十字结构光的物体表面三维重建方法和系统 |
| CN106705879B (zh) * | 2015-08-06 | 2019-09-17 | 赛特株式会社 | 利用三维形状测量仪的物体测量方法及用于该测量方法的三维形状测量仪 |
| CN105290621B (zh) * | 2015-10-12 | 2017-07-11 | 深圳市海目星激光科技有限公司 | 一种基于视觉引导的高速高精度极耳切割方法和设备 |
| CN106109015A (zh) * | 2016-08-18 | 2016-11-16 | 秦春晖 | 一种头戴式医疗系统及其操作方法 |
| CN106959078B (zh) * | 2017-02-28 | 2019-07-30 | 苏州凡目视觉科技有限公司 | 一种用于三维轮廓测量的轮廓测量方法 |
| CN110709669B (zh) * | 2017-06-07 | 2021-09-07 | 索尼半导体解决方案公司 | 信息处理装置和方法 |
| CN107941170A (zh) * | 2017-12-15 | 2018-04-20 | 武汉夕睿光电技术有限公司 | 道路检测装置及系统 |
| CN108428251A (zh) * | 2018-03-09 | 2018-08-21 | 深圳市中捷视科科技有限公司 | 一种基于机器视觉技术激光结构光自动标定方法 |
| CN108332684A (zh) * | 2018-04-23 | 2018-07-27 | 中国科学院苏州生物医学工程技术研究所 | 一种基于结构光照明显微技术的三维轮廓测量方法 |
| CN109269952A (zh) * | 2018-10-16 | 2019-01-25 | 东南大学 | 一种丝状物料在气流床矩形提升管中三维颗粒浓度信息的测量装置和方法 |
| CN110375674A (zh) * | 2019-07-02 | 2019-10-25 | 东莞理工学院 | 一种精密制造设备的视觉检测系统 |
| CN110702002B (zh) * | 2019-07-04 | 2020-12-11 | 天津大学 | 一种曲面透光件多层同步测量方法 |
| CN110988145A (zh) * | 2019-10-21 | 2020-04-10 | 深圳先进技术研究院 | 扫描超声检测系统 |
| CN113465536A (zh) * | 2021-06-30 | 2021-10-01 | 皖江工学院 | 基于摄像头引导的激光云台及其工作方法 |
| CN113390357B (zh) * | 2021-07-08 | 2022-06-07 | 南京航空航天大学 | 一种基于双目多线结构光的铆钉齐平度测量方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102589476A (zh) * | 2012-02-13 | 2012-07-18 | 天津大学 | 高速扫描整体成像三维测量方法 |
| CN102661956A (zh) * | 2012-04-24 | 2012-09-12 | 浙江大学 | 超光滑表面缺陷检测系统及其畸变校正方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5161024B2 (ja) * | 2008-10-10 | 2013-03-13 | 富士フイルム株式会社 | 3次元形状計測用撮影装置および方法並びにプログラム |
| JP4797109B2 (ja) * | 2010-03-26 | 2011-10-19 | Necエンジニアリング株式会社 | 三次元形状計測装置及び三次元形状計測方法 |
-
2013
- 2013-08-30 CN CN201310390389.0A patent/CN103438832B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102589476A (zh) * | 2012-02-13 | 2012-07-18 | 天津大学 | 高速扫描整体成像三维测量方法 |
| CN102661956A (zh) * | 2012-04-24 | 2012-09-12 | 浙江大学 | 超光滑表面缺陷检测系统及其畸变校正方法 |
Non-Patent Citations (1)
| Title |
|---|
| 《基于线结构光的复合式视觉测量技术的研究》;张梅风;《中国优秀硕士学位论文全文数据库 信息科技辑》;20090215;I138-487 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103438832A (zh) | 2013-12-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103438832B (zh) | 基于线结构光的三维影像测量方法 | |
| CN102435770B (zh) | 一种单相机三维体视粒子图像测速系统 | |
| CN106767527B (zh) | 一种三维轮廓的光学混合检测方法 | |
| CN105203034B (zh) | 一种基于单目摄像头三维测距模型的测高测面积方法 | |
| US9454821B2 (en) | One method of depth perception based on binary laser speckle images | |
| CN103971406B (zh) | 基于线结构光的水下目标三维重建方法 | |
| CN103759671B (zh) | 一种牙模三维表面数据非接触扫描方法 | |
| Massot-Campos et al. | Underwater laser-based structured light system for one-shot 3D reconstruction | |
| CN107578464A (zh) | 一种基于线激光扫描的传送带工件三维轮廓测量方法 | |
| CN107167118B (zh) | 一种基于非编码平行多线的稳定实时激光测量方法 | |
| CN103499297A (zh) | 一种基于ccd的高精度测量方法 | |
| CN111260720A (zh) | 一种基于深度学习方法的目标高度测定系统 | |
| CN104567728A (zh) | 激光视觉轮廓测量系统及测量方法、立体靶标 | |
| CN106289099B (zh) | 一种单目视觉系统及基于该系统的三维尺寸快速测量方法 | |
| CN111091076B (zh) | 基于立体视觉的隧道限界数据测量方法 | |
| CN105069789B (zh) | 基于编码网格模板的结构光动态场景深度获取方法 | |
| CN108362228B (zh) | 一种基于双光机的光刀光栅混合式三维测量装置及测量方法 | |
| KR20160121509A (ko) | 2개의 카메라로부터의 곡선의 세트의 구조형 광 매칭 | |
| CN106595519A (zh) | 一种基于激光mems投影的柔性三维轮廓测量方法及装置 | |
| CN106996748A (zh) | 一种基于双目视觉的轮径测量方法 | |
| CN101825441A (zh) | 一种光电铁路轨距-站台间距测量方法 | |
| CN103852060A (zh) | 一种基于单目视觉的可见光图像测距方法 | |
| TWI659390B (zh) | 應用於物件檢測之攝影機與雷射測距儀的數據融合方法 | |
| CN104236479A (zh) | 一种线结构光三维测量系统及3d纹理图像构造算法 | |
| Zhang et al. | Accuracy improvement in laser stripe extraction for large-scale triangulation scanning measurement system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160120 Termination date: 20170830 |