CN114371168B - 一种基于反射照度计算的机器视觉照明设计方法 - Google Patents
一种基于反射照度计算的机器视觉照明设计方法 Download PDFInfo
- Publication number
- CN114371168B CN114371168B CN202111682069.3A CN202111682069A CN114371168B CN 114371168 B CN114371168 B CN 114371168B CN 202111682069 A CN202111682069 A CN 202111682069A CN 114371168 B CN114371168 B CN 114371168B
- Authority
- CN
- China
- Prior art keywords
- light source
- illumination
- led lamp
- reflection
- illuminance
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/84—Systems specially adapted for particular applications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/55—Specular reflectivity
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明属于机器视觉测量领域,提出了一种基于反射照度计算的机器视觉照明设计方法。该方法首先基于相机与同轴LED环形光源照明的几何位置关系计算包括光源距离、光源角度、数量、光源位置等参数的LED光源入射照度,然后基于双向光照反射模型计算了包含材料反射系数、入射角和反射角等参数的LED光源反射照度,给定光源的特征参数以及待测物体的特征参数,通过计算分析光源反射照度均匀度得到清晰成像时对应的光源最佳光源距离及照明角度等照明参数。该方法兼具高效、准确、通用性高等优势,对机器视觉的检测方法具有良好的应用价值。
Description
技术领域
本发明属于机器视觉测量领域,涉及一种基于反射照度计算的机器视觉照明设计方法。
背景技术
机器视觉测量技术因具有非接触、高精高效、测量设备体积小、成本低等优势,目前被广泛应用到航空航天、汽车、电子、计算机等重要工业制造领域。但在上述领域的实际应用中,受现场环境和待测零件表面反光特性不均的影响,致使在成像过程中难以采集清晰的图像特征。在机器视觉测量中,光源的合理布置和选择是高质量图像采集的重要保障。现有的光源照明设计方法主要是针对光源对物体照射的辐照度进行计算,通过入射照明的均匀度对照明结果进行评价,然而该方法计算参数中没有考虑不同待测物体本身的反射特征,此类方法对图像质量的提高有限,且不具备普适性。
针对机器视觉中光源照明的优化设计,桂林电子科技大学的甘勇等人于2021年在《应用光学》第42卷第4期发表了文章《光源角度对金属表面图像边缘的影响分析》,对光源的最佳照明角度进行了仿真分析与实验研究,通过改变光源角度,通过理论计算与仿真计算出了最佳的光源角度,并结合Sobel算法对最佳角度下的图像和其他角度下的图像边缘质量进行了对比。但是该方法采用了单一变量法,只改变了光源的照明角度,对于光源与相机的放置位置没有考虑,难以满足其他的照明设计方案。
发明内容
本发明为了弥补现有光源设计方法的不足,发明了一种基于反射照度计算的机器视觉照明设计方法,其目的是针对机器视觉中待测参数的不同特征,计算出同轴多角度环形LED阵列光源的最佳照明参数以实现高质量图像的获取,为基于机器视觉几何量参数的高精测量提供保障。首先基于相机与同轴LED环形光源照明的几何位置关系计算了包括光源距离、光源角度、数量、光源位置等参数的LED光源入射照度,然后基于双向光照反射模型Bidirectional Reflectance Distribution Function,BRDF光照反射模型计算了包含材料反射系数、入射角和反射角等参数的LED光源反射照度,给定光源的特征参数以及待测物体的特征参数,通过计算分析光源反射照度均匀度得到清晰成像时对应的光源最佳光源距离及照明角度等照明参数。该方法兼具高效、准确、通用性高等优势,对机器视觉的检测方法具有良好的应用价值。
本发明采用的技术方案是一种基于反射照度计算的视觉照明设计方法,该方法首先建立视觉测量系统的坐标系,相机1和待测物体4同轴,以相机1和待测物体4中心连线为z轴,以垂直于z轴的平面为xoy平面,以z轴与xoy平面交点为原点o;所设计的光源为多奇数层同心LED环形阵列,中间层环形阵列以上的称为上环形阵列,中间层环形阵列以下的称为下环形阵列,具体步骤如下:
第一步、计算LED光源入射照度;
单个的LED灯珠辐射照度值为:
式(1)中,L为LED灯珠的发光亮度,A为LED灯珠的发射面积,E(x,y,z)为辐射照度,(x,y,z)是被照明目标点的坐标,(xi,yi,zi)是单个LED的坐标,m由下式定义:
式(2)中,θ1/2是LED灯珠的半衰角,在环形光源阵列中,单个LED灯珠的坐标为:
式(3)中,R为LED灯珠所在位置到原点的距离,β为LED灯珠所在位置到原点的连线与z轴夹角,为LED灯珠与y轴夹角:
式(4)中,N为一个环形光源阵列总的LED灯珠个数,因此一个环形光源阵列的照度为所有LED灯珠照度之和的表达为:
式(5)中,n1为一分区,n2为LED发光的分区总数;当只有一个分区亮时,n1=n2=1,j为一个分区每个环形阵列对应的LED灯珠个数;根据第n层LED与中间第一层LED之间的几何关系,第n层LED灯珠的坐标为:
式(6)中,α=90°-β;n为上环形阵列距离中间环形阵列的层数,d为上下LED灯珠之间的距离,为上环形阵列中LED灯珠的坐标,/>下环形阵列中LED灯珠坐标,最终多层环形光源的入射照度值为:
式(6)中,Ej(x,y,z)为每个环形光源阵列的照度,En为所有环形光源阵列的照度;
第二步、计算待测物体特征的反射照度;
双向光照反射模型中的镜面反射光为:
式(8)中,Lo为光源的反射照度值,Ei为入射光的照度值,F为菲涅尔反射函数,D为微表面法线分布函数,G为几何衰减函数,θi为入射光线与表面法线的夹角,θo为表面法线与反射光线的夹角,其中D和G为微表面理论对应的函数,取1,F函数为下式:
F=RF(0°)+(1-RF(0°))(1-cosθi)5
式(9)中,其中η为材质的折射率;光源的反射照度为:
式(10)中,θi和θo的几何关系为:
式(11)中,α为待测表面与水平面的夹角,最终计算得到LED反射照度为:
第三步、计算反射照明均匀度;
反射照度平均值与最大值的比值即反射照明均匀度作为成像清晰的评价指标,反射照明均匀度为:
式(13)中,x∈[a,b],y∈[c,d]为待测表面的范围,m为x轴对应的行数,n为y轴对应的列数,乘积mn为待测表面所有待测点的总数;mn为待测表面所有待测点的总数;
第四步、计算最佳的照明参数;
根据待测物体4表面的特征参数和LED灯珠特征参数,计算不同光源距离R和照明角度β对应的所有反射照明均匀度值μ(x,y),最大μ(x,y)对应的R和β即为最佳的光源距离R和照明角度β。
要测量横刃长度5时,需要先照亮该表面即待照明表面8,接着采集图像后用图像边缘提取算法提取横刃边缘,最后进行长度的计算。
本发明的有益效果:本发明充分考虑了实际测量情况下待测物体表面反射率及光线入射和反射对相机成像的影响,有效引入了基于BRDF光照反射模型中的镜面反射函数。在给定待测表面的范围后,提出了一种基于反射照度计算的视觉照明设计方法,使用该方法设计光源,采集的图像经过边缘提取后,边缘清晰且连续,验证了本方法的有效性和准确性。本方法提出的基于反射照度计算的视觉照明设计方法可以准确快速地实现对光源最佳照明参数设计工作,对实际视觉测量中光源照明设计具有较好的应用价值。该方法是一种具有实际视觉测量应用价值的计算方法。计算方法简单高效,计算准确性高,通用性好。
附图说明
图1是本发明提出的视觉照明设计方法中基于反射照度的计算方法流程图;
图2(a)是LED同轴照明视觉测量系统示意图;
图2(b)是LED灯珠位置示意图;
图3是同轴多角度环形阵列光源LED几何位置示意图。
图4是使用本发明提出的照明设计方法计算出的照明参数后采集到的图像。
图5是使用本发明提出的照明设计方法经过边缘处理后提取到的待测物体尺寸测量所需的边缘信息。
图中:1为相机,2为镜头,3为LED光源,4为待测物体,5为横刃长度,6为钻芯厚度,7为微齿槽深度,8为待照明表面。
具体实施方式
下面结合附图和实施例对本发明作进一步阐述。
本实施例选用微齿钻头端面参数视觉测量中对横刃长度的图像采集时半衰角为θ1/2=30°的LED照明参数进行计算。
其中,每个半衰角为θ1/2=30°的LED灯珠的选用参数为:LED的发光亮度L=10000Wm-2sr-1,LED的发光面积A=πmm2,LED环形阵列总个数为N=120个,分区为8个分区,只选用一分区n1=n2=1,一分区的每个环形阵列的LED灯珠个数为j=9个,每个环形阵列中相邻LED灯珠的夹角为第一个LED灯珠与y轴的夹角/>总共的环形阵列层数n=11层,每层LED间隔距离为d=2mm。光源到坐标原点的距离首先取最小值R=45mm,LED光线到原点的连线与竖直方向的夹角取最值β=30°。
待测微齿钻头的横刃范围x∈[0,3]mm,y∈[0,2]mm,z=0mm,横刃范围所有待测点以0.1mm间隔的总个数为mn=651个,钻头的材料的折射率近似为η=0.71,横刃表面与xoy平面的夹角近似为α=15°。
如图1所示视觉照明设计方法流程图,首先按照图2建立系统坐标系,接着把上述物体特征参数作为输入量,LED特征参数作为常量,接着使用本发明提出的基于反射照度计算的视觉照明设计方法,计算方法的具体步骤如下;
第一步、计算LED光源入射照度;
将给定的R,β,代入式(3)首先计算一个环形阵列中一个LED灯珠坐标为xi=22.50mm,yi=0mm,zi=38.97mm,接着将θ1/2代入式(2)得m=3,则第一个LED灯珠到点(0,0,0)的LED入射照明函数值根据式(1)计算得E=5.793×10-1w/mm2,第二个LED灯珠与第一个LED灯珠的夹角根据式(4)计算为/>到第9个LED灯珠夹角计算为/>将上述每个LED灯珠辐射照度值代入式(5)得Ej=0.542w/mm2,根据式(6)计算上层LED灯珠坐标为最终一分区多层LED光源环形总入射辐照度利用式(7)为En(0,0,0)=5.214w/mm2。重复上述计算过程计算651个点的入射辐照度。
第二步、计算LED光源反射照度
首先计算点(0,0,0)的反射照度值,根据式(11)计算得β0=54.65°,α=15°,θi=39.65°,接着利用式(9)计算菲涅尔系数RF(39.65°)=0.71,根据菲涅尔系数和式(8)计算镜面反射系数fs=0.2396,最后计算LED反射照度值Lo=1.443×10-2w/mm2。重复上述计算过程计算651个点的反射照度值。最后根据式(12)计算整个待测范围内的反射照度值为
第三步、计算反射照明均匀度
根据上述计算结果和式(13)可以计算得到R=45,β=30°时反射照明均匀度μ=0.7730,重复计算上述步骤,计算R∈[45,150],β∈[15°,75°]对应的不同的反射照明均匀度值,最后利用式(14)找到最佳的LED照明参数为R=100mm,β=34.38°。
在给定待测物体表面范围后,本发明提出的基于反射照度计算的视觉照明设计方法,对光源参数进行了设计,使用设计后的光源进行图像采集,后经过图像边缘提取后的提取的钻头横刃结果如图5所示。从图中可以看出,本发明提出的基于反射照度计算的视觉照明设计方法得到的图像及边缘检测图清晰且连续,满足后续测量实际长度的需要,证明了本发明所提出的方法是正确和有效的。
本发明提出了一种基于反射照度计算的视觉照明设计方法,并充分考虑了真实测量情况下待测物体表面反射率及光线入射和反射对相机成像的影响,有效引入了基于BRDF光照反射模型中的镜面反射函数。在给定待测物体表面的范围后,本发明提出的基于反射照度计算的视觉照明设计方法可以得到视觉测量中光源的最佳照明参数,进一步使用设计后的光源可以得到特征边缘清晰的图像,对实际视觉测量中光源的设计具有较好的应用价值。
Claims (1)
1.一种基于反射照度计算的机器视觉照明设计方法,其特征在于,该方法首先建立视觉测量系统的坐标系,相机(1)和待测物体(4)同轴,以相机(1)和待测物体(4)中心连线为z轴,以垂直于z轴的平面为xoy平面,以z轴与xoy平面交点为原点o;所设计的光源为多奇数层同心LED环形阵列,中间层环形阵列以上的称为上环形阵列,中间层环形阵列以下的称为下环形阵列,具体步骤如下:
第一步、计算LED光源入射照度;
单个的LED灯珠辐射照度值为:
式(1)中,L为LED灯珠的发光亮度,A为LED灯珠的发射面积,E(x,y,z)为辐射照度,(x,y,z)是被照明目标点的坐标,(xi,yi,zi)是单个LED的坐标,m由下式定义:
式(2)中,θ1/2是LED灯珠的半衰角,在环形光源阵列中,单个LED灯珠的坐标为:
zi=Rcosβ
式(3)中,R为LED灯珠所在位置到原点的距离,β为LED灯珠所在位置到原点的连线与z轴夹角,为LED灯珠与y轴夹角:
式(4)中,N为一个环形光源阵列总的LED灯珠个数,因此一个环形光源阵列的照度为所有LED灯珠照度之和的表达为:
式(5)中,n1为一分区,n2为LED发光的分区总数;当只有一个分区亮时,n1=n2=1,j为一个分区每个环形阵列对应的LED灯珠个数;根据第n层LED与中间第一层LED之间的几何关系,第n层LED灯珠的坐标为:
式(6)中,α=90°-β;n为上环形阵列距离中间环形阵列的层数,d为上下LED灯珠之间的距离,为上环形阵列中LED灯珠的坐标,/>下环形阵列中LED灯珠坐标,最终多层环形光源的入射照度值为:
式(6)中,Ej(x,y,z)为每个环形光源阵列的照度,En为所有环形光源阵列的照度;
第二步、计算待测物体特征的反射照度;
双向光照反射模型中的镜面反射光为:
式(8)中,Lo为光源的反射照度值,Ei为入射光的照度值,F为菲涅尔反射函数,D为微表面法线分布函数,G为几何衰减函数,θi为入射光线与表面法线的夹角,θo为表面法线与反射光线的夹角,其中D和G为微表面理论对应的函数,取1,F函数为下式:
F=RF(0°)+(1-RF(0°))(1-cosθi)5
式(9)中,其中η为材质的折射率;光源的反射照度为:
式(10)中,θi和θo的几何关系为:
式(11)中,α为待测表面与水平面的夹角,最终计算得到LED反射照度为:
第三步、计算反射照明均匀度;
反射照度平均值与最大值的比值即反射照明均匀度作为成像清晰的评价指标,反射照明均匀度为:
式(13)中,x∈[a,b],y∈[c,d]为待测表面的范围,m为x轴对应的行数,n为y轴对应的列数,乘积mn为待测表面所有待测点的总数;mn为待测表面所有待测点的总数;
第四步、计算最佳的照明参数;
根据待测物体(4)表面的特征参数和LED灯珠特征参数,计算不同光源距离R和照明角度β对应的所有反射照明均匀度值μ(x,y),最大μ(x,y)对应的R和β即为最佳的光源距离R和照明角度β。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111682069.3A CN114371168B (zh) | 2021-12-31 | 2021-12-31 | 一种基于反射照度计算的机器视觉照明设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111682069.3A CN114371168B (zh) | 2021-12-31 | 2021-12-31 | 一种基于反射照度计算的机器视觉照明设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114371168A CN114371168A (zh) | 2022-04-19 |
CN114371168B true CN114371168B (zh) | 2023-10-13 |
Family
ID=81142854
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111682069.3A Active CN114371168B (zh) | 2021-12-31 | 2021-12-31 | 一种基于反射照度计算的机器视觉照明设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114371168B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5764874A (en) * | 1994-10-31 | 1998-06-09 | Northeast Robotics, Inc. | Imaging system utilizing both diffuse and specular reflection characteristics |
CN1769834A (zh) * | 2004-10-21 | 2006-05-10 | 株式会社米姿托约 | 用于控制一致精确和高速检查视觉系统的模糊限制基系统和方法 |
WO2013111134A1 (en) * | 2012-01-24 | 2013-08-01 | Yechezkal Evan Spero | Detector controlled illuminating system |
CN103857982A (zh) * | 2011-08-15 | 2014-06-11 | 纽约市哥伦比亚大学理事会 | 使用漫射结构光实现机器视觉的系统和方法 |
CN109901344A (zh) * | 2019-04-01 | 2019-06-18 | 合刃科技(武汉)有限公司 | 环形同轴光源产生装置、机器视觉系统及应用方法 |
CN110543741A (zh) * | 2019-09-20 | 2019-12-06 | 中南大学 | 一种用于视觉检测的直接照明光源设计方法 |
CN110596122A (zh) * | 2019-09-20 | 2019-12-20 | 中南大学 | 一种用于视觉检测的漫反射光源系统及设计方法 |
-
2021
- 2021-12-31 CN CN202111682069.3A patent/CN114371168B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5764874A (en) * | 1994-10-31 | 1998-06-09 | Northeast Robotics, Inc. | Imaging system utilizing both diffuse and specular reflection characteristics |
CN1769834A (zh) * | 2004-10-21 | 2006-05-10 | 株式会社米姿托约 | 用于控制一致精确和高速检查视觉系统的模糊限制基系统和方法 |
CN103857982A (zh) * | 2011-08-15 | 2014-06-11 | 纽约市哥伦比亚大学理事会 | 使用漫射结构光实现机器视觉的系统和方法 |
WO2013111134A1 (en) * | 2012-01-24 | 2013-08-01 | Yechezkal Evan Spero | Detector controlled illuminating system |
CN109901344A (zh) * | 2019-04-01 | 2019-06-18 | 合刃科技(武汉)有限公司 | 环形同轴光源产生装置、机器视觉系统及应用方法 |
CN110543741A (zh) * | 2019-09-20 | 2019-12-06 | 中南大学 | 一种用于视觉检测的直接照明光源设计方法 |
CN110596122A (zh) * | 2019-09-20 | 2019-12-20 | 中南大学 | 一种用于视觉检测的漫反射光源系统及设计方法 |
Non-Patent Citations (3)
Title |
---|
"光源角度对金属表面图像边缘的影响分析";甘勇等;《应用光学》;第42卷(第4期);全文 * |
千分尺检定装置照明系统的设计;孔明;徐志玲;陈旺达;厉志飞;金骥;;机床与液压(第19期);全文 * |
基于光照约束的AGV视觉导引非均匀光照增强方法;杨旭;楼佩煌;武星;喻俊;;机械科学与技术(第12期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114371168A (zh) | 2022-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9370799B2 (en) | Method and system for optically inspecting a manufactured part at a single inspection station having a measurement axis | |
US9575013B2 (en) | Non-contact method and system for inspecting a manufactured part at an inspection station having a measurement axis | |
US9697596B2 (en) | Method and system for optically inspecting parts | |
EP2252856B1 (en) | Method of determination of glass surface shape and optical distortion by reflected optical imaging | |
CN108760765A (zh) | 一种基于侧视相机拍摄的表面损伤缺陷检测装置及方法 | |
CN102435512B (zh) | 一种板料成形极限曲线的视觉测量方法 | |
US20130258046A1 (en) | Method and system for optically inspecting outer peripheral surfaces of parts | |
US20140346094A1 (en) | High-speed, triangulation-based, 3-d method and system for inspecting manufactured parts and sorting the inspected parts | |
CN102410811A (zh) | 一种弯管参数的测量方法和系统 | |
CN102853786A (zh) | 平整度检测装置和方法 | |
CN105424726A (zh) | 基于机器视觉的发光面板检测方法 | |
US20170307541A1 (en) | Method and system for optically inspecting headed manufactured parts | |
CN1808056A (zh) | 利用二维检测器测量薄膜特性的装置及测量方法 | |
CN103246007A (zh) | 导光装置、前光模块以及反射式显示器 | |
CN108346141A (zh) | 单边侧入光式导光板缺陷提取方法 | |
US20170307538A1 (en) | Method and system for optically inspecting headed manufactured parts | |
CN114371168B (zh) | 一种基于反射照度计算的机器视觉照明设计方法 | |
CN115375636A (zh) | 一种动力电池模组全尺寸检测方法及设备 | |
CN104848808A (zh) | 一种表面粗糙度检测方法和设备 | |
CN201811716U (zh) | 陶瓷板检测设备 | |
CN109345521B (zh) | 一种基于皮革表面缺陷检测的led自由曲面阵列 | |
CN114414487A (zh) | 一种二维与三维融合的圆钢表面成像系统 | |
Reynolds et al. | Theory and applications of a surface inspection technique using double-pass retroreflection | |
CN204027528U (zh) | 一种视觉检测装置 | |
CN102003941B (zh) | 基于视觉的大型系泊链五环长测量方法及其测量装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |