连铸坯三维测量多线结构光的中心条纹快速提取

连铸坯三维测量多线结构光的中心条纹快速提取 Fast extraction center stripe of multi-line-structured light for three-dimensional measurement of casting slab 1 first-author 孙传富（1999-），男，安徽阜阳人，硕士研究生，主要从事机器视觉、光学三维测量方面的研究。E-mail： scf@shu.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39926830&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39926840&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39926835&type= 22.01333237 32.00399780 孙传富 （1999-），男，安徽阜阳人，硕士研究生，主要从事机器视觉、光学三维测量方面的研究。E-mail： scf@shu.edu.cn scf@shu.edu.cn 1 2 2 1 corresp E-mail： zjzhang@staff.shu.edu.cn E-mail： zjzhang@staff.shu.edu.cn 张之江（1968-），男，上海人，博士，教授，博士生导师，主要从事计算机视觉、图像处理方面的研究。E-mail：zjzhang@staff.shu.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39926845&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39926855&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39926849&type= 22.01332855 32.00400543 张之江 （1968-），男，上海人，博士，教授，博士生导师，主要从事计算机视觉、图像处理方面的研究。E-mail： zjzhang@staff.shu.edu.cn zjzhang@staff.shu.edu.cn 上海大学 特种光纤与光接入网重点实验室，上海 200444 上海大学 特种光纤与光接入网重点实验室，上海 200444 Key Laboratory of Specialty Fiber Optics and Optical Access Networks， Shanghai University， Shanghai 200444， China Key Laboratory of Specialty Fiber Optics and Optical Access Networks， Shanghai University， Shanghai 200444， China 中冶宝钢技术服务有限公司，上海 201999 中冶宝钢技术服务有限公司，上海 201999 Metallurgical Baosteel Technology Service Co. LTD.， Shanghai 201999， China Metallurgical Baosteel Technology Service Co. LTD.， Shanghai 201999， China 在利用多线结构光系统对大尺度连铸坯进行三维测量的过程中，线结构光中心条纹的提取精度是影响三维测量精度的重要因素。提出了一种对传统灰度质心方法进行优化的多线结构光中心条纹提取方法。对连铸坯表面激光横截面不满足高斯分布的线结构光条纹，利用图像背景差分的方式去除环境噪声，并根据线结构光条纹与背景间的灰度变化信息确定线结构光的边界，同时提取线结构光的感兴趣区域。根据光条在梯度方向上的灰度积分比例计算灰度质心法的自适应灰度阈值，利用灰度质心方法处理感兴趣区域提取出光条中心点，最后结合附近半径为5个像素的邻域内所定位的中心点进行质心重提取，获得连铸坯表面光条的亚像素中心坐标。现场测量结果表明：对反射特性不均匀的连铸坯表面的激光条纹光条中心进行提取，最终连铸坯边缘轨迹点的三维测量结果标准偏差在2 mm以内，该方法具有精度高、速度快、鲁棒性强的特点。 In the three-dimensional （3D） measurement of large-scale casting slabs using a multi-line structured light system， the accuracy of extracting the central stripe of the structured light is an important factor affecting the final 3D measurement accuracy. In this study， a multi-line structured light center stripe extraction method based on the optimization of the traditional gray gravity method is proposed. For linear structured light stripes projected onto the slab surface whose laser cross-section does not fit a Gaussian distribution， firstly， the environmental noise is removed by analyzing the image background difference. The boundary of linear structured light is determined according to the change in grayscale value between the linear structured light stripe and the background， and the region of interest of structured light is extracted. Then， the adaptive gray threshold of the gray gravity method is calculated according to the gray integral proportion of the light stripe in the gradient direction. The gray gravity method is used to process the region of interest and extract the approximate center point of the light stripe. Finally， the sub-pixel center coordinates of the light stripe are obtained by optimizing the extraction combined with the gravity of the roughly located center point in the neighborhood with a radius of 5 pixels. The experimental results show that this method has high precision， high speed， and strong robustness in extracting the center of the laser stripe on the surface of a casting slab with uneven reflection characteristics. The standard deviation of the 3D measurement results of the final casting slab edge trajectory points is within 2 mm. 三维测量 多线结构光 质心法 光条中心提取 连铸坯 three-dimensional measurement multi-line-structured light gravity method stripe center extraction casting slab 1 引 言 1 在连铸坯的生产过程中，由钢水凝固产生的连铸坯受到高温的影响，通常会产生许多表面质量问题，不利于后续的生产加工与处理，因此，需要对连铸坯边缘进行倒角来提高连铸坯的质量。连铸坯的倒角工序是在连铸坯的边缘位置利用火焰吹扫枪切削边缘，需要一套高精度的测量系统来测量连铸坯边缘所在的位置，进而配合执行机构上的火焰清理完成连铸坯的倒角工序。 近年来，结构光三维测量系统因测量速度快、精度高而得到了广泛应用 ［ 1 1 ］ 。结构光测量方法主要分为面结构光测量与线结构光测量。其中，面结构光测量是利用投影仪将预先编码的图案投影到被测物体上，并通过相机记录包含三维形状信息的畸变图案 ［ 2 2 ］ ，然后利用相移算法等方式提取相位图以重建对象的轮廓 ［ 3 3 ］ 。线结构光测量系统主要由线结构光发生器以及工业相机组成。基于光学三角法的原理，由线结构光发生器向待测物体表面投射线结构光，线结构光条纹受到待测物体表面轮廓的调制会产生形变。再由工业相机采集线结构光图像并分析线结构光条纹的形变信息，进而获得待测物体的三维信息 ［ 4 4 ］ 。 在连铸坯测量现场中，被测连铸坯尺度较大，单个结构光发生器无法实现高精度测量；且测量环境中存在高温、粉尘等影响，若选用功耗较大的多个投影设备进行现场测量，难以保证系统的稳定性且成本较高，因此，通常选用功率较低的多个线结构光发生器进行连铸坯的三维测量。 在线结构光视觉测量中，线结构光条纹中心的提取尤为重要，直接影响最终的测量精度 ［ 5 5 ］ 。传统的线结构光中心条纹提取方法根据提取精度的不同主要分为像素级和亚像素级两大类。其中，像素级提取方法有极值法 ［ 6 6 ］ 、阈值法 ［ 7 7 ］ 以及骨架细化法 ［ 8 8 ］ 等。这些方法获得的光条中心精度较低，在高精度视觉测量中应用较少。亚像素级方法包括灰度质心法 ［ 9 9 ］ 、Steger算法 ［ 10 10 - 11 11 10-11 ］ 等。其中，灰度质心法是利用图像中光条分布区域内的灰度信息计算质心并将它作为线结构光的中心，该方法能够减少条纹灰度分布不对称引起的误差且计算速度较快。Steger算法是利用Hessian矩阵计算线结构光条纹的法线方向，将法线方向上的极值点作为最终的线结构光中心。该方法精度高、鲁棒性好，但其计算量较大，速度相对较慢 ［ 12 12 ］ 。 传统方法在线结构光中心提取上的速度及精度均有不足。为了提高提取速度及精度，文献［ 13 13 ］在传统重心方法上结合轮廓跟踪算法对图像中的线结构光部分进行分析，显著降低了线结构光图像光条中心计算的复杂度。文献［ 14 14 ］提出一种自适应阈值的内部推进算法进行光条中心提取，在鲁棒性、速度、精度上均具有良好的效果。针对粗糙金属表面上不符合高斯分布的线结构光中心提取，文献［ 15 15 ］对粗糙金属表面的反射特性进行了分析，并根据反射特性设计了一系列的图像增强与分割方法对线结构光图像进行预处理，并利用灰度质心法进行光条中心的提取。该方法在图像预处理上具有良好的效果，但在中心点提取中直接使用传统灰度质心法精度略有不足。文献［ 16 16 ］根据金属表面的不均匀反射特性提出利用区域增长的方式消除激光散斑噪声的影响，并利用灰度质心法的结果作为初值，结合Sobel算子计算光条的梯度向量，在梯度方向利用双线性插值的方式得到光条中心的亚像素坐标。该方法计算过程较为复杂，在处理多线结构光条纹时用时较长，不满足快速测量的需求。 本文提出了一种多线结构光中心条纹提取方法，首先对连铸坯图像进行背景差分等处理以获取线结构光图像，利用线结构光边缘定位光条在图像中的位置并提取感兴趣区域（Region of Interest， ROI），以传统灰度质心方法获得光条中心点并对该中心点进行优化，最后得到高精度的线结构光条纹的亚像素中心坐标。此方法应用于连铸坯现场三维测量中取得了良好的效果。 2 图像分析与处理 1 2.1　连铸坯表面反射特性分析 2 连铸坯是由钢水直接浇铸产生的粗加工工件，其表面平整性较差，且存在颗粒较大的钢渣等残留物。当线结构光投射到连铸坯表面时，表面的不平整使得投射在其表面的激光条纹发生无规则散射。根据粗糙表面反射理论，激光以一定角度入射到粗糙物体表面时，入射光发生散射 ［ 15 15 ］ 。该部分的散射分为镜面反射以及无规则散射，前者遵循几何光学反射定律；后者则形成光斑状的散射中心，从中心向两端延伸的散射光带，对光条成像的影响较大 ［ 15 15 - 17 17 15-17 ］ 。 连铸坯三维测量系统在最初的结构设计阶段，考虑到现场测量环境以及相机与激光模组的安装位置条件，以及待测连铸坯与安装位置的距离，为保证线结构光投影在对应的相机视场内，投影线结构光以约45°的角度进行安装。连铸坯表面线结构光条纹的横截面并不能满足标准高斯分布，而是随入射角度发生偏移 ［ 18 18 ］ ，如 图1 图1 所示。除此之外，线结构光发生器与投射物体表面对焦不清晰，投射表面不平整性引起反射等情况，导致成像时激光灰度值变化曲线并不具有单峰性，会出现 图1 图1 （c）中的陡峭变化。 2.2　图像背景去噪 2 根据上述分析，相机成像平面上的光强噪声是由环境背景光、线结构光投影方向以及连铸坯表面的反射特性等多方面因素叠加产生的。工业现场测量环境中，光线变化较为常见，背景波动较为明显，所以图像噪声较大。连铸坯三维测量系统中，采集连铸坯的原始图像与线结构光投影图像之间的时间间隔较短，因此认为现场环境光对图像的影响是不变的 ［ 19 19 ］ 。连铸坯三维测量系统采集的连铸坯部分图像如 图2 图2 所示。 线结构光中心提取过程中，图像的有效信息为图像中线结构光条纹区域的灰度分布，而现场测量中环境光是以噪声的形式存在。如 图3 图3 所示， 图3 图3 （a）是不带有线结构光条纹的连铸坯灰度曲面，图像中明显存在大量的噪声。 图3 图3 （b）图为带有线结构光条纹的灰度曲面，将带有线结构光条纹的图像与不带有线结构光条纹的图像通过背景差分的方式可以有效去除环境噪声的影响 ［ 19 19 ］ 。 图3 图3 （c）为通过背景差分处理后的灰度曲面，可以明显看出在无线结构光条纹分布的区域，图像噪声明显降低。 背景差分法有效抑制了图像中的环境噪声，如 图4 图4 所示。 图4 图4 （a）为背景差分之后的图像，可以看出图像中环境背景噪声较大程度上受到抑制，线结构光条纹部分的特征保留。 图4 图4 （b）为 图4 图4 （a）中第4条激光线的片段图像。 图4 图4 （c）和 图4 图4 （d）分别对应 图4 图4 （b）在两个光条截面处的灰度统计图，可以看出，由于连铸坯表面的不平整性导致结构光无规则散射，光条投影效果依旧较差。需要进一步根据图像特征分离有效的ROI区域，并通过光条中心提取算法获取准确的光条中心点。 2.3　光条区域定位 2 在光条中心提取过程中，为避免随机噪声的影响同时加快处理速度，通常只对带有线结构光条纹的区域进行处理。如 图4 图4 （a）所示，图像包含多条激光线，需要首先区分图像中的各条激光线，再对单条激光线进行激光中心条纹提取，才能获得较为精确的激光中心条纹坐标。为获取多线结构光在图像中的分布信息，选取图像中单行的灰度数据进行分析。 图5 图5 为 图4 图4 （a）经背景差分后第2 560行的灰度统计图，可以看出，在图像中的8条线结构光条纹处图像灰度值变化较为明显，与背景灰度值间存在较大差异。可根据该特点进行结构光条纹的光条区域定位。 为确定单个激光条纹所在的ROI，需要对背景差分得到的图像进行处理，主要包括高斯模糊和自适应多阈值分割。选用5×5的高斯卷积核与原图进行卷积，再利用对图像亮度及对比度具有较强鲁棒性的OTSU算法对高斯模糊后的图像进行阈值分割 ［ 20 20 ］ 。阈值分割处理后的图像根据 <math id="M1"><mi>x</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39926933&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39926930&type= 2.20133328 2.79399991 方向上其像素值的变化确定激光线所在的大致区域。 根据激光条纹横截面灰度值呈现突变分布的性质，选取合适的光条纹灰度阈值以确定激光条纹在图像中的具体位置。在激光线的正常分布区域，取光条横截面上灰度值大于光条纹灰度阈值的两个临界像素点分别对应每条激光线的左右边缘。并根据激光条纹呈直线分布的特点，选取左右边缘特征点分别拟合直线，根据两条直线之间的位置关系便可依次有序分隔定位激光线。具体实现过程如下： （1） 选取合适的间隔值 <math id="M2"><mi>d</mi><mi>i</mi><mi>s</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39926942&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39926937&type= 4.48733330 2.79399991 沿图像的 <math id="M3"><mi>y</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39926950&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39926946&type= 2.11666679 3.72533321 方向遍历图像像素点，若该像素点在 <math id="M4"><mi>x</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927024&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927019&type= 2.20133328 2.79399991 方向上的前一个像素点灰度值 <math id="M5"><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>n</mi><mo>⋅</mo><mi mathvariant="normal">d</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">s</mi></mrow></mfenced></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927035&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927029&type= 12.61533356 3.89466691 小于激光线低阈值 <math id="M6"><mi>l</mi><mo>_</mo><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927274&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927271&type= 8.89000034 3.30200005 ，且下一个像素点灰度值 <math id="M7"><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>+</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>n</mi><mo>⋅</mo><mi mathvariant="normal">d</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">s</mi></mrow></mfenced></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927054&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927049&type= 12.61533356 3.89466691 大于激光线阈值 <math id="M8"><mi>h</mi><mo>_</mo><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927280&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927277&type= 9.65200043 3.30200005 ，则像素点 <math id="M9"><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>n</mi><mo>⋅</mo><mi mathvariant="normal">d</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">s</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927006&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927001&type= 14.30866718 3.97933316 判定为激光线左边缘点，表示为： <math display="block" id="M10"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>n</mi><mo>⋅</mo><mi mathvariant="normal">d</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">s</mi></mrow></mfenced></mrow></msub><mo><</mo><mi>l</mi><mo>_</mo><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi></mtd></mtr><mtr><mtd><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>+</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>n</mi><mo>⋅</mo><mi mathvariant="normal">d</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">s</mi></mrow></mfenced></mrow></msub><mo>></mo><mi>h</mi><mo>_</mo><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927015&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927011&type= 28.87133408 10.32933331 . 若该像素点在 <math id="M11"><mi>x</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927024&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927019&type= 2.20133328 2.79399991 方向上的上一个像素点的灰度值 <math id="M12"><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>n</mi><mo>⋅</mo><mi mathvariant="normal">d</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">s</mi></mrow></mfenced></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927035&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927029&type= 12.61533356 3.89466691 大于激光线阈值 <math id="M13"><mi>h</mi><mo>_</mo><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927280&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927277&type= 9.65200043 3.30200005 ，且下一个像素点灰度值 <math id="M14"><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>+</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>n</mi><mo>⋅</mo><mi mathvariant="normal">d</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">s</mi></mrow></mfenced></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927054&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927049&type= 12.61533356 3.89466691 小于激光线低阈值 <math id="M15"><mi>l</mi><mo>_</mo><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927274&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927271&type= 8.89000034 3.30200005 ，则判定像素点 <math id="M16"><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>n</mi><mo>⋅</mo><mi>d</mi><mi>i</mi><mi>s</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927070&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927065&type= 14.56266594 3.97933316 为激光线右边缘点，表示为： <math display="block" id="M17"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>n</mi><mo>⋅</mo><mi mathvariant="normal">d</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">s</mi></mrow></mfenced></mrow></msub><mo>></mo><mi>h</mi><mo>_</mo><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi></mtd></mtr><mtr><mtd><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>+</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>n</mi><mo>⋅</mo><mi mathvariant="normal">d</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">s</mi></mrow></mfenced></mrow></msub><mo><</mo><mi>l</mi><mo>_</mo><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927081&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927077&type= 28.87133408 10.32933331 . （2） 由于同一间隔行中连续两个左右边缘点定位一条线结构光条纹，则每一行可定位图像投影线结构光数量线结构光条纹。实际环境中，由于投影结构光在图像中较弱导致左右边缘所定位的结构光数量小于图像投影线的结构光数量，因此，选取定位线结构光条纹数量等于图像投影线结构光数量的左右边缘点并利用最小二乘法进行线性回归拟合直线。将定位激光线数量不等于图像投影线结构光数量的左右边缘点分别进行不同的直线匹配，计算边缘点到直线的欧式距离，以最小欧式距离的直线与边缘点进行匹配，完成线结构光边缘匹配，实现光条纹的完整定位。 （3） 单条激光线ROI的确定。根据上一步骤中所确定的单条激光线的左右边缘点，选取合适的左右边缘距离进行ROI提取与拼接，实现单条激光线ROI区域提取。 2.4　中心提取 2 在背景差分的图像中进行光条区域定位，将多线结构光图像分割成多个单条激光线的ROI区域。由于连铸坯表面的复杂反射特性，光条在不同位置处的宽度不均匀。为获取准确的光条中心亚像素坐标，本文在传统灰度质心法的基础上，在光条纹法线方向上，采用自适应灰度阈值结合光条中心点邻域信息对线结构光条纹中心进行优化，获得高精度的光条中心。步骤如下： （1） 利用传统灰度质心法获取光条行方向上的中心点； （2） 根据行方向中心点求解光条梯度方向上的自适应灰度阈值，重新利用该阈值计算光条中心点； （3） 利用质心方法对中心点进行重提取，获取高精度的光条亚像素中心点。 2 . 4 . 1　传统灰度质心法 3 光条中心提取时，首先利用灰度质心方法处理ROI提取出行方向的中心点。灰度质心法是对预处理之后图像线结构光截面上的多个像素及其灰度信息加权计算求得条纹的亚像素中心，计算公式如下 ［ 21 21 ］ ： <math display="block" id="M18"><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi></mrow></msub><mfenced separators="|"><mrow><msub><mrow><mi>x</mi></mrow><mrow><mi>p</mi></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>p</mi></mrow></msub></mrow></mfenced><mo>=</mo><mfenced separators="|"><mrow><mfrac><mrow><msub><mo>∑</mo><mrow><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></mfenced></mrow></msub><mo>></mo><msub><mrow><mi>t</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></msub><mrow><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></mfenced></mrow></msub><mo>⋅</mo><mi>i</mi></mrow></mrow><mrow><msub><mo>∑</mo><mrow><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></mfenced></mrow></msub><mo>></mo><msub><mrow><mi>t</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></msub><mrow><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></mfenced></mrow></msub></mrow></mrow></mfrac><mo>,</mo><mi>j</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927091&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927087&type= 47.83666611 13.97000027 ， 其中： <math id="M19"><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi></mrow></msub><mfenced separators="|"><mrow><msub><mrow><mi>x</mi></mrow><mrow><mi>p</mi></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>p</mi></mrow></msub></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927100&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927096&type= 15.49400043 4.82600021 是亚像素中心坐标， <math id="M20"><mi>I</mi><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927109&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927104&type= 9.48266602 4.74133301 是图像在第 <math id="M21"><mi>i</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927118&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927114&type= 1.10066664 2.79399991 行第 <math id="M22"><mi>j</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927129&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927123&type= 1.60866666 3.64066648 列的像素灰度值， <math id="M23"><msub><mrow><mi>t</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927137&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927132&type= 2.28600001 3.72533321 表示光条纹阈值。 2 . 4 . 2　自适应灰度阈值 3 在传统灰度质心方法中，光条纹阈值一定程度上影响了光条中心的计算结果。连铸坯测量中，当现场环境光较强或是线结构光投影较弱时，背景差分方式获得的图像受到现场环境光照的影响，光条纹区域的灰度值较小，使用固定阈值的灰度质心法会导致部分光条中心点的缺失。因此，本文利用自适应灰度质心法的阈值来解决此类问题。 以行方向上的中心点为中心，利用Sobel算子计算一定范围像素点的梯度方向 ［ 22 22 ］ ，并统计梯度方向相同的像素点数量，以数量最多的梯度方向作为光条在该区域内的梯度方向。沿着梯度方向上的像素点作为光条截面的像素点，通过光条截面像素点的灰度积分比例获取自适应灰度阈值。灰度光条截面上的像素点自适应灰度阈值为： <math display="block" id="M24"><mi mathvariant="normal">e</mi><mi mathvariant="normal">m</mi><mi mathvariant="normal">p</mi><mo>=</mo><mfrac><mrow><msub><mo>∑</mo><mrow><msub><mrow><mi>I</mi></mrow><mrow><mo stretchy="false">(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo stretchy="false">)</mo></mrow></msub><mo>></mo><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow></msub><mrow><mfenced separators="|"><mrow><msub><mrow><mi>I</mi></mrow><mrow><mo stretchy="false">(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo stretchy="false">)</mo></mrow></msub><mo>-</mo><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow></msub></mrow></mfenced></mrow></mrow><mrow><mo>∑</mo><msub><mrow><mi>I</mi></mrow><mrow><mo stretchy="false">(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo stretchy="false">)</mo></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927147&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927142&type= 40.38600159 11.93799973 ， 式中： <math id="M25"><mi mathvariant="normal">e</mi><mi mathvariant="normal">m</mi><mi mathvariant="normal">p</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927156&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927151&type= 6.51933336 3.64066648 表示灰度积分比例，为经验值； <math id="M26"><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></mfenced></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927165&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927160&type= 5.24933338 4.23333359 表示光条截面像素点的灰度值； <math id="M27"><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927176&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927172&type= 3.13266683 3.72533321 表示自适应灰度阈值。 在获得光条截面方向上的自适应灰度阈值后，再通过灰度质心法计算光条的中心点，公式如下： <math display="block" id="M28"><msub><mrow><mover accent="true"><mi>C</mi><mo>¯</mo></mover></mrow><mrow><mi>p</mi></mrow></msub><mfenced separators="|"><mrow><msub><mrow><mover accent="true"><mi>x</mi><mo>¯</mo></mover></mrow><mrow><mi>p</mi></mrow></msub><mo>,</mo><msub><mrow><mover accent="true"><mi>y</mi><mo>¯</mo></mover></mrow><mrow><mi>p</mi></mrow></msub></mrow></mfenced><mo>=</mo><mfenced separators="|"><mrow><mfrac><mrow><msub><mo>∑</mo><mrow><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></mfenced></mrow></msub><mo>></mo><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow></msub></mrow></msub><mrow><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></mfenced></mrow></msub><mo>⋅</mo><mi>i</mi></mrow></mrow><mrow><msub><mo>∑</mo><mrow><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></mfenced></mrow></msub><mo>></mo><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow></msub></mrow></msub><mrow><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></mfenced></mrow></msub></mrow></mrow></mfrac><mo>,</mo><mfrac><mrow><msub><mo>∑</mo><mrow><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></mfenced></mrow></msub><mo>></mo><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow></msub></mrow></msub><mrow><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></mfenced></mrow></msub><mo>⋅</mo><mi>j</mi></mrow></mrow><mrow><msub><mo>∑</mo><mrow><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></mfenced></mrow></msub><mo>></mo><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow></msub></mrow></msub><mrow><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></mfenced></mrow></msub></mrow></mrow></mfrac></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927184&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927179&type= 69.34200287 13.97000027 ， 式中： <math id="M29"><msub><mrow><mover accent="true"><mi>C</mi><mo>¯</mo></mover></mrow><mrow><mi>p</mi></mrow></msub><mfenced separators="|"><mrow><msub><mrow><mover accent="true"><mi>x</mi><mo>¯</mo></mover></mrow><mrow><mi>p</mi></mrow></msub><mo>,</mo><msub><mrow><mover accent="true"><mi>y</mi><mo>¯</mo></mover></mrow><mrow><mi>p</mi></mrow></msub></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927191&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927187&type= 15.49400043 4.82600021 表示光条中心点。 2 . 4 . 3　光条中心重提取 3 为提高线结构光条纹中心点精度，进一步优化获得较为平滑的中心条纹点，对由灰度质心法得到的线结构光条纹中心点进行优化。 如 图6 图6 所示，在图像坐标系 <math id="M30"><mi>o</mi><mo>-</mo><mi>x</mi><mi>y</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927199&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927196&type= 8.04333401 3.72533321 内由灰度质心法求出各点的初始结果， <math id="M31"><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi></mrow></msub><mfenced separators="|"><mrow><msub><mrow><mi>x</mi></mrow><mrow><mi>p</mi></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>p</mi></mrow></msub></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927204&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927201&type= 15.49400043 4.82600021 为其中一点。以 <math id="M32"><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927211&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927208&type= 3.72533321 3.72533321 为中心， <math id="M33"><msub><mrow><mi>r</mi></mrow><mrow><mi>p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927315&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927313&type= 2.62466669 3.72533321 为距离范围，选取范围内共 <math id="M34"><mfenced separators="|"><mrow><mn mathvariant="normal">2</mn><msub><mrow><mi>r</mi></mrow><mrow><mi>p</mi></mrow></msub><mo>+</mo><mn mathvariant="normal">1</mn></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927221&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927218&type= 13.97000027 4.82600021 个点重新计算质心 <math id="M35"><msub><mrow><mover accent="true"><mi>C</mi><mo>^</mo></mover></mrow><mrow><mi>p</mi></mrow></msub><mfenced separators="|"><mrow><msub><mrow><mover accent="true"><mi>x</mi><mo>^</mo></mover></mrow><mrow><mi>p</mi></mrow></msub><mo>,</mo><msub><mrow><mover accent="true"><mi>y</mi><mo>^</mo></mover></mrow><mrow><mi>p</mi></mrow></msub></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927226&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927223&type= 15.49400043 4.99533367 作为优化后的结果值，即有： <math display="block" id="M36"><msub><mrow><mover accent="true"><mi>C</mi><mo>^</mo></mover></mrow><mrow><mi>p</mi></mrow></msub><mfenced separators="|"><mrow><msub><mrow><mover accent="true"><mi>x</mi><mo>^</mo></mover></mrow><mrow><mi>p</mi></mrow></msub><mo>,</mo><msub><mrow><mover accent="true"><mi>y</mi><mo>^</mo></mover></mrow><mrow><mi>p</mi></mrow></msub></mrow></mfenced><mo>=</mo><mfenced separators="|"><mrow><mfrac><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>x</mi><mo>=</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>p</mi></mrow></msub><mo>-</mo><msub><mrow><mi>r</mi></mrow><mrow><mi>p</mi></mrow></msub></mrow><mrow><msub><mrow><mi>x</mi></mrow><mrow><mi>p</mi></mrow></msub><mo>+</mo><msub><mrow><mi>r</mi></mrow><mrow><mi>p</mi></mrow></msub></mrow></munderover></mstyle><mrow><mi>x</mi><mo>⋅</mo><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced></mrow></msub></mrow></mrow><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>x</mi><mo>=</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>p</mi></mrow></msub><mo>-</mo><msub><mrow><mi>r</mi></mrow><mrow><mi>p</mi></mrow></msub></mrow><mrow><msub><mrow><mi>x</mi></mrow><mrow><mi>p</mi></mrow></msub><mo>+</mo><msub><mrow><mi>r</mi></mrow><mrow><mi>p</mi></mrow></msub></mrow></munderover></mstyle><mrow><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced></mrow></msub></mrow></mrow></mfrac><mo>,</mo><mfrac><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>x</mi><mo>=</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>p</mi></mrow></msub><mo>-</mo><msub><mrow><mi>r</mi></mrow><mrow><mi>p</mi></mrow></msub></mrow><mrow><msub><mrow><mi>x</mi></mrow><mrow><mi>p</mi></mrow></msub><mo>+</mo><msub><mrow><mi>r</mi></mrow><mrow><mi>p</mi></mrow></msub></mrow></munderover></mstyle><mrow><mi>y</mi><mo>⋅</mo><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced></mrow></msub></mrow></mrow><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>x</mi><mo>=</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>p</mi></mrow></msub><mo>-</mo><msub><mrow><mi>r</mi></mrow><mrow><mi>p</mi></mrow></msub></mrow><mrow><msub><mrow><mi>x</mi></mrow><mrow><mi>p</mi></mrow></msub><mo>+</mo><msub><mrow><mi>r</mi></mrow><mrow><mi>p</mi></mrow></msub></mrow></munderover></mstyle><mrow><msub><mrow><mi>I</mi></mrow><mrow><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced></mrow></msub></mrow></mrow></mfrac></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927230&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927228&type= 69.25733185 22.18266678 . 最终得到准确的线结构光中心点亚像素坐标 <math id="M37"><msub><mrow><mover accent="true"><mi>C</mi><mo>^</mo></mover></mrow><mrow><mi>p</mi></mrow></msub><mfenced separators="|"><mrow><msub><mrow><mover accent="true"><mi>x</mi><mo>^</mo></mover></mrow><mrow><mi>p</mi></mrow></msub><mo>,</mo><msub><mrow><mover accent="true"><mi>y</mi><mo>^</mo></mover></mrow><mrow><mi>p</mi></mrow></msub></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927249&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927246&type= 15.49400043 4.99533367 。 3 实 验 1 3.1　实验环境 2 实验环境为中冶宝钢连铸坯清理吹扫工作现场，连铸坯三维测量系统实物如 图7 图7 所示。系统中，工业相机使用大恒水星ME2P-2621-15U3M相机，相机分辨率为5 120×5 120 pixel，像素物理尺寸为 <math id="M38"><mn mathvariant="normal">2.5</mn><mtext> </mtext><mi mathvariant="normal">μ</mi><mi mathvariant="normal">m</mi><mo>×</mo><mn mathvariant="normal">2.5</mn><mtext> </mtext><mi mathvariant="normal">μ</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927255&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927252&type= 27.77066803 3.64066648 ；相机所装配的镜头为LM12FC24M型号的Kowa镜头，焦距为12 mm。激光器是由深圳森泰达公司生产的KYL650N200-X1670单线激光器，其功率为50 mW，波长为650 nm，激光光源距离待测连铸坯约1.7 m。 3.2　图像ROI提取与自适应阈值 2 为了测试本文所提出的ROI提取方法的有效性，实验中使用 图4 图4 （a）背景差分图像对连铸坯表面的线结构光进行ROI提取。根据大量现场采集图像得到的经验性数据，间隔值 <math id="M39"><mi mathvariant="normal">d</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">s</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927268&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927266&type= 4.31799984 2.79399991 取10，激光线低阈值 <math id="M40"><mi>l</mi><mo>_</mo><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927274&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927271&type= 8.89000034 3.30200005 取10，高阈值 <math id="M41"><mi>h</mi><mo>_</mo><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927280&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927277&type= 9.65200043 3.30200005 取30。经过图像背景差分后，根据该阈值对亮条纹以及暗条纹均可以快速定位出连铸坯表面线结构光条纹的位置。最终提取效果如 图8 图8 所示，对差分图像中的8条激光线自左到右依次提取片段。所获得的线结构光ROI片段均能完整地定位出线结构光在图像中的准确位置。 在单线结构光条纹ROI内进行自适应阈值计算，选取灰度积分比例 <math id="M42"><mi mathvariant="normal">e</mi><mi mathvariant="normal">m</mi><mi mathvariant="normal">p</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927290&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927288&type= 6.51933336 3.64066648 为0.7。设置自适应阈值的初值为0，并设定一个较小的步长进行迭代求解，当灰度积分比例大于等于经验值时，则可求得自适应阈值。所得到的自适应灰度阈值效果如 图9 图9 所示。从 图9 图9 可以明显观察到，自适应阈值可以准确地根据图像像素灰度值分离光条与图像背景，抑制背景中散斑噪声的影响。 3.3　光条中心线提取鲁棒性测试 2 为了证明本文所提出的光条提取方法的鲁棒性，采用不同方法与本文方法进行比较，得到的误差对比如 图10 图10 所示。其中，Steger方法、GGM方法以及本文提出的方法得到的是亚像素级别的光条中心坐标，极值法与几何中心法得到的是像素级别的中心点坐标。如 表1 表1 所示，极值法与几何中心法的像素中心点标准差明显高于其他3种亚像素中心提取方法。亚像素提取方法中，传统灰度质心法提取出的光条中心点标准差高于本文所提出的方法与Steger方法。本文提出的方法以及Steger方法处理连铸坯表面激光条纹得到的光条纹中心点的平均像素值与标准差较为接近，均可获得平滑性良好的光条中心点。 Method Mean pixel value Standard deviation Steger Method ［ 11 11 ］ 60.334 0.123 Extreme value ［ 6 6 ］ 60.235 0.970 GGM ［ 9 9 ］ 60.341 0.365 Proposed method 60.345 0.112 Geometric center ［ 16 16 ］ 60.000 0.707 pixel 进一步将本文所提出的方法应用于连铸坯三维测量系统所采集的其他图像中，无论在光条光强较弱还是光条过曝的情况下，本文所提出的方法均可以提取到精度高且较为平滑的条纹中心点，如 图11 图11 所示，进一步验证了本文方法的鲁棒性。 3.4　光条中心线提取速度测试 2 在速度测试实验中，用于图像处理的计算机的CPU型号为Intel i5-1135G7 2.40 GHz，RAM为16 G。软件环境为Visual Studio 2015。对实验使用现场所采集到的7组不同连铸坯图像进行测试，其中图像的尺寸为5 120×5 120 pixel。对于每组图像，使用传统灰度质心法、Steger方法以及本文提出的方法提取图像中完整的8个光条的中心条纹。其中，本文方法的 <math id="M43"><msub><mrow><mi>r</mi></mrow><mrow><mi>p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927315&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927313&type= 2.62466669 3.72533321 值取5。 在程序计算的复杂度上，传统灰度质心法的计算过程最为简单，本文方法是基于灰度质心法结果的优化，因此运行时间比传统灰度质心法更长；而Steger方法在计算过程中涉及到高斯核与图像中每个像素的卷积过程，因此运行速度最慢。所测试的数据在不同方法下光条中心点提取算法的运行时长如 表2 表2 所示。其中，GGM方法与本文方法利用8个线程对分离出的单个光条图像进行同步处理。从 表2 表2 可以看出，实际算法运行过程中灰度质心法算法的运行时长明显低于其他两种方法，Steger方法速度最慢，本文提出的方法居中。本文方法的单组图像平均处理时长约为300 ms，运行时长为传统灰度质心法的3倍左右，显著低于Steger方法，满足连铸坯三维测量系统的应用需求。 No. GGM ［ 9 9 ］ Steger method ［ 11 11 ］ Proposed method 1 103 5 136 346 2 95 5 139 292 3 89 5 157 284 4 112 5 655 318 5 100 5 136 288 6 99 5 211 288 7 100 5 080 286 ms 3.5　现场应用 2 为验证本文提出方法在工业测量现场应用中的可靠性，采用多线结构光中心提取算法结合搭建的连铸坯三维测量系统对现场长约11 m的连铸坯进行了三维测量。 连铸坯三维测量系统中共有7台线性排列的工业相机，除1，2号相机外，相邻两台工业相机之间视场相切。7台相机共采集视场内54条线结构光的图像，对投射到连铸坯表面的线结构光进行三维重建获取连铸坯上下边缘位置的三维点坐标，进而配合现场的控制火焰吹扫枪执行机构完成连铸坯的倒角工序。 连铸坯三维测量系统中，相机标定方法采用单目标定 ［ 23 23 ］ 算法得到相机内参，再通过公共视场的双目标定方法，得到多个相机间的外参矩阵。完成相机坐标系间的坐标转换，使其可统一到同一坐标系下。其中，相机单目标定的重投影误差约为0.005 pixel，双目标定的重投影误差约为0.1 pixel。光平面标定方法是利用交比不变法获取激光器发射的光平面在相机坐标系下的三维平面方程 ［ 24 24 ］ ，光平面的标定误差约为0.003 mm。 在获得图像上的光斑中心的二维像素坐标后，可通过相机内参获得在相机坐标系下的一条射线，与光平面相交，联立方程求解得到唯一点，即为所求的三维点坐标。 图12 图12 （a）为7台工业相机采集到的完整连铸坯图像。以 图12 图12 （a）中4号相机为例，部分光条中心提取效果如 图12 图12 （b）所示。将7台工业相机视场内的线结构光三维测量结果统一至1号相机坐标系下，形成的三维点云如 图12 图12 （c）所示。 1;2;3; 表3 表3 为连铸坯三维测量系统测得的 图12 图12 （a）中连铸坯表面线结构光上下边缘在上下两个执行机构坐标系下的部分三维坐标。根据最终测量结果中标准偏差进行系统性能评估，对上下边缘三维坐标数据中 x 坐标进行线性回归拟合直线，并计算标准偏差。最终测得标准偏差在2 mm以内，满足工业测量现场连铸坯三维坐标的测量要求。 2 No. 3 Upper edge points 3 Lower edge points x y z x y z 1 -491.200 238.83 1 619.94 -543.88 251.27 1 618.43 2 -495.130 575.36 1 633.69 -542.01 582.07 1 634.33 3 -496.750 747.65 1 636.70 -540.40 766.51 1 637.59 4 -497.510 966.847 1 641.35 -538.40 978.24 1 641.18 5 -498.730 1 213.17 1 647.47 -536.99 1 225.99 1 646.80 6 -498.686 1 480.00 1 654.00 -533.46 1 454.73 1 651.55 7 -500.940 1 717.11 1 659.15 -532.02 1 716.72 1 657.85 8 -500.909 1 880.56 1 661.92 -530.42 1 878.53 1 660.78 … … … … … … … Standard deviation 3 1.52 3 1.73 mm 4 结 论 1 针对连铸坯表面不平整特性导致的投影线结构光不规则散射等问题，本文提出了一种对连铸坯表面图像处理及多线结构光中心条纹提取的方法。通过分析背景差分后图像的灰度特性，可以准确提取线结构光条纹在图像中的位置，并利用自适应阈值法获取完整的线结构光条纹的ROI。在光条中心提取中，该方法相较于传统的灰度质心法，获得的条纹亚像素坐标更为准确，且速度相较于Steger算法更快。本文方法在传统灰度质心法上进行改进，对连铸坯三维测量系统图像中心的提取效果较好，能够提取到亚像素光条纹中心点，最终测量结果的标准偏差在2 mm以内，适用于连铸坯的三维现场测量。 参考文献 LV S Z ， SUN Q ， ZHANG Y Y ， et al . Projector distortion correction in 3D shape measurement using a structured-light system by deep neural networks ［J］. Optics Letters ， 2020 ， 45 （ 1 ）： 204 - 207 . doi: 10.1364/ol.45.000204 http://dx.doi.org/10.1364/ol.45.000204 Projector distortion correction in 3D shape measurement using a structured-light system by deep neural networks HE X Y ， ZHENG D L ， QIAN K M ， et al . Quaternary gray-code phase unwrapping for binary fringe projection profilometry ［J］. Optics and Lasers in Engineering ， 2019 ， 121 ： 358 - 368 . doi: 10.1016/j.optlaseng.2019.04.009 http://dx.doi.org/10.1016/j.optlaseng.2019.04.009 Quaternary gray-code phase unwrapping for binary fringe projection profilometry LV S Z ， SUN Q ， ZHANG Y Y ， et al . Monotonicity analysis of absolute phase unwrapping by geometric constraint in a structured light system ［J］. Optics Express ， 2020 ， 28 （ 7 ）： 9885 . doi: 10.1364/oe.386646 http://dx.doi.org/10.1364/oe.386646 Monotonicity analysis of absolute phase unwrapping by geometric constraint in a structured light system 王利 ， 陈念年 ， 巫玲 ， 等 . 高噪声背景下激光条纹亚像素中心的提取 ［J］. 应用光学 ， 2016 ， 37 （ 2 ）： 321 - 326 . doi: 10.5768/jao201637.0207004 http://dx.doi.org/10.5768/jao201637.0207004 高噪声背景下激光条纹亚像素中心的提取 WANG L ， CHEN N N ， WU L ， et al . Extraction of laser stripe sub-pixel center in high-noise background ［J］. Journal of Applied Optics ， 2016 ， 37 （ 2 ）： 321 - 326 . （in Chinese） . doi: 10.5768/jao201637.0207004 http://dx.doi.org/10.5768/jao201637.0207004 Extraction of laser stripe sub-pixel center in high-noise background 张瑞瑛 ， 周萍 ， 冯煦 ， 等 . 大视场下线结构光光条中心的快速提取 ［J］. 应用光学 ， 2010 ， 31 （ 3 ）： 432 - 436 . doi: 10.3969/j.issn.1002-2082.2010.03.019 http://dx.doi.org/10.3969/j.issn.1002-2082.2010.03.019 大视场下线结构光光条中心的快速提取 ZHANG R Y ， ZHOU P ， FENG X ， et al . Rapid extraction of line-structured light stripe in large field of view ［J］. Journal of Applied Optics ， 2010 ， 31 （ 3 ）： 432 - 436 . （in Chinese） . doi: 10.3969/j.issn.1002-2082.2010.03.019 http://dx.doi.org/10.3969/j.issn.1002-2082.2010.03.019 Rapid extraction of line-structured light stripe in large field of view EPSTEIN A ， TESSLER N ， EINZIGER P D . Analytical extraction of the recombination zone location in organic light-emitting diodes from emission pattern extrema ［J］. Optics Letters ， 2010 ， 35 （ 20 ）： 3366 - 3368 . doi: 10.1364/ol.35.003366 http://dx.doi.org/10.1364/ol.35.003366 Analytical extraction of the recombination zone location in organic light-emitting diodes from emission pattern extrema SAKUMA K ， HASEGAWA S ， TAKAHASHI H ， et al . Holographic laser sweeper for in-process debris removal ［J］. Applied Physics B ， 2015 ， 119 （ 3 ）： 533 - 538 . doi: 10.1007/s00340-015-6011-6 http://dx.doi.org/10.1007/s00340-015-6011-6 Holographic laser sweeper for in-process debris removal 潘硕 ， 邵新杰 ， 李晓磊 ， 等 . 基于块匹配的激光条纹亚像素中心提取 ［J］. 激光与红外 ， 2021 ， 51 （ 6 ）： 814 - 819 . doi: 10.3969/j.issn.1001-5078.2021.06.020 http://dx.doi.org/10.3969/j.issn.1001-5078.2021.06.020 基于块匹配的激光条纹亚像素中心提取 PAN SH ， SHAO X J ， LI X L ， et al . Sub-pixel extraction of laser stripe center based on block-matching ［J］. Laser ＆ Infrared ， 2021 ， 51 （ 6 ）： 814 - 819 . （in Chinese） . doi: 10.3969/j.issn.1001-5078.2021.06.020 http://dx.doi.org/10.3969/j.issn.1001-5078.2021.06.020 Sub-pixel extraction of laser stripe center based on block-matching LI Y H ， ZHOU J B ， HUANG F S ， et al . Sub-pixel extraction of laser stripe center using an improved gray-gravity method ［J］. Sensors （Basel， Switzerland） ， 2017 ， 17 （ 4 ）： 814 . doi: 10.3390/s17040814 http://dx.doi.org/10.3390/s17040814 Sub-pixel extraction of laser stripe center using an improved gray-gravity method 南方 ， 李大华 ， 高强 ， 等 . 改进Steger算法的自适应光条纹中心提取 ［J］. 激光杂志 ， 2018 ， 39 （ 1 ）： 85 - 88 . doi: 10.14016/j.cnki.jgzz.2018.01.085 http://dx.doi.org/10.14016/j.cnki.jgzz.2018.01.085 改进Steger算法的自适应光条纹中心提取 NAN F ， LI D H ， GAO Q ， et al . Implementation of adaptive light stripe center extraction of improved Steger algorithm ［J］. Laser Journal ， 2018 ， 39 （ 1 ）： 85 - 88 . （in Chinese） . doi: 10.14016/j.cnki.jgzz.2018.01.085 http://dx.doi.org/10.14016/j.cnki.jgzz.2018.01.085 Implementation of adaptive light stripe center extraction of improved Steger algorithm STEGER C . An unbiased detector of curvilinear structures ［J］. IEEE Transactions on Pattern Analysis and Machine Intelligence ， 1998 ， 20 （ 2 ）： 113 - 125 . doi: 10.1109/34.659930 http://dx.doi.org/10.1109/34.659930 An unbiased detector of curvilinear structures 符钦伟 ， 杜宣 ， 张燕 ， 等 . 半透明胶条表面光条中心线提取 ［J］. 光学 精密工程 ， 2021 ， 29 （ 4 ）： 701 - 709 . doi: 10.37188/OPE.20212904.0701 http://dx.doi.org/10.37188/OPE.20212904.0701 半透明胶条表面光条中心线提取 FU Q W ， DU X ， ZHANG Y ， et al . Extraction of laser stripe center on translucent glue strip surface ［J］. Opt. Precision Eng. ， 2021 ， 29 （ 4 ）： 701 - 709 . （in Chinese） . doi: 10.37188/OPE.20212904.0701 http://dx.doi.org/10.37188/OPE.20212904.0701 Extraction of laser stripe center on translucent glue strip surface 李伟明 ， 彭国 ， 高兴宇 ， 等 . 线激光光条中心快速提取算法 ［J］. 中国激光 ， 2020 ， 47 （ 3 ）： 192 - 199 . doi: 10.3788/cjl202047.0304002 http://dx.doi.org/10.3788/cjl202047.0304002 线激光光条中心快速提取算法 LI W M ， PENG G ， GAO X Y ， et al . Fast extraction algorithm for line laser strip centers ［J］. Chinese Journal of Lasers ， 2020 ， 47 （ 3 ）： 192 - 199 . （in Chinese） . doi: 10.3788/cjl202047.0304002 http://dx.doi.org/10.3788/cjl202047.0304002 Fast extraction algorithm for line laser strip centers 李伟明 ， 梅枫 ， 胡增 ， 等 . 线激光光条中心内部推进提取算法 ［J］. 中国激光 ， 2021 ， 48 （ 11 ）： 105 - 117 . doi: 10.3788/cjl202148.1104002 http://dx.doi.org/10.3788/cjl202148.1104002 线激光光条中心内部推进提取算法 LI W M ， MEI F ， HU Z ， et al . Internal propulsion algorithm for extracting center of line laser stripe ［J］. Chinese Journal of Lasers ， 2021 ， 48 （ 11 ）： 105 - 117 . （in Chinese） . doi: 10.3788/cjl202148.1104002 http://dx.doi.org/10.3788/cjl202148.1104002 Internal propulsion algorithm for extracting center of line laser stripe 赵博华 ， 王伯雄 ， 张金 ， 等 . 粗糙金属表面光条中心提取方法 ［J］. 光学 精密工程 ， 2011 ， 19 （ 9 ）： 2138 - 2145 . doi: 10.3788/ope.20111909.2138 http://dx.doi.org/10.3788/ope.20111909.2138 粗糙金属表面光条中心提取方法 ZHAO B H ， WANG B X ， ZHANG J ， et al . Extraction of laser stripe center on rough metal surface ［J］. Opt. Precision Eng. ， 2011 ， 19 （ 9 ）： 2138 - 2145 . （in Chinese） . doi: 10.3788/ope.20111909.2138 http://dx.doi.org/10.3788/ope.20111909.2138 Extraction of laser stripe center on rough metal surface WANG H F ， WANG Y F ， ZHANG J J ， et al . Laser stripe center detection under the condition of uneven scattering metal surface for geometric measurement ［J］. IEEE Transactions on Instrumentation and Measurement ， 2020 ， 69 （ 5 ）： 2182 - 2192 . doi: 10.1109/tim.2019.2921440 http://dx.doi.org/10.1109/tim.2019.2921440 Laser stripe center detection under the condition of uneven scattering metal surface for geometric measurement 文晶娅 . 激光测距仪光学系统设计及目标漫反射特性研究 ［D］. 武汉 ： 华中科技大学 ， 2007 . 激光测距仪光学系统设计及目标漫反射特性研究 WEN J Y . Optical System Design of Laser Rangefinder and Research on Diffuse Reflectivity of Reflective Target ［D］. Wuhan ： Huazhong University of Science and Technology ， 2007 . （in Chinese） Optical System Design of Laser Rangefinder and Research on Diffuse Reflectivity of Reflective Target SUN Q C ， CHEN J ， LI C J . A robust method to extract a laser stripe centre based on grey level moment ［J］. Optics and Lasers in Engineering ， 2015 ， 67 ： 122 - 127 . doi: 10.1016/j.optlaseng.2014.11.007 http://dx.doi.org/10.1016/j.optlaseng.2014.11.007 A robust method to extract a laser stripe centre based on grey level moment TSAI D M ， LAI S C . Independent component analysis-based background subtraction for indoor surveillance ［J］. IEEE Transactions on Image Processing ， 2009 ， 18 （ 1 ）： 158 - 167 . doi: 10.1109/tip.2008.2007558 http://dx.doi.org/10.1109/tip.2008.2007558 Independent component analysis-based background subtraction for indoor surveillance GAO H P ， XU G L . Robust analysis and laser stripe center extraction for rail images ［J］. Applied Sciences ， 2021 ， 11 （ 5 ）. doi: 10.3390/app11052038 http://dx.doi.org/10.3390/app11052038 Robust analysis and laser stripe center extraction for rail images 贾振元 ， 樊超楠 ， 刘巍 ， 等 . 双目视觉测量中等匹配点的光条中心提取 ［J］. 光学 精密工程 ， 2016 ， 24 （ 7 ）： 1582 - 1591 . doi: 10.3788/ope.20162407.1582 http://dx.doi.org/10.3788/ope.20162407.1582 双目视觉测量中等匹配点的光条中心提取 JIA ZH Y ， FAN CH N ， LIU W ， et al . Extraction of laser stripe centers based on equal matching points for binocular vision measurement ［J］. Opt. Precision Eng. ， 2016 ， 24 （ 7 ）： 1582 - 1591 . （in Chinese） . doi: 10.3788/ope.20162407.1582 http://dx.doi.org/10.3788/ope.20162407.1582 Extraction of laser stripe centers based on equal matching points for binocular vision measurement LI Z K ， MA L P ， LONG X L ， et al . Hardware-oriented algorithm for high-speed laser centerline extraction based on hessian matrix ［J］. IEEE Transactions on Instrumentation and Measurement ， 2021 ， 70 ： 1 - 14 . doi: 10.1109/tim.2021.3081163 http://dx.doi.org/10.1109/tim.2021.3081163 Hardware-oriented algorithm for high-speed laser centerline extraction based on hessian matrix STROBL K H ， HIRZINGER G . More accurate pinhole camera calibration with imperfect planar target ［C］. 2011 IEEE International Conference on Computer Vision Workshops . Barcelona， Spain . IEEE ， 2011 ： 1068 - 1075 . doi: 10.1109/iccvw.2011.6130369 http://dx.doi.org/10.1109/iccvw.2011.6130369 More accurate pinhole camera calibration with imperfect planar target 许丽 ， 张之江 . 结构光测量系统的误差传递分析 ［J］. 光学 精密工程 ， 2009 ， 17 （ 2 ）： 306 - 313 . doi: 10.3321/j.issn:1004-924X.2009.02.012 http://dx.doi.org/10.3321/j.issn:1004-924X.2009.02.012 结构光测量系统的误差传递分析 XU L ， ZHANG ZH J . Error propagation analysis of structured light system ［J］. Opt. Precision Eng. ， 2009 ， 17 （ 2 ）： 306 - 313 . （in Chinese） . doi: 10.3321/j.issn:1004-924X.2009.02.012 http://dx.doi.org/10.3321/j.issn:1004-924X.2009.02.012 Error propagation analysis of structured light system 10.37188/OPE.20233103.0380 TP391 A 1004-924X（2023）03-0380-13 孙传富,彭涛,陆永刚等.连铸坯三维测量多线结构光的中心条纹快速提取[J].光学精密工程,2023,31(03):380-392. 孙传富,彭涛,陆永刚等.连铸坯三维测量多线结构光的中心条纹快速提取[J].光学精密工程,2023,31(03):380-392. SUN Chuanfu,PENG Tao,LU Yonggang,et al.Fast extraction center stripe of multi-line-structured light for three-dimensional measurement of casting slab[J].Optics and Precision Engineering,2023,31(03):380-392. SUN Chuanfu,PENG Tao,LU Yonggang,et al.Fast extraction center stripe of multi-line-structured light for three-dimensional measurement of casting slab[J].Optics and Precision Engineering,2023,31(03):380-392. 国家自然科学基金面上项目（No.61572307） 2022-08-03 2022-09-06 2023-02-10 2023-02-15 光学精密工程 3 31