CN114565565A - 一种视觉测量目标中心亚像素定位方法 - Google Patents
一种视觉测量目标中心亚像素定位方法 Download PDFInfo
- Publication number
- CN114565565A CN114565565A CN202210128326.7A CN202210128326A CN114565565A CN 114565565 A CN114565565 A CN 114565565A CN 202210128326 A CN202210128326 A CN 202210128326A CN 114565565 A CN114565565 A CN 114565565A
- Authority
- CN
- China
- Prior art keywords
- pixel
- formula
- edge
- center
- target
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000005259 measurement Methods 0.000 title claims abstract description 30
- 238000004364 calculation method Methods 0.000 claims abstract description 32
- 238000012216 screening Methods 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000005457 optimization Methods 0.000 claims description 5
- 238000013329 compounding Methods 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 235000021384 green leafy vegetables Nutrition 0.000 claims description 3
- 238000004613 tight binding model Methods 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 abstract description 15
- 238000003384 imaging method Methods 0.000 abstract description 11
- 230000000007 visual effect Effects 0.000 abstract description 4
- 238000009826 distribution Methods 0.000 description 7
- 238000004140 cleaning Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 241000974840 Ellipes Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000000499 gel Substances 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20021—Dividing image into blocks, subimages or windows
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20092—Interactive image processing based on input by user
- G06T2207/20104—Interactive definition of region of interest [ROI]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computing Systems (AREA)
- Quality & Reliability (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种视觉测量目标中心亚像素定位方法,即一种边缘模糊中心过饱和状态的目标光斑高精度中心定位方法,解决了离焦成像系统无法使用边缘拟合中心定位算法,以及高斯曲面拟合无法准确拟合过饱和光斑的问题。本发明通过对目标光斑边缘区域准确定位并筛选,得到中心过饱和光斑的有效计算区域,采用高斯权重函数与距离权重函数相复合,得到每个有效像素对中心定位提供的影响因子,根据加权质心定位算法的思想,结合像素影响因子计算得到目标中心定位。本发明精度高,抗干扰能力强,受噪声影响小,适用于多种视觉测量环境。
Description
技术领域
本发明涉及计算机视觉、视觉测量技术领域,具体为一种视觉测量目标中心亚像素定位方法,实现了针对目标边缘模糊且亮度饱和情况下的高精度中心定位。
背景技术
目标中心定位算法被广泛应用于机器视觉中的三维重建、视觉测量和航空航天器姿态测量等诸多领域。目标定位的准确程度对最终测量结果起到了决定性的作用。在视觉测量中,常用的目标有以下三种:红外主动光源、人工逆反射标志和被测物形貌特征点。红外主动光源以其独有的稳定性和良好的抗干扰能力,被广泛的使用于各种视觉测量系统中。本发明主要针对红外光源图像目标中心定位时遇到的亮度饱和问题,提供一种基于目标光斑边缘区域计算目标中心亚像素坐标的方法。
当前,视觉测量领域的目标识别主要分为两大类。
第一类是基于边缘信息的椭圆拟合方法,该方法适用于对焦成像系统中的人工逆反射标志目标中心定位。这一类测量系统实时性差,对于10米以上的大型工件只能采用离线测量,且逆反射标志受环境光影响严重,不适于户外测量。基于边缘的中心定位算法在实际应用中存在以下不足:1)测量前的准备工作及测量完成后的清理工作复杂;由于使用人工逆反射标志,所以在测量前要将反射标志粘贴于被测工件表面,为了保证测量精度,往往与要大量的反光标志。这同时为测量完成的清理工作带来了很多麻烦。2)应用场景受限;若要达到工业测量级别的精度,就需要保证目标边缘成像清晰且稳定,这就使得测量只能在室内相对稳定的环境中进行。3)测量过程非常耗时;为了保证成像清晰,拍摄距离通常保持在3至5米范围,对于大型工件就需要测量人员从不同角度拍摄成百上千的图像,以得到理想精度的测量结果。4)抗干扰能力差;人工逆反射标志的成像受环境光影响严重,并且拍摄角度过大也会导致成像不清晰。
第二类是基于目标灰度分布进行拟合,主要包括灰度质心法、加权灰度法、高斯曲面拟合法等算法。基于灰度分布的方法多在离焦成像系统中使用,所用的目标为红外主动光源。这一类测量系统可在线测量超大型工件,且能够应对复杂的室外测量环境。基于灰度分布的中心定位算法也存在这以下不足:1)目标光斑亮度难统一;测量过程中不同目标之间距离远,导致相同曝光下成像的亮度不同,容易出现远处目标亮度太小的同时近处目标亮度过高导致成像过曝,从而使得目标中心定位不准确。2)灰度分布要求高;基于灰度的方法在定位时所有的信息均来自像素灰度分布,同一个目标,灰度分布轻微变化有可能导致中心定位结果差别非常大。
发明内容
本发明目的是提供一种边缘模糊中心过饱和状态的目标光斑高精度中心定位方法,用于解决离焦成像系统无法使用边缘拟合中心定位算法,以及高斯曲面拟合无法准确拟合过饱和光斑的问题。
本发明是采用如下技术方案实现的:
一种视觉测量目标中心亚像素定位方法,包括如下步骤:
步骤一、确定目标光斑兴趣区域。
步骤二、采用双调和插值方法对兴趣区域进行像素细分,所采用的插值公式为:
Ip=GpW (1)
其中,下标p为所求插值点编号;式(1)中Gp为1×k阶格林函数矩阵,计算公式为:
Gp=[d01 d02 ... d0k] (2)
式(2)中的元素d为插值核函数,具体计算公式为:
W=G-1Z (4)
其中,Z为k×1阶属性向量,由所有像素值组成Z=[z1 z2 ... zk]T
G为k×k格林函数矩阵,如:
式(5)矩阵中的元素d计算公式如式(3);
将式(2)至式(5)代入式(1)中计算出插值点的像素值。
步骤三、首先使用Sobel算子得到所有像素x方向和y方向的梯度和幅值;若某像素梯度幅值在其法向上为极大值,则判断该像素为边缘像素,否则置零;记边缘像素坐标为(xi,yi),i=1,...,k,然后根据如下公式计算每个边缘像素的补偿值:
其中,i为选定的边缘像素索引号,|gx|和|gy|表示x方向和y方向梯度值,||g(x,y)||为像素(x,y)的梯度幅值;将式(6)得到补偿值与边缘像素坐标相加得到亚像素边缘坐标(xi+ηx,yi+ηy),这些点是目标光斑边缘区域梯度变化最大的位置。
步骤四、对步骤三中得到的边缘点椭圆拟合,然后选取[0.9rmax,1.1rmax]范围内的像素组成待选像素集合{PixelTBS}。
步骤五、对目标光斑兴趣区域的灰度值进行灰度投票统计,得到灰度值统计数列{GrayStatic}。
步骤六、对步骤五中得到数列{GrayStatic}求导,并根据如下约束选取最大统计值对应的灰度值Imax1和第二大统计值对应的灰度值Imax2:
|Imax1-Imax2|>150 (7)。
步骤七、在步骤六筛选得到的Imax1和Imax2附近找统计数量百分比大于1%的灰度值,分别与Imax1、Imax2组成背景阈值数列{Th_Background_GrayValues}和中心区域阈值数列{Th_Center_GrayValues};采用如下计算公式得到背景阈值和中心区域阈值:
步骤八、对步骤四中得到的待选像素集合{PixelTBS}使用步骤七得到的两个阈值进行筛选,所有在阈值范围内的像素(Th_Ibg<Im<Th_Icen)组成集合{Valid_Pixel},参与后续中心坐标计算,而不在阈值范围内的像素灰度值置零。
步骤九、利用步骤八中有效像素{Valid_Pixel}与步骤三中使用公式(6)计算得到的边缘点在梯度方向的距离计算每个有效像素的距离权重因子Wdi(i=1,...,k);距离权重因子的计算公式为:
其中,dist为像素i到亚像素边缘的距离。
步骤十、采用二维离散高斯函数对步骤八得到的有效像素集合{Valid_Pixel}拟合,根据得到的拟合参数计算每个有效像素的灰度权重因子Wgi(i=1,...,k);采用的二维高斯拟合函数为:
使用式(11)作为目标函数进行非线性优化,求得Wgi(i=1,...,k)最优解;具体计算公式为:
其中,Ii为像素i=1,...,k的灰度值。
步骤十一、将像素灰度值Ii和步骤十计算的灰度权重因子Wgi进行复合,得到新的影响因子λi,复合公式为:
步骤十二、将步骤十一计算得出的影响因子代入式(13)计算目标光斑中心位置(xc,yc):
本发明使用高精度双调和插值算法对目标光斑进行高精度快速插值,将光斑边缘定位和中心定位提升至亚像素级别,通过对目标光斑边缘区域准确定位并筛选,得到中心过饱和光斑的有效计算区域,采用高斯权重函数与距离权重函数相复合,得到每个有效像素对中心定位提供的影响因子,根据加权质心定位算法的思想,结合像素影响因子可以计算得到更接近真实目标中心的定位坐标,解决了由于目标亮度差别大,部分光斑中心存在过饱和像素,从而导致目标中心定位存在截断误差的问题。
本发明设计合理,对于大直径光斑的中心定位精度明显优于高斯曲面拟合算法,且抗干扰能力强受噪声影响较小,鲁棒性较好,适用于多种视觉测量环境。
附图说明
图1表示本发明实施例中的光斑兴趣区域。
图2表示本发明实施例中的Sobel算子。
图3表示本发明实施例与现有技术的结果对比示意图。
具体实施方式
以下结合一个具体实例对本发明进行详细解释,并将现有文献[2]和高斯曲面拟合法做横向比较以说明本发明的有益效果。
本发明实施例使用仿真光斑。仿真相机参数选用Basler A102f工业相机,分辨率为1392×1040,镜头焦距f=12mm。使用MATLAB科学计算软件做相关算法验证。具体步骤如下:
一、确定目标光斑兴趣区域,所用方法包括但不限于图像滤波、阈值分割、目标形状筛选等方法,最终获得的目标光斑兴趣区域如图1所示。
本实施例具体为:人工设定8个光斑中心名义值,采用高斯光斑模型结合离焦成像模型在同一幅图像上生成8个目标光斑,在图像上添加信噪比为50dB的高斯噪声,重复添加100组,生成100组图像,显而易见,每组图像都包含8个仿真目标光斑。将图像分别进行二值化后,根据matlab的连通域信息确定8个光斑大致的中心坐标
二、对于一个目标光斑,其包含的灰度信息十分有限,为了最大程度利用目标光斑的灰度信息,使用像素细分技术对两个相邻像素进行灰度值插值,尽可能减小量化误差。本发明使用双调和插值方法(具体参考文献[1]:Deng X,Tang Z.Moving surface splineinterpolation based on Green’s function[J].Mathematical geosciences,2011,43(6):663-680.)对兴趣区域进行像素细分,该方法对比广泛使用的三次样条插值具有更高的插值精度和更短的计算时间。所采用的插值公式为:
Ip=GpW (1)
下标p为所求插值点编号。
式(1)中Gp为1×k阶格林函数矩阵,计算公式为:
Gp=[d01 d02 ... d0k] (2)
式(2)中的元素d为插值核函数,具体计算公式为:
W=G-1Z (4)
Z为k×1阶属性向量,由所有像素值组成Z=[z1 z2 ... zk]T;
G为k×k格林函数矩阵,形如:
矩阵中的元素d计算公式如式(3);
将式(2)至式(5)代入式(1)可以计算出插值点的像素值。
本实施例具体为:对于第p张图象的第c个目标光斑,选取每个坐标周围(31×31)pixel的区域作为计算窗口window,将窗口左上角像素坐标记为该坐标值为目标光斑中心坐标的全局偏移量。对窗口window内的光斑像素灰度值双调和插值,得到(121×121)pixel=14641pixel的待计算区域roi_interp。
三、为了选取合适的边缘像素进行中心定位,需要对光斑目标进行亚像素边缘提取,确定有效边缘区域。现有技术文献[2](石达顺,唐朝.基于统计高斯拟合的圆形光斑中心定位方法[J].测控技术,2020,39(7):51-56.)使用Canny算法提取圆形光斑边缘,但由于本发明所针对目标光斑为离焦成像,无法取得清晰的边缘成像,故使用文献[2]中的方法会导致边缘偏差较大。本发明首先使用Sobel算子得到所有像素x方向和y方向的梯度和幅值,图2例举了竖直方向和水平方向的Sobel算子。若某像素梯度幅值在其法向上为极大值,则判断该像素为边缘像素,否则置零。记边缘像素坐标为(xi,yi),i=1,...,k,然后根据如下公式计算每个边缘像素的补偿值:
其中,i为选定的边缘像素索引号,|gx|和|gy|表示x方向和y方向梯度值,||g(x,y)||为像素(x,y)的梯度幅值。将式(6)得到补偿值与边缘像素坐标相加得到亚像素边缘坐标(xi+ηx,yi+ηy),这些点是目标光斑边缘区域梯度变化最大的位置。
本实施例具体为:使用Sobel算子计算roi_interp内所有像素的梯度幅值,取像素梯度方向的极大值作为边缘像素坐标整数部分,记为其中ei=1,K,Nedge表示该区域内共有Nedge个像素是边缘像素。对于边缘像素ei,根据公式(6)计算其补偿值(ηx,ηy),将补偿值与边缘像素坐标整数部分相加,得到最终的亚像素边缘坐标
四、对步骤三中得到的边缘点椭圆拟合,然后选取[0.9rmax,1.1rmax]范围内的像素组成待选像素集合{PixelTBS}。
本实施例具体为:使用最小二乘拟合方法拟合区域roi_interp内所有边缘点,拟合函数为椭圆函数ax2+by2+cxy+dx+ey+f=0,得到椭圆的参数(a,b,c,d,e,f)从而计算出中心(xellip,yellip)和长轴半径rmax,计算区域roi_interp内所有像素与椭圆中心(xellip,yellip)的距离选择|dj-rmax|<0.1rmax的像素放入待选像素集合{PixelTBS}。
五、对目标光斑兴趣区域的灰度值进行灰度投票统计,得到灰度值统计数列{GrayStatic}。
六、对步骤五中得到数列{GrayStatic}求导,并根据如下约束选取最大统计值对应的灰度值Imax1和第二大统计值对应的灰度值Imax2:
|Imax1-Imax2|>150 (7)
七、在步骤六筛选得到的Imax1和Imax2附近找统计数量百分比大于1%的灰度值,分别与Imax1、Imax2组成背景阈值数列{Th_Background_GrayValues}和中心区域阈值数列{Th_Center_GrayValues}。采用如下计算公式得到背景阈值和中心区域阈值:
本实施例具体为:统计roi_interp内像素的灰度值数量,由于图像传感器为8位,所以灰度值为0—255的整数,统计间隔1灰度值,得到数列{GrayStatic}={N1,...,N255}。降序排列{GrayStatic}将灰度值计数最多的灰度值Imax1记为第一个阈值,在剩余的灰度值计数中寻找满足|Imax1-Imax2|>150的第二个灰度值计数最多的灰度值Imax2,假设Imax1<Imax2。挑选灰度值计数Nvl>0.01×14641≈146,(vl=1,...,255)的所有灰度值Ivl,将|Imax1-Ivl|<20的灰度值Ivl与Imax1一起记为背景阈值数列{Th_Background_GrayValues},将|Imax2-Ivl|<20的灰度值Ivl与Imax2一起记为中心区域阈值数列{Th_Center_GrayValues}。根据式(8)寻找背景阈值和中心阈值。
八、对步骤四中得到的待选像素集合{PixelTBS}使用步骤七得到的两个阈值进行筛选,所有在阈值范围内的像素(Th_Ibg<Im<Th_Icen)组成集合{Valid_Pixel},参与后续中心坐标计算,而不在阈值范围内的像素灰度值置零。
本实施例具体为:对于所有(xj,yj)∈{PixelTBS},如果Th_Ibg<I(xj,yj)<Th_Icen,则将像素(xj,yj)记为有效像素(xsj,ysj),组成有效像素集合{Valid_Pixel}。
九、利用步骤八中有效像素{Valid_Pixel}与步骤三中使用公式(6)计算得到的边缘点在梯度方向的距离计算每个有效像素的距离权重因子Wdi(i=1,...,k)。主要为了消除以下两个因素的影响:第一,环境中不可避免的杂散光,以及测量设备电子器件产生的噪声等因素的影响,目标光斑中比较接近背景区域的像素受噪声影响较大,在计算目标中心时,这些像素所占权重应较小;第二,目标光斑中心区域像素过饱和,由于感光元器件互相之间存在“串扰”现象,饱和像素将对相邻像素产生一定影响,故这些受影响像素参与中心计算的权重也应更小。距离权重因子的计算公式为:
dist为像素i到亚像素边缘的距离。
本实施例具体为:对于有效像素集合{Valid_Pixel}中的像素(xsj,ysj),连接(xsj,ysj)与椭圆中心(xellip,yellip)该线段(或线段延长线)交椭圆圆周于点(xint,yint),计算(xint,yint)与(xsj,ysj)之间的距离得到像素(xsj,ysj)的距离权重因子:
十、采用二维离散高斯函数对步骤八得到的有效像素集合{Valid_Pixel}拟合,现有技术文献[2]直接使用高斯函数拟合灰度分布来计算目标光斑中心,但其所用光斑直径较大且边缘清晰,针对本发明所使用的离焦光斑效果不佳。本发明根据得到的拟合参数可以计算每个有效像素的灰度权重因子Wgi(i=1,...,k),一方面避免现有技术文献[3](郝颖明,束枫等.目标中心定位方法对位姿测量结果的影响分析[J].仪器仪表学报,2007,28(4):1-6.)中指出的“高斯曲面拟合对大直径目标光斑精度及鲁棒性较差”的问题,另一方面可滤除随机噪声对后续中心定位的影响。采用的二维高斯拟合函数为:
使用式(11)作为目标函数进行非线性优化,可求得Wgi(i=1,...,k)最优解,由于步骤九中提到的原因,本发明在优化目标函数中加入距离权值因子,保证梯度越大的像素对拟合结果影响越大,越靠近边缘和饱和位置的像素对拟合结果影响越小,具体计算公式为:
其中,Ii为像素i=1,...,k的灰度值。
本实施例具体为:使用二维高斯函数对效像素集合{Valid_Pixel}中的像素进行曲面拟合,二维高斯函数为:其中Wgsi为像素(xsj,ysj)的所求灰度权重因子,设参数初值为1,采用LM非线性优化方法求取参数最优值。目标函数为:其中k为有效像素个数,Wdsj为计算得到的像素(xsj,ysj)的距离权重因子。
十一、采用加权质心法的思想,将像素灰度值Ii和步骤十计算的灰度权重因子Wgi进行复合,得到新的影响因子λi,复合公式为:
十三、将步骤十二的结果存入寄存器,然后对当前图像的下一个目标光斑重复步骤二至步骤十二,确定下一个目标光斑的图像坐标,当全部完成8个目标光斑的中心定位后按同样的方法对其余图像上的目标光斑进行中心定位,直到100组图像全部完成。
对上述过程在MATLAB软件中编程实现算法验证,完成100组图像的总用时为11.44秒。相比之下,高斯曲面拟合算法总用时为12.52秒,文献[2]提供的方法总用时为15.41秒,可以发现三种方法总用时差别不大。
为了验证本发明方法的实施精度,将三种方法对同一个图像光斑的定位结果展示于图3,光斑名义值以黑色十字表示,本发明所使用的目标光斑中心定位准确度明显优于其他两种方法,且100结果的重复性也更好。
总之,本发明通过对目标光斑边缘区域准确定位并筛选,得到中心过饱和光斑的有效计算区域,采用高斯权重函数与距离权重函数相复合,得到每个有效像素对中心定位提供的影响因子,根据加权质心定位算法的思想,结合像素影响因子计算得到目标中心定位。本发明精度高,抗干扰能力强,受噪声影响小,适用于多种视觉测量环境。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照本发明实施例进行了详细说明,本领域的普通技术人员应当理解,本发明方法的保护范围不局限应用于红外靶灯识别,基于本方法的原理发明的可控光源的识别方法,无论是等同变换或改进,均不应排除在本发明的保护范围之外。
Claims (1)
1.一种视觉测量目标中心亚像素定位方法,其特征在于:包括如下步骤:
步骤一、确定目标光斑兴趣区域;
步骤二、采用双调和插值方法对兴趣区域进行像素细分,所采用的插值公式为:
Ip=GpW (1)
其中,下标p为所求插值点编号;式(1)中Gp为1×k阶格林函数矩阵,计算公式为:
Gp=[d01 d02...d0k] (2)
式(2)中的元素d为插值核函数,具体计算公式为:
W=G-1Z (4)
其中,Z为k×1阶属性向量,由所有像素值组成Z=[z1 z2...zk]TG为k×k格林函数矩阵,如:
式(5)矩阵中的元素d计算公式如式(3);
将式(2)至式(5)代入式(1)中计算出插值点的像素值;
步骤三、使用Sobel算子得到所有像素x方向和y方向的梯度和幅值;若某像素梯度幅值在其法向上为极大值,则判断该像素为边缘像素,否则置零;记边缘像素坐标为(xi,yi),i=1,...,k,然后根据如下公式计算每个边缘像素的补偿值:
其中,i为选定的边缘像素索引号,|gx|和|gy|表示x方向和y方向梯度值,||g(x,y)||为像素(x,y)的梯度幅值;将式(6)得到补偿值与边缘像素坐标相加得到亚像素边缘坐标(xi+ηx,yi+ηy),这些点是目标光斑边缘区域梯度变化最大的位置;
步骤四、对步骤三中得到的边缘点椭圆拟合,然后选取[0.9rmax,1.1rmax]范围内的像素组成待选像素集合{PixelTBS};
步骤五、对目标光斑兴趣区域的灰度值进行灰度投票统计,得到灰度值统计数列{GrayStatic};
步骤六、对步骤五中得到数列{GrayStatic}求导,并根据如下约束选取最大统计值对应的灰度值Imax1和第二大统计值对应的灰度值Imax2:
|Imax1-Imax2|>150 (7)
步骤七、在步骤六筛选得到的Imax1和Imax2附近找统计数量百分比大于1%的灰度值,分别与Imax1、Imax2组成背景阈值数列{Th_Background_GrayValues}和中心区域阈值数列{Th_Center_GrayValues};采用如下计算公式得到背景阈值和中心区域阈值:
步骤八、对步骤四中得到的待选像素集合{PixelTBS}使用步骤七得到的两个阈值进行筛选,所有在阈值范围内的像素(Th_Ibg<Im<Th_Icen)组成集合{Valid_Pixel},参与后续中心坐标计算,而不在阈值范围内的像素灰度值置零;
步骤九、利用步骤八中有效像素{Valid_Pixel}与步骤三中使用公式(6)计算得到的边缘点在梯度方向的距离计算每个有效像素的距离权重因子Wdi(i=1,...,k);距离权重因子的计算公式为:
其中,dist为像素i到亚像素边缘的距离;
步骤十、采用二维离散高斯函数对步骤八得到的有效像素集合{Valid_Pixel}拟合,根据得到的拟合参数计算每个有效像素的灰度权重因子Wgi(i=1,...,k);采用的二维高斯拟合函数为:
使用式(11)作为目标函数进行非线性优化,求得Wgi(i=1,...,k)最优解;具体计算公式为:
其中,Ii为像素i=1,...,k的灰度值;
步骤十一、将像素灰度值Ii和步骤十计算的灰度权重因子Wgi进行复合,得到新的影响因子λi,复合公式为:
步骤十二、将步骤十一计算得出的影响因子代入式(13)计算目标光斑中心位置(xc,yc):
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210128326.7A CN114565565A (zh) | 2022-02-11 | 2022-02-11 | 一种视觉测量目标中心亚像素定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210128326.7A CN114565565A (zh) | 2022-02-11 | 2022-02-11 | 一种视觉测量目标中心亚像素定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114565565A true CN114565565A (zh) | 2022-05-31 |
Family
ID=81712997
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210128326.7A Pending CN114565565A (zh) | 2022-02-11 | 2022-02-11 | 一种视觉测量目标中心亚像素定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114565565A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115205317A (zh) * | 2022-09-15 | 2022-10-18 | 山东高速集团有限公司创新研究院 | 一种桥梁监测光电标靶图像光斑中心点提取方法 |
CN115393441A (zh) * | 2022-10-27 | 2022-11-25 | 长春理工大学 | 一种基于液晶光闭环系统的光斑高精度提取定位方法 |
-
2022
- 2022-02-11 CN CN202210128326.7A patent/CN114565565A/zh active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115205317A (zh) * | 2022-09-15 | 2022-10-18 | 山东高速集团有限公司创新研究院 | 一种桥梁监测光电标靶图像光斑中心点提取方法 |
CN115393441A (zh) * | 2022-10-27 | 2022-11-25 | 长春理工大学 | 一种基于液晶光闭环系统的光斑高精度提取定位方法 |
CN115393441B (zh) * | 2022-10-27 | 2023-01-20 | 长春理工大学 | 一种基于液晶光闭环系统的光斑高精度提取定位方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112797915B (zh) | 一种线结构光测量系统的标定方法、标定装置、以及系统 | |
CN104331896B (zh) | 一种基于深度信息的系统标定方法 | |
CN114565565A (zh) | 一种视觉测量目标中心亚像素定位方法 | |
CN109523595B (zh) | 一种建筑工程直线棱角间距视觉测量方法 | |
CN111563878B (zh) | 一种空间目标定位方法 | |
CN106600648A (zh) | 标定相机内参数和畸变系数的立体编码靶标及其标定方法 | |
CN109976391B (zh) | 光学聚焦镜头的焦点感光元件三维空间视觉伺服定位系统及方法 | |
CN108168541B (zh) | 一种改进的亚像元星点质心定位方法 | |
CN110889829A (zh) | 一种基于鱼眼镜头的单目测距方法 | |
CN113012234B (zh) | 基于平面变换的高精度相机标定方法 | |
CN107589069B (zh) | 一种物体碰撞恢复系数的非接触式测量方法 | |
CN110261069B (zh) | 一种用于光学镜头的检测方法 | |
CN111508011A (zh) | 一种飞行时间相机的深度数据校准方法 | |
CN114018932A (zh) | 基于矩形标定物的路面病害指标测量方法 | |
US20050206874A1 (en) | Apparatus and method for determining the range of remote point light sources | |
CN114359405A (zh) | 一种离轴沙姆3d线激光相机的标定方法 | |
CN109974618A (zh) | 多传感器视觉测量系统的全局标定方法 | |
CN115201883A (zh) | 一种运动目标视频定位测速系统及方法 | |
CN116625258A (zh) | 链条间距测量系统及链条间距测量方法 | |
CN109506629B (zh) | 一种水下核燃料组件检测装置旋转中心标定的方法 | |
CN113008158B (zh) | 多线激光轮胎花纹深度测量方法 | |
CN114199160A (zh) | 基于二值编码光栅散焦投影的电路板元器件几何检测方法 | |
CN112132891A (zh) | 一种扩大标定空间的方法 | |
CN110440792A (zh) | 基于小天体不规则度的导航信息提取方法 | |
CN109242911A (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 |