CN103234475B - 一种基于激光三角测量法的亚像素表面形态检测方法 - Google Patents
一种基于激光三角测量法的亚像素表面形态检测方法 Download PDFInfo
- Publication number
- CN103234475B CN103234475B CN201210490109.9A CN201210490109A CN103234475B CN 103234475 B CN103234475 B CN 103234475B CN 201210490109 A CN201210490109 A CN 201210490109A CN 103234475 B CN103234475 B CN 103234475B
- Authority
- CN
- China
- Prior art keywords
- laser
- sub
- image
- method based
- pix
- 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.)
- Active
Links
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提供了一种基于激光三角测量法的亚像素表面形态检测方法,该方法通过测量零件不同位置轮廓切片,叠加得到整个零件轮廓。本发明充分利用图像的彩色信息,在激光条纹是红色,背景是黑色的先验知识指导下,提取RGB图像的红色通道分量,减小噪声影响,提高高斯拟合条纹中心精度。通过对同一位置多次拍摄叠加激光条纹图像提高精度,降低了对激光条纹宽度和亮度的要求,降低了激光器的成本。通过先粗略确定条纹中心位置然后采用固定数量的像素拟合确定亚像素中心的方法避免了确定激光条纹边界而采用复杂算法,降低了算法复杂度。
Description
技术领域
本发明涉及非接触测量领域,特别涉及一种基于激光三角测量法的亚像素精度表面形态测量和分析检测方法。
背景技术
近年来非接触式测量越来越多的受到重视,与传统接触式表面形态测量和分析方法相比,激光非接触式测量具有独特优势:1、传统测量需要接触工件表面,不但会造成测量设备的磨损和损害,还容易划伤工件表面,而激光非接触方法避免了对被测物体造成划痕和磨损,尤其适用于各种柔软材料、易腐蚀材料和传统方式无法检测的表面形态测量和分析;2、接触式表面形态测量和分析设备需要工作过程中多次校准,而对于特定的测量环境,激光非接触测量设备只需一次性设置和校准,终身受用,节约维护和再标定的时间和金钱成本;3、非接触测量设备可以安装在生产线上,实现在线测量。
目前激光非接触测量主要使用的激光器分为点激光和一字线激光。一字线激光由于采样面积大,在测量精度和使用范围上较点激光有明显优势。但目前的使用一字线激光器的激光三角测量法的很大程度上精度取决于激光器条纹的粗细、均匀性和亮度分布情况。使用性能更优的激光器来提高测量精度的方法增加了设备的成本。
在激光三角测量法中,激光条纹的中心代表工件表面轮廓。高斯拟合法为确定激光条纹中心的主要方法,但目前的方法有以下几个不足:1、不同表面激光散射情况不同,仅使用单幅激光条纹图像拟合的情况下,为了提高照射亮度,很容易产生过曝光情况,导致估计条纹中心产生错误;2、通过多次测量来提高测量精度,并没有消除每次的测量误差;3、处理激光条纹图像时直接把彩色图转成灰度图处理,并没有利用图像的彩色信息;4、通过确定激光条纹的边缘范围来确定高斯拟合所需像素点数量,这种方法受到了边缘范围精确度限制,并增加了计算时间。
发明内容
本发明的目的在于针对上述存在问题和不足,提出一种基于激光三角测量法的亚像素表面形态检测方法,该方法具有很好的实用性,使用普通的一字线激光器就能实现零件表面形态如高度、粗糙度、平面度、轮廓的准确检测,大大降低了激光器的使用成本。
本发明为了完成其技术目的所采用的技术方案是:一种基于激光三角测量法的亚像素表面形态检测方法,使用一字线激光器产生的激光进行三角测量,包括以下步骤:
A、采用一字线激光器产生的激光照射零件;
B、采用彩色CCD相机对零件的同一位置采集2-10幅激光条纹图片;并将采集到的所有激光条纹图片叠加形成一副叠加图像;
C、提取叠加图像的红色通道图像并对其进行中值滤波,滤除噪声;
D、沿进行中值滤波后的红色通道图像列方向搜索每一列像素极大值点;
E、在每一列像素极大值点两侧各取5-10个像素点进行高斯曲线拟合,提取亚像素精度的激光条纹中心线,形成亚像素表面形态图像。
进一步的,上述的基于激光三角测量法的亚像素表面形态检测方法中:在步骤A之前,还需要对测量测量设备进行校准和对测量***进行标定;
校准设备包括使一字线激光器和水平工作台表面垂直,CCD摄像机主光轴与激光器在一个平面内,与水平表面垂线成45~60度角;
对***进行标定包括得到真实空间高度与CCD拍摄图像高度之间的关系。
进一步的,上述的基于激光三角测量法的亚像素表面形态检测方法中:在步骤B中采用彩色CCD相机对零件照像前还有进行带通滤光的步骤。
进一步的,上述的基于激光三角测量法的亚像素表面形态检测方法中:所述的步骤C中:滤波采用N*N(N=3,5,7)的矩形模板。
进一步的,上述的基于激光三角测量法的亚像素表面形态检测方法中:所使用的激光器为红光激光器。
进一步的,上述的基于激光三角测量法的亚像素表面形态检测方法中:所述的红光激光器的波长为635~660nm。
本发明提出的一种基于激光三角测量法的亚像素表面形态检测方法,与已有的方法具有以下优点:
1、充分利用图像的彩色信息,在激光条纹是红色,背景是黑色的先验知识指导下,提取RGB图像的红色通道分量,减小噪声影响,提高高斯拟合条纹中心精度。
2、通过对同一位置多次拍摄叠加激光条纹图像提高精度,降低了对激光条纹宽度和亮度的要求,降低了激光器的成本。
3、通过先粗略确定条纹中心位置然后采用固定数量的像素拟合确定亚像素中心的方法避免了确定激光条纹边界而采用复杂算法,降低了算法复杂度。
下面通过结合具体实施例对本发明进行进一步的说明。
附图说明
附图1、三角测量法的示意图。
具体实施方式
本实施例是一种基于激光三角测量法的亚像素表面形态检测方法,包括如下步骤:
1、***安装如图1所示采用深圳华用科技型号为HV500GC的500万彩色像素相机,相机镜头焦距为35mm,在镜头前安装中心波长为650nm,带宽为40nm的带通滤光片,使用波长为650nm的激光器。使一字线激光器和水平工作台表面垂直,CCD摄像机主光轴与激光器在一个平面内,与水平工作台表面垂线成45~60度角。
图1中A、B 分别为工件表面处于不同高度时形成激光光斑的位置;A'、B'分别为A、B 经过CCD 镜头后成像之像点;θ为CCD 镜头光轴与垂直方向的夹角;ϕ为成像角;D 为激光光斑的位移量,即工件表面高度变化值;δ为激光光斑所成像点的位移量。
2、***安装好后,调节光圈,避免相机像素点过曝光。用相机对工件的同一位置连续拍摄5幅大小为的激光条纹图片,M、N分别表示图像的长和宽,每幅图片记为,x,y分别为图像的行和列的坐标分量。使用公式(1)将5幅激光条纹图像叠加成一副叠加图像。
(1)
3、叠加图像是彩色图像,每个像素点有RGB三个分量,这些分量可以组织成一个式2列向量形式:
(2)
其中是红色图像中像素亮度,其他两个分量表示绿色图像和蓝色图像的亮度。提取彩色的叠加图像的红色亮度分量得到得到红色通道图像。
4、采用3*3的方形模板对红色通道图像进行中值滤波,滤去噪声,得到中值滤波后的红色通道图像。
5、按列搜索中值滤波后的红色通道图像的每列像素的极大值,记录每列极大值位置(i=1,2……M),获得条纹中心的粗略位置。
6、采用固定数量的像素点进行高斯拟合求取条纹亚像素精度的中心,在每一列极大值位置两侧各选择K(5<K<10)个像素点进行高斯拟合求取条纹中心。具体实现过程如下:
对于第m列,令, 表示用来拟合高斯曲线的2K+1个像素点的灰度值。
假设待拟合的高斯灰度曲线为:
(4)
其中u为条纹的亮度中心坐标。
下面采用最小二乘法实现高斯曲线拟合。
对(4)式两边取对数,将指数分布转化成二次曲线分布:
(5)
令=,,,
则(5)式可以写为:
(6)
利用最小二乘法求出,激光条纹中心位置为:
(7)
本实施例中的步骤均在计算机中通过编程的方法实现,用MATLAB证明方法的有效性、正确性,***编程编程采用Visual C++ 和OpenCV具体实现。
Claims (6)
1.一种基于激光三角测量法的亚像素表面形态检测方法,该方法通过测量零件不同位置轮廓切片,叠加得到整个零件轮廓;其特征在于:包括以下步骤:
A、采用一字线激光器产生的激光照射零件;
B、采用彩色CCD相机对零件的同一位置采集2-10幅大小为M×N的激光条纹图片fi(x,y),x,y分别为图像的行和列的坐标分量;M、N分别表示图像的长和宽的像素数;并将采集到的所有激光条纹图片叠加形成一副叠加图像g(x,y),
s为彩色CCD相机对零件的同一位置采集的图片数;
C、提取叠加图像的红色通道图像r(x,y)并对其进行中值滤波,滤除噪声,得到中值滤波后的红色通道图像s(x,y);
D、沿进行中值滤波后的红色通道图像s(x,y)列方向搜索每一列像素极大值点;记录位置每列极大值位置ni,i=1,2……M,获得条纹中心的粗略位置;
E、在每一列像素极大值点ni两侧各取K,5<K<10个像素点进行高斯曲线拟合,提取亚像素精度的激光条纹中心线,形成亚像素表面形态图像;包括:
对于第m列,令f(x,m),x∈[ni-K,ni+K]表示用来拟合高斯曲线的2K+1个像素点的灰度值,
假设待拟合的高斯灰度曲线为:
其中u为条纹的亮度中心坐标,σ为标准偏差;
对待拟合的高斯灰度曲线式两边取对数,将指数分布转化成二次曲线分布:
令F(x,m)=lnf(x,m),
二次曲线分布可以写为:
F(x,m)=a0x2+a1x+a2
利用最小二乘法求出a0,a1,激光条纹中心位置为:
2.根据权利要求1所述的基于激光三角测量法的亚像素表面形态检测方法,其特征在于:在步骤A之前,还需要对测量设备进行校准和对测量***进行标定;
校准设备包括使一字线激光器和水平工作台表面垂直,CCD摄像机主光轴与激光器在一个平面内,与水平表面垂线成45~60度角;
对***进行标定包括得到真实空间高度与CCD拍摄图像高度之间的关系。
3.根据权利要求1所述的基于激光三角测量法的亚像素表面形态检测方法,其特征在于:在步骤B中采用彩色CCD相机对零件照像前还有进行带通滤光的步骤。
4.根据权利要求1所述的基于激光三角测量法的亚像素表面形态检测方法,其特征在于:所述的步骤C中:滤波采用N*N ,N=3,5,7的矩形模板。
5.根据权利要求1至4中任一所述的基于激光三角测量法的亚像素表面形态检测方法,其特征在于:所使用的激光器为红光激光器。
6.根据权利要求5所述的基于激光三角测量法的亚像素表面形态检测方法,其特征在于:所述的红光激光器的波长为635nm或者650nm或者660nm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210490109.9A CN103234475B (zh) | 2012-11-27 | 2012-11-27 | 一种基于激光三角测量法的亚像素表面形态检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210490109.9A CN103234475B (zh) | 2012-11-27 | 2012-11-27 | 一种基于激光三角测量法的亚像素表面形态检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103234475A CN103234475A (zh) | 2013-08-07 |
CN103234475B true CN103234475B (zh) | 2017-02-08 |
Family
ID=48882527
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210490109.9A Active CN103234475B (zh) | 2012-11-27 | 2012-11-27 | 一种基于激光三角测量法的亚像素表面形态检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103234475B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105241386B (zh) * | 2015-10-09 | 2018-01-05 | 中国水产科学研究院东海水产研究所 | 水下圆筒状网箱模型测量方法 |
CN106228542A (zh) * | 2016-07-13 | 2016-12-14 | 苏州光图智能科技有限公司 | 高速激光投影线峰值探测方法 |
CN110057325B (zh) * | 2019-04-26 | 2020-06-23 | 湖南大学 | 一种基于成像仿真的表面粗糙度检测方法及计算设备 |
CN111829457B (zh) * | 2020-07-09 | 2022-06-10 | 中国科学院光电技术研究所 | 基于结构光照明显微***的超薄膜器件三维形貌检测方法 |
CN111798519A (zh) * | 2020-07-21 | 2020-10-20 | 广东博智林机器人有限公司 | 激光条纹中心的提取方法、装置、电子设备及存储介质 |
CN112241964B (zh) * | 2020-09-22 | 2022-12-27 | 天津大学 | 一种用于线结构光非接触测量的光条中心提取方法 |
CN114266734B (zh) * | 2021-12-03 | 2023-04-07 | 深圳市桑瑞生物科技有限公司 | 植物生长阶段的识别方法、装置、设备及存储介质 |
CN116188459B (zh) * | 2023-04-23 | 2023-09-15 | 合肥金星智控科技股份有限公司 | 一种用于皮带撕裂检测的线激光快速识别方法及*** |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5905567A (en) * | 1994-03-25 | 1999-05-18 | Cae Newnes Ltd. | Method and apparatus for optimizing sub-pixel resolution in a triangulation based distance measuring device |
US6046812A (en) * | 1997-05-29 | 2000-04-04 | Korea Atomic Energy Research Institute | Shape-measuring laser apparatus using anisotropic magnification optics |
JP2000283753A (ja) * | 1999-03-31 | 2000-10-13 | Fuji Heavy Ind Ltd | ステレオ画像による測距装置 |
EP1510973A2 (en) * | 2003-08-29 | 2005-03-02 | Samsung Electronics Co., Ltd. | Method and apparatus for image-based photorealistic 3D face modeling |
CN101178812A (zh) * | 2007-12-10 | 2008-05-14 | 北京航空航天大学 | 一种结构光光条中心线提取的混合图像处理方法 |
CN102322822A (zh) * | 2011-08-08 | 2012-01-18 | 西安交通大学 | 一种三频彩色条纹投影三维测量方法 |
CN102572304A (zh) * | 2011-12-13 | 2012-07-11 | 广东威创视讯科技股份有限公司 | 一种图像叠加处理方法及装置 |
CN102495026B (zh) * | 2011-11-23 | 2013-08-28 | 天津大学 | 一种用于线激光扫描视觉测量***的光带中心线提取方法 |
-
2012
- 2012-11-27 CN CN201210490109.9A patent/CN103234475B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5905567A (en) * | 1994-03-25 | 1999-05-18 | Cae Newnes Ltd. | Method and apparatus for optimizing sub-pixel resolution in a triangulation based distance measuring device |
US6046812A (en) * | 1997-05-29 | 2000-04-04 | Korea Atomic Energy Research Institute | Shape-measuring laser apparatus using anisotropic magnification optics |
JP2000283753A (ja) * | 1999-03-31 | 2000-10-13 | Fuji Heavy Ind Ltd | ステレオ画像による測距装置 |
EP1510973A2 (en) * | 2003-08-29 | 2005-03-02 | Samsung Electronics Co., Ltd. | Method and apparatus for image-based photorealistic 3D face modeling |
CN101178812A (zh) * | 2007-12-10 | 2008-05-14 | 北京航空航天大学 | 一种结构光光条中心线提取的混合图像处理方法 |
CN102322822A (zh) * | 2011-08-08 | 2012-01-18 | 西安交通大学 | 一种三频彩色条纹投影三维测量方法 |
CN102495026B (zh) * | 2011-11-23 | 2013-08-28 | 天津大学 | 一种用于线激光扫描视觉测量***的光带中心线提取方法 |
CN102572304A (zh) * | 2011-12-13 | 2012-07-11 | 广东威创视讯科技股份有限公司 | 一种图像叠加处理方法及装置 |
Non-Patent Citations (3)
Title |
---|
实时采集多幅图像叠加平均去噪的研究与实现;由春秋;《哈尔滨师范大学自然科学学报》;20071231;第23卷(第3期);第32-35页 * |
彩色三维激光扫描成像***的研究;陶立;《中国优秀硕士学位论文全文数据库 信息科技辑》;20041215(第04期);第1140-507页 * |
激光三角法在物体三维轮廓测量中的应用;汤强晋;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20070415(第04期);第B002-119页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103234475A (zh) | 2013-08-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103234475B (zh) | 一种基于激光三角测量法的亚像素表面形态检测方法 | |
CN110322702B (zh) | 一种基于双目立体视觉***的车辆智能测速方法 | |
US9965870B2 (en) | Camera calibration method using a calibration target | |
CN105550670B (zh) | 一种目标物体动态跟踪与测量定位方法 | |
CN107133983B (zh) | 成捆圆钢端面双目视觉***与空间定位及计数方法 | |
TWI394431B (zh) | 評價立體影像顯示面板的方法及系統 | |
CN108760767A (zh) | 基于机器视觉的大尺寸液晶屏缺陷检测方法 | |
CN109919975B (zh) | 一种基于坐标标定的广域监控运动目标关联方法 | |
CN110472722A (zh) | 基于机器视觉技术的钢坯计数装置和方法 | |
CN108257137A (zh) | 一种基于视觉光斑自动判读的角度测量方法及*** | |
CN108833912A (zh) | 一种摄像机机芯光轴中心和视场角的测量方法及*** | |
CN107356202A (zh) | 一种激光扫描测量***目标自动照准方法 | |
CN109919856B (zh) | 基于双目视觉的沥青路面构造深度检测方法 | |
US11562478B2 (en) | Method and system for testing field of view | |
CN105139384B (zh) | 缺陷胶囊检测的方法和装置 | |
CN106404682A (zh) | 一种土壤颜色识别方法 | |
CN109341668A (zh) | 基于折射投影模型和光束追踪法的多相机测量方法 | |
CN102853777A (zh) | 基于图像处理的制动间隙宽度测量方法和设备 | |
CN107589551A (zh) | 一种多孔径偏振成像装置及*** | |
CN111179335A (zh) | 一种基于双目视觉的立木测定方法 | |
CN106709955A (zh) | 基于双目立体视觉的空间坐标系标定***和方法 | |
CN108288065A (zh) | 一种基于图像分析的四轮定位仪检测方法 | |
CN107360418A (zh) | 一种全景相机颜色一致性评判***及方法 | |
CN111121643A (zh) | 一种道路宽度的测量方法及*** | |
CN105391998B (zh) | 微光夜视仪分辨率自动检测方法和装置 |
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 |