[go: up one dir, main page]

CN115014248A - A method for identification and flatness determination of laser projection lines - Google Patents

A method for identification and flatness determination of laser projection lines Download PDF

Info

Publication number
CN115014248A
CN115014248A CN202210650078.2A CN202210650078A CN115014248A CN 115014248 A CN115014248 A CN 115014248A CN 202210650078 A CN202210650078 A CN 202210650078A CN 115014248 A CN115014248 A CN 115014248A
Authority
CN
China
Prior art keywords
laser
line
plane
image
pixel
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.)
Granted
Application number
CN202210650078.2A
Other languages
Chinese (zh)
Other versions
CN115014248B (en
Inventor
杜万和
刘明言
杨敬辉
韩玲玲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Polytechnic University
Original Assignee
Shanghai Polytechnic University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Polytechnic University filed Critical Shanghai Polytechnic University
Priority to CN202210650078.2A priority Critical patent/CN115014248B/en
Publication of CN115014248A publication Critical patent/CN115014248A/en
Application granted granted Critical
Publication of CN115014248B publication Critical patent/CN115014248B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/30Measuring arrangements characterised by the use of optical techniques for measuring roughness or irregularity of surfaces

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

The invention relates to the technical field of laser detection, in particular to a method for identifying a laser projection line and judging the flatness, which vertically shoots a laser line image through a camera, firstly extracts the color, the outline and the characteristics of the laser projection line, carries out filtering treatment on the laser projection line to reserve effective pixels of the laser main line, uses the effective pixels, then carries out plane fitting on all effective pixels by using a least square method, calculates the current plane discrete degree, simultaneously simulates a three-dimensional image when the levelness deviates by 20 seconds, calculates the plane discrete degree when the levelness deviates by 20 seconds, and contrasts and judges whether the flatness of the laser projection line is qualified or not; the method can detect the brightness and the flatness of the laser projection line, and compares the brightness and the flatness with standard parameters to judge the qualification degree of the laser head.

Description

一种激光投射线的识别与平面度判定方法A method for identification and flatness determination of laser projection lines

技术领域technical field

本发明涉及激光检测技术领域,尤其涉及一种激光投射线的识别与平面度判定方法。The invention relates to the technical field of laser detection, in particular to a method for identifying and flatness determination of a laser projection line.

背景技术Background technique

激光投线仪,又被称作是激光标线仪或激光水准仪,其实是在普通水准仪望远镜筒上安装并固定了激光装置而制成的一类测量仪器,在使用的过程中,激光投线仪通过发射激光束,使激光束通过棱镜导光系统形成激光面以投射出水平和铅垂的激光线,最终实现测量的目的。其中,具有十二线多功能贴地贴墙激光投线仪是当下最流行的,而判断质量的优劣主要通过激光头来衡量。因此,为了确保工厂能提供合格的激光投线仪,需要在出厂前对激光头进行合格率的检测,以便将不合格的激光头排除掉。Laser line projection instrument, also known as laser marking instrument or laser level, is actually a type of measuring instrument made by installing and fixing a laser device on the telescope barrel of an ordinary level. The instrument emits a laser beam, and makes the laser beam pass through the prism light guide system to form a laser surface to project horizontal and vertical laser lines, and finally achieve the purpose of measurement. Among them, the multi-functional floor-to-wall laser line projection instrument with twelve lines is the most popular at present, and the quality of judgment is mainly measured by the laser head. Therefore, in order to ensure that the factory can provide a qualified laser line projection instrument, it is necessary to test the pass rate of the laser head before leaving the factory, so as to exclude the unqualified laser head.

激光头射出的激光水平偏移角度是当前判断激光头是否合格的重要因素,水平偏移角度在20秒角度以内为合格,大于20秒为不合格。现有的激光头来料检测均为人工检测,人工通过手动转台将激光头调至摆放到固定的位置,通过观察画面中的激光线的相对位置判断是否合格。由于存在人的主观判断、人眼疲劳和手动调整速度慢等因素,可能会导致检测中的误差,同时激光长时间直射到人眼中会发生灼伤等危害,并且检测标准只有一项,检测项目单一。The horizontal offset angle of the laser emitted by the laser head is an important factor in judging whether the laser head is qualified or not. The horizontal offset angle within 20 seconds is qualified, and if it is greater than 20 seconds, it is unqualified. The existing laser head incoming inspection is all manual inspection. Manually adjust the laser head to a fixed position through a manual turntable, and judge whether it is qualified or not by observing the relative position of the laser line in the screen. Due to factors such as human subjective judgment, human eye fatigue and slow manual adjustment, errors in detection may occur. At the same time, if the laser is directly exposed to human eyes for a long time, it will cause burns and other hazards, and there is only one detection standard and a single detection item. .

因此,急需一种技术来解决该问题。Therefore, a technique is urgently needed to solve this problem.

发明内容SUMMARY OF THE INVENTION

本发明的目的在于克服上述现有技术的问题,提供了一种激光投射线的识别与平面度判定方法,通过分析一个平面内4张图像并进行平面拟合等处理,实现判断整个平面的激光线平面度,克服了了以往检测方法中通过对比法进行判断激光线偏移位置的单一性,并且制订检验标准,提高了检测的科学性和准确性。The purpose of the present invention is to overcome the above-mentioned problems of the prior art, and to provide a method for identifying and judging the flatness of a laser projection line. The line flatness overcomes the singleness of judging the offset position of the laser line by the comparison method in the previous detection method, and formulates the inspection standard, which improves the scientificity and accuracy of the detection.

上述目的是通过以下技术方案来实现:The above purpose is achieved through the following technical solutions:

一种激光投射线的识别与平面度判定方法,包括如下步骤:A method for identifying and judging flatness of a laser projection line, comprising the following steps:

步骤(1)以待检测激光头的光源点为中心原点,在所述待检测激光头的四周等距离布置4个平行光管,相邻所述平行光管间的夹角为90°,所述平行光管中分划板上分别设置有十字刻度线,4个所述十字刻度线的中心点与所述光源点位于同一水平面;获取每个所述分划板上十字刻度线图像,得4张刻度线图像;Step (1) takes the light source point of the laser head to be detected as the center origin, and arranges 4 parallel light pipes at equal distances around the laser head to be detected, and the included angle between the adjacent parallel light pipes is 90°, so Cross tick marks are respectively provided on the reticle in the collimator light pipe, and the center points of the 4 cross tick marks are located on the same horizontal plane as the light source point; the image of the reticle tick marks on each of the reticle is obtained to obtain 4 tick images;

步骤(2)对步骤(1)中4张所述刻度线图像分别提取骨架,记录所述十字刻度线中心点的像素坐标和两个最小刻度的像素间距;Step (2) extracts skeletons respectively to 4 described tick marks images in step (1), and records the pixel coordinates of the center point of the cross tick marks and the pixel spacing of two minimum scales;

步骤(3)点亮所述待检测激光头,获取投射在4个所述分划板上的激光线图像,得4张激光线样本图像;Step (3) lighting the to-be-detected laser head, acquiring the laser line images projected on the 4 described reticles, and obtaining 4 laser line sample images;

步骤(4)对步骤(3)中4张所述激光线样本图像提取特征颜色,得激光线图像;Step (4) extracts characteristic colors from 4 described laser line sample images in step (3) to obtain laser line images;

步骤(5)对步骤(4)中所述激光线图像进行分割,提取骨架特征,获取激光线中心线的图像;Step (5) segmenting the laser line image described in step (4), extracting skeleton features, and obtaining an image of the center line of the laser line;

步骤(6)对步骤(4)中所述激光线图像中激光线的外轮廓区域提取亚像素精密轮廓,并对线段的XLD做近似计算直线计算,得到完整的激光轮廓线;Step (6) extracts the sub-pixel precise contour from the outer contour area of the laser line in the laser line image described in the step (4), and performs approximate calculation straight line calculation to the XLD of the line segment to obtain a complete laser contour line;

步骤(7)对步骤(3)所述激光线样本图像的激光线目标区域图像进行分割,计算特征区域的平均灰度值,表示每段激光线的亮度;Step (7) segmenting the laser line target area image of the laser line sample image described in step (3), calculating the average gray value of the characteristic area, indicating the brightness of each laser line;

步骤(8)对步骤(3)中的4个所述激光线样本图像进行阈值分割,将有效像素点,进行坐标转换,以步骤(2)中所述十字刻度线中心点的像素坐标为原点,将像素坐标改为位置坐标,并带入三维空间坐标系,使用拟合空间平面算法,计算平面的离散程度,表示平面度σ;In step (8), threshold segmentation is performed on the four laser line sample images in step (3), and the effective pixel points are converted into coordinates, and the pixel coordinates of the center point of the cross tick mark in step (2) are used as the origin. , change the pixel coordinate to the position coordinate, and bring it into the three-dimensional space coordinate system, and use the fitting space plane algorithm to calculate the discrete degree of the plane, indicating the flatness σ;

步骤(9)对步骤(3)中激光线的区域图像内的所有像素点,模拟偏差20秒时的三维图像,计算偏差20秒时的平面度σ20,对比步骤(8)中的平面度σ,Step (9) For all the pixels in the area image of the laser line in step (3), simulate the three-dimensional image when the deviation is 20 seconds, calculate the flatness σ 20 when the deviation is 20 seconds, and compare the flatness in step (8) σ,

若σ<σ20,则为合格;If σ<σ 20 , it is qualified;

若σ≥σ20,则为不合格。If σ≥σ 20 , it is unqualified.

进一步地,所述步骤(2)具体为:Further, the step (2) is specifically:

设S(X)表示X的骨架,则所述刻度线图像的骨架表达式为:Let S(X) represent the skeleton of X, then the skeleton expression of the tick mark image is:

Figure BDA0003685839540000021
Figure BDA0003685839540000021

Figure BDA0003685839540000022
Figure BDA0003685839540000022

其中,Sn(X)为X的第n个骨架子集,N是

Figure BDA0003685839540000023
运算将X腐蚀成空集前的最后一次迭代次数,
Figure BDA0003685839540000024
Figure BDA0003685839540000025
表示连续n次用B对X进行腐蚀,即where Sn(X) is the nth skeleton subset of X, and N is
Figure BDA0003685839540000023
The last number of iterations before the operation erodes X into an empty set,
Figure BDA0003685839540000024
Figure BDA0003685839540000025
Indicates that X is eroded with B for n consecutive times, that is,

Figure BDA0003685839540000026
Figure BDA0003685839540000026

从Sn(X)中区分出水平刻度线和竖直刻度线,经对其角度进行判断,设degn为第n个骨架的角度,若degn>5°,则为竖直刻度线,用SVn(X)表示,若degn<5°,则为水平刻度线,用SHn(X)表示;遍历所有刻度线,取最长的竖直刻度线SHmax(X)和最长的水平刻度线SVmax(X),记录他们的交点坐标(u0i,v0i),其中(u0,v0)为十字刻度线中心点的像素坐标,i为相机编号。Distinguish the horizontal scale line and the vertical scale line from S n (X), after judging its angle, let deg n be the angle of the nth skeleton, if deg n >5°, it is the vertical scale line, It is represented by SV n (X), if deg n < 5°, it is a horizontal scale line, represented by SH n (X); traverse all the scale lines, take the longest vertical scale line SH max (X) and the longest The horizontal scale line SV max (X), record their intersection coordinates (u 0i , v 0i ), where (u 0 , v 0 ) is the pixel coordinate of the center point of the cross scale line, and i is the camera number.

进一步地,步骤(4)中所述激光线样本图像提取特征颜色,具体为:在RGB颜色模型中R、G、B颜色分量互不影响,因激光主线颜色为绿色,将G分量分离并转化成灰度图,有利于进行颜色特征的提取。Further, the characteristic color of the laser line sample image extraction described in step (4) is specifically: in the RGB color model, the R, G, B color components do not affect each other, and because the laser main line color is green, the G component is separated and converted. A grayscale image is useful for color feature extraction.

进一步地,在所述步骤(7)中,设定Ni为区域内像素点的个数,则计算激光亮度Li的公式为:Further, in the step (7), set Ni as the number of pixel points in the area, then the formula for calculating the laser brightness Li is:

Figure BDA0003685839540000031
Figure BDA0003685839540000031

其中,i为相机编号,j为每个像素的编号。where i is the camera number and j is the number of each pixel.

进一步地,所述步骤(8)具体包括如下步骤:Further, the step (8) specifically includes the following steps:

步骤(8-1)建立空间坐标系,以(0,0,0)为对称中点,设设置于所述平行光管后端用于拍照的相机距离光源点的距离为D,则每个方向上的激光平面距离到对称中点的距离为D;以每个分划板中心交点(u0i,v0i)为坐标原点,对区域内的像素点(uj,vj)进行坐标转换,引入转为相对与原点的位置坐标(xij,yij,zij),其中(u0,v0)为十字刻度线中心点的像素坐标;Step (8-1) establishes a space coordinate system, takes (0, 0, 0) as the symmetrical midpoint, and assumes that the distance between the camera and the light source point arranged at the rear end of the collimated light pipe for taking pictures is D, then each The distance from the laser plane in the direction to the symmetrical midpoint is D; take the center intersection (u 0i , v 0i ) of each reticle as the coordinate origin, and perform coordinate transformation on the pixel points (u j , v j ) in the area , introduce the position coordinates (x ij , y ij , z ij ) relative to the origin, where (u 0 , v 0 ) are the pixel coordinates of the center point of the cross tick mark;

步骤(8-2)建立1号相机对应的正面图坐标,如下式:Step (8-2) Establish the front view coordinates corresponding to the No. 1 camera, as follows:

x1j=Dx 1j =D

Figure BDA0003685839540000032
Figure BDA0003685839540000032

Figure BDA0003685839540000033
Figure BDA0003685839540000033

建立2号相机对应的左面图坐标,如下式:Establish the coordinates of the left image corresponding to the No. 2 camera, as follows:

Figure BDA0003685839540000034
Figure BDA0003685839540000034

y2j=-Dy 2j = -D

Figure BDA0003685839540000035
Figure BDA0003685839540000035

建立3号相机对应的后面图坐标,如下式:Establish the coordinates of the rear image corresponding to camera 3, as follows:

x3j=-Dx 3j = -D

Figure BDA0003685839540000036
Figure BDA0003685839540000036

Figure BDA0003685839540000041
Figure BDA0003685839540000041

建立4号相机对应的右面图坐标,如下式:Establish the coordinates of the right image corresponding to camera No. 4, as follows:

Figure BDA0003685839540000042
Figure BDA0003685839540000042

y4j=-Dy 4j = -D

Figure BDA0003685839540000043
Figure BDA0003685839540000043

步骤(8-3)通过三维点云,将所有相机获取到的有效像素点(xij,yij,zij)放置到三维坐标系中,利用最小二乘法进行平面拟合,计算平面离散程度,得平面度σ。Step (8-3) Through the three-dimensional point cloud, the effective pixels (x ij , y ij , z ij ) obtained by all cameras are placed in the three-dimensional coordinate system, and the least squares method is used to perform plane fitting to calculate the degree of plane dispersion. , the flatness σ is obtained.

进一步地,所述步骤(8-3)具体为:Further, described step (8-3) is specifically:

通过三维点云,将所有相机获取到的阈值内的像素点(xij,yij,zij)放置到三维坐标系中,定义空间中的平面方程为Ax+By+Cz=0,此处定义标准空间平面方程为式如下:Through the three-dimensional point cloud, the pixel points (x ij , y ij , z ij ) within the threshold obtained by all cameras are placed in the three-dimensional coordinate system, and the plane equation in the defined space is Ax+By+Cz=0, here The standard space plane equation is defined as the following formula:

p(x,y,z)=ax-by+cz+1=0p(x, y, z)=ax-by+cz+1=0

式中x,y,z分别为3个坐标轴,a,b,c分别为空间平面方程的3个系数;In the formula, x, y, z are the three coordinate axes, a, b, c are the three coefficients of the space plane equation;

设点集M{(x1,y1,z1),(x2,y2,z2),...,(xn,yn,zn)},其最佳拟合平面满足下式:Let the point set M{(x 1 , y 1 , z 1 ), (x 2 , y 2 , z 2 ), ..., (x n , y n , z n )}, its best fitting plane satisfy The following formula:

Figure BDA0003685839540000044
Figure BDA0003685839540000044

式中:p(x,y,z)=0,则:In the formula: p(x, y, z)=0, then:

Figure BDA0003685839540000045
Figure BDA0003685839540000045

分别对a、b、c求0偏导数,且满足

Figure BDA0003685839540000046
成立,整理后得下式:Find 0 partial derivatives with respect to a, b, and c respectively, and satisfy
Figure BDA0003685839540000046
is established, the following formula is obtained after sorting:

Figure BDA0003685839540000047
Figure BDA0003685839540000047

Assume

Figure BDA0003685839540000048
Figure BDA0003685839540000048

则:but:

QX=KQX=K

解得平面方程系数为:The coefficients of the solved plane equation are:

X=Q-1KX=Q -1 K

将求解的a,b,c带入平面方程,在三维空间坐标系画出平面Ax+By+Cz=0。Bring the solved a, b, c into the plane equation, and draw the plane Ax+By+Cz=0 in the three-dimensional space coordinate system.

平面度表示方法用距离标准差表示,如下所示:The flatness representation method is expressed in distance standard deviation as follows:

Figure BDA0003685839540000051
Figure BDA0003685839540000051

其中,di为有效像素点到拟合平面的距离;davg为有效像素点到拟合平面的平均距离。Among them, d i is the distance from the effective pixel point to the fitting plane; davg is the average distance from the effective pixel point to the fitting plane.

有益效果beneficial effect

本发明所提供的一种激光投射线的识别与平面度判定方法,通过提取激光投射线的颜色、形状和纹理特征,运用平面拟合的方法进行处理,同时加入了激光亮度和激光杂光条纹检测,基于这种多特征融合的方式,使得检测激光头更为精准。本方法节约了人工调整的时间,并且通过对不同激光头进行大量的数据采集,得到科学的标准参数,采用算法精准区分出合格与不合格的激光头,为提高激光头来料入库的检测做出一定贡献。The method for identifying and judging the flatness of a laser projection line provided by the present invention extracts the color, shape and texture features of the laser projection line, uses the method of plane fitting for processing, and adds laser brightness and laser stray light stripes at the same time. Detection, based on this multi-feature fusion method, makes the detection of the laser head more accurate. This method saves the time of manual adjustment, and obtains scientific standard parameters by collecting a large amount of data for different laser heads, and uses an algorithm to accurately distinguish qualified and unqualified laser heads, in order to improve the detection of incoming materials for laser heads. make some contribution.

附图说明Description of drawings

图1为本发明所述一种激光投射线的识别与平面度判定方法的流程图;1 is a flowchart of a method for identifying and flatness determination of a laser projection line according to the present invention;

图2为本发明所述一种激光投射线的识别与平面度判定方法中分划板上十字刻度线图像示意图;2 is a schematic diagram of an image of a cross scale line on a reticle in a method for identifying and flatness determination of a laser projection line according to the present invention;

图3为本发明所述一种激光投射线的识别与平面度判定方法中预处理后提取的分划板上十字刻度线图像示意图;3 is a schematic diagram of an image of a cross tick mark on a reticle extracted after preprocessing in a method for identifying and flatness judging of a laser projection line according to the present invention;

图4为本发明所述一种激光投射线的识别与平面度判定方法中激光投射线轮廓提取示意图;4 is a schematic diagram of contour extraction of a laser projection line in a method for identifying and flatness determination of a laser projection line according to the present invention;

图5为本发明所述一种激光投射线的识别与平面度判定方法中4张激光投射线图像拟合平面图;5 is a fitting plan view of four laser projection line images in a method for identifying and flatness determination of a laser projection line according to the present invention;

图6为本发明所述一种激光投射线的识别与平面度判定方法中最终分段显示亮度结果示意图。FIG. 6 is a schematic diagram of final segmented display brightness results in a method for identifying and flatness determination of a laser projection line according to the present invention.

具体实施方式Detailed ways

下面根据附图和实施例对本发明作进一步详细说明。所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。The present invention will be described in further detail below according to the accompanying drawings and embodiments. The described embodiments are only some, but not all, embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

如图1和2所示,一种激光投射线的识别与平面度判定方法,包括如下步骤:As shown in Figures 1 and 2, a method for identifying and judging flatness of a laser projection line includes the following steps:

步骤(1)以待检测激光头的光源点为中心原点,在所述待检测激光头的四周等距离布置4个平行光管,相邻所述平行光管间的夹角为90°,所述平行光管中分划板上分别设置有十字刻度线,4个所述十字刻度线的中心点与所述光源点位于同一水平面;获取每个所述分划板上十字刻度线图像,得4张刻度线图像;具体的,本实施例中所述平行光管中设置有分划板、光源和目镜,在目镜的后端设置有用于拍照的相机,4个平行光管分别对应1号相机、2号相机、3号相机和4号相机;处理分划板图像时,需要打开光源,在明亮环境下获取分划板图像(即刻度线图像);结合Halcon18库实时处理分划板和激光投射线图像,获取激光投射线的亮度信息并计算平面度,激光投射线图像样本来自黑暗环境下平行拍摄的相机;Step (1) takes the light source point of the laser head to be detected as the center origin, and arranges 4 parallel light pipes at equal distances around the laser head to be detected, and the included angle between the adjacent parallel light pipes is 90°, so Cross tick marks are respectively provided on the reticle in the collimator light pipe, and the center points of the 4 cross tick marks are located on the same horizontal plane as the light source point; the image of the reticle tick marks on each of the reticle is obtained to obtain 4 scale line images; specifically, the collimating light pipe described in this embodiment is provided with a reticle, a light source and an eyepiece, and a camera for taking pictures is provided at the rear end of the eyepiece, and the four collimating light pipes correspond to No. 1 respectively. Camera, Camera 2, Camera 3, and Camera 4; when processing reticle images, you need to turn on the light source and obtain reticle images (ie tick images) in a bright environment; combine the Halcon18 library to process reticle and reticle images in real time. Laser projection line image, obtain the brightness information of the laser projection line and calculate the flatness, the laser projection line image sample is from a camera that is photographed in parallel in a dark environment;

本实施例中待检测激光头通过一个三维角位转台进行角度的调节,以便将其发射的激光线很好的投射到分划板上;所述三维角位转台通过三个伺服电机控制,包括控制X向转角的第一伺服电机、控制Y向转角的第二伺服电机和控制夹具(用于夹持或固定待检测激光)做旋转的第三伺服电机,三个伺服电机于伺服控制连接。In this embodiment, the angle of the laser head to be detected is adjusted through a three-dimensional angular turntable, so that the laser line emitted by it can be well projected onto the reticle; the three-dimensional angular turntable is controlled by three servo motors, including The first servo motor that controls the X-direction rotation angle, the second servo motor that controls the Y-direction rotation angle, and the third servo motor that controls the fixture (used to clamp or fix the laser to be detected) for rotation, and the three servo motors are connected to the servo control.

步骤(2)对步骤(1)中4张所述刻度线图像分别提取骨架,记录所述十字刻度线中心点的像素坐标和两个最小刻度的像素间距;本实施例中所述十字刻度线中心点的像素坐标记为(u0,v0),记两个最小刻度的像素间距为Dmin,实际分划板一格刻度大小为0.1mm,即0.1mm对应的像素距离是DminIn step (2), the skeletons are respectively extracted from the four tick marks in step (1), and the pixel coordinates of the center point of the tick marks and the pixel spacing of the two minimum scales are recorded; the tick marks in this embodiment are The pixel coordinates of the center point are marked as (u 0 , v 0 ), and the pixel spacing between the two smallest scales is recorded as D min , and the actual reticle scale size is 0.1 mm, that is, the pixel distance corresponding to 0.1 mm is D min ;

具体的,设S(X)表示X的骨架,则所述刻度线图像的骨架表达式为:Specifically, let S(X) represent the skeleton of X, then the skeleton expression of the tick mark image is:

Figure BDA0003685839540000061
Figure BDA0003685839540000061

Figure BDA0003685839540000062
Figure BDA0003685839540000062

其中,Sn(X)为X的第n个骨架子集,N是

Figure BDA0003685839540000063
运算将X腐蚀成空集前的最后一次迭代次数,
Figure BDA0003685839540000064
Figure BDA0003685839540000065
表示连续n次用B对X进行腐蚀,即where Sn(X) is the nth skeleton subset of X, and N is
Figure BDA0003685839540000063
The last number of iterations before the operation erodes X into an empty set,
Figure BDA0003685839540000064
Figure BDA0003685839540000065
Indicates that X is eroded with B for n consecutive times, that is,

Figure BDA0003685839540000066
Figure BDA0003685839540000066

从Sn(X)中区分出水平刻度线和竖直刻度线,经对其角度进行判断,设degn为第n个骨架的角度,若degn>5°,则为竖直刻度线,用SVn(X)表示,若degn<5°,则为水平刻度线,用SHn(X)表示;遍历所有刻度线,取最长的竖直刻度线SHmax(X)和最长的水平刻度线SVmax(X),记录他们的交点坐标(u0i,v0i),其中(u0,v0)为十字刻度线中心点的像素坐标,i为相机编号;Distinguish the horizontal scale line and the vertical scale line from S n (X), after judging its angle, let deg n be the angle of the nth skeleton, if deg n >5°, it is the vertical scale line, It is represented by SV n (X), if deg n < 5°, it is a horizontal scale line, represented by SH n (X); traverse all the scale lines, take the longest vertical scale line SH max (X) and the longest The horizontal scale line SV max (X), record their intersection coordinates (u 0i , v 0i ), where (u 0 , v 0 ) is the pixel coordinate of the center point of the cross scale line, and i is the camera number;

步骤(3)点亮所述待检测激光头,获取投射在4个所述分划板上的激光线图像,得4张激光线样本图像;所述激光线样本图像来自黑暗环境下平行拍摄的相机;Step (3) lighting the laser head to be detected, acquiring the laser line images projected on the 4 reticles, and obtaining 4 laser line sample images; the laser line sample images are from parallel shooting in a dark environment camera;

步骤(4)对步骤(3)中4张所述激光线样本图像提取特征颜色,得激光线图像;Step (4) extracts characteristic colors from 4 described laser line sample images in step (3) to obtain laser line images;

具体的,RGB颜色模型中的R、G、B颜色分量互不影响,这有利于对图像的颜色和亮度进行处理。因激光主线颜色为绿色,故将G分量分离并转化成灰度图,这样将更加有利于进行颜色特征的提取。Specifically, the R, G, and B color components in the RGB color model do not affect each other, which is beneficial for processing the color and brightness of the image. Since the color of the main line of the laser is green, the G component is separated and converted into a grayscale image, which is more conducive to the extraction of color features.

步骤(5)对步骤(4)中所述激光线图像进行分割,提取骨架特征,获取激光线中心线的图像;Step (5) segmenting the laser line image described in step (4), extracting skeleton features, and obtaining an image of the center line of the laser line;

步骤(6)对步骤(4)中所述激光线图像中激光线的外轮廓区域提取亚像素精密轮廓,并对线段的XLD做近似计算直线计算,得到完整的激光轮廓线;Step (6) extracts the sub-pixel precise contour from the outer contour area of the laser line in the laser line image described in the step (4), and performs approximate calculation straight line calculation to the XLD of the line segment to obtain a complete laser contour line;

步骤(7)对步骤(3)所述激光线样本图像的激光线目标区域图像进行分割,计算特征区域的平均灰度值,表示每段激光线的亮度;Step (7) segmenting the laser line target area image of the laser line sample image described in step (3), calculating the average gray value of the characteristic area, indicating the brightness of each laser line;

在本步骤中中,设定Ni为区域内像素点的个数,则计算激光亮度Li的公式为:In this step, set Ni as the number of pixels in the area, then the formula for calculating the laser brightness Li is:

Figure BDA0003685839540000071
Figure BDA0003685839540000071

其中,i为相机编号,j为每个像素的编号。where i is the camera number and j is the number of each pixel.

步骤(8)对步骤(3)中的4个所述激光线样本图像进行阈值分割,将有效像素点,进行坐标转换,以步骤(2)中所述十字刻度线中心点的像素坐标为原点,将像素坐标改为位置坐标,并带入三维空间坐标系,使用拟合空间平面算法,计算平面的离散程度,表示平面度σ;In step (8), threshold segmentation is performed on the four laser line sample images in step (3), and the effective pixel points are converted into coordinates, and the pixel coordinates of the center point of the cross tick mark in step (2) are used as the origin. , change the pixel coordinate to the position coordinate, and bring it into the three-dimensional space coordinate system, and use the fitting space plane algorithm to calculate the discrete degree of the plane, indicating the flatness σ;

具体的,本步骤具体包括如下步骤:Specifically, this step specifically includes the following steps:

步骤(8-1)建立空间坐标系,以(0,0,0)为对称中点,设设置于所述平行光管后端用于拍照的相机距离光源点的距离为D,则每个方向上的激光平面距离到对称中点的距离为D;以每个分划板中心交点(u0i,v0i)为坐标原点,对区域内的像素点(uj,vj)进行坐标转换,引入转为相对与原点的位置坐标(xij,yij,zij),其中(u0,v0)为十字刻度线中心点的像素坐标;Step (8-1) establishes a space coordinate system, takes (0, 0, 0) as the symmetrical midpoint, and assumes that the distance between the camera and the light source point arranged at the rear end of the collimated light pipe for taking pictures is D, then each The distance from the laser plane in the direction to the symmetrical midpoint is D; take the center intersection (u 0i , v 0i ) of each reticle as the coordinate origin, and perform coordinate transformation on the pixel points (u j , v j ) in the area , introduce the position coordinates (x ij , y ij , z ij ) relative to the origin, where (u 0 , v 0 ) are the pixel coordinates of the center point of the cross tick mark;

步骤(8-2)建立1号相机对应的正面图坐标,如下式:Step (8-2) Establish the front view coordinates corresponding to the No. 1 camera, as follows:

x1j=Dx 1j =D

Figure BDA0003685839540000072
Figure BDA0003685839540000072

Figure BDA0003685839540000073
Figure BDA0003685839540000073

建立2号相机对应的左面图坐标,如下式:Establish the coordinates of the left image corresponding to the No. 2 camera, as follows:

Figure BDA0003685839540000081
Figure BDA0003685839540000081

y2j=-Dy 2j = -D

Figure BDA0003685839540000082
Figure BDA0003685839540000082

建立3号相机对应的后面图坐标,如下式:Establish the coordinates of the rear image corresponding to camera 3, as follows:

x3j=-Dx 3j = -D

Figure BDA0003685839540000083
Figure BDA0003685839540000083

Figure BDA0003685839540000084
Figure BDA0003685839540000084

建立4号相机对应的右面图坐标,如下式:Establish the coordinates of the right image corresponding to camera No. 4, as follows:

Figure BDA0003685839540000085
Figure BDA0003685839540000085

y4j=-Dy 4j = -D

Figure BDA0003685839540000086
Figure BDA0003685839540000086

步骤(8-3)通过三维点云,将所有相机获取到的有效像素点(xij,yij,zij)放置到三维坐标系中,利用最小二乘法进行平面拟合,计算平面离散程度,得平面度σ。Step (8-3) Through the three-dimensional point cloud, the effective pixels (x ij , y ij , z ij ) obtained by all cameras are placed in the three-dimensional coordinate system, and the least squares method is used to perform plane fitting to calculate the degree of plane dispersion. , the flatness σ is obtained.

进一步地,所述步骤(8-3)具体为:Further, described step (8-3) is specifically:

通过三维点云,将所有相机获取到的阈值内的像素点(xij,yij,zij)放置到三维坐标系中,定义空间中的平面方程为Ax+By+Cz=0,此处定义标准空间平面方程为式如下:Through the three-dimensional point cloud, the pixel points (x ij , y ij , z ij ) within the threshold obtained by all cameras are placed in the three-dimensional coordinate system, and the plane equation in the defined space is Ax+By+Cz=0, here The standard space plane equation is defined as the following formula:

p(x,y,z)=ax-by+cz+1=0p(x, y, z)=ax-by+cz+1=0

式中x,y,z分别为3个坐标轴,a,b,c分别为空间平面方程的3个系数;In the formula, x, y, z are the three coordinate axes, a, b, c are the three coefficients of the space plane equation;

设点集M{(x1,y1,z1),(x2,y2,z2),...,(xn,yn,zn)},其最佳拟合平面满足下式:Let the point set M{(x 1 , y 1 , z 1 ), (x 2 , y 2 , z 2 ), ..., (x n , y n , z n )}, its best fitting plane satisfy The following formula:

Figure BDA0003685839540000087
Figure BDA0003685839540000087

式中:p(x,y,z)=0,则:In the formula: p(x, y, z)=0, then:

Figure BDA0003685839540000088
Figure BDA0003685839540000088

分别对a、b、c求0偏导数,且满足

Figure BDA0003685839540000089
成立,整理后得下式:Find 0 partial derivatives with respect to a, b, and c respectively, and satisfy
Figure BDA0003685839540000089
is established, the following formula is obtained after sorting:

Figure BDA0003685839540000091
Figure BDA0003685839540000091

Assume

Figure BDA0003685839540000092
Figure BDA0003685839540000092

则:but:

QX=KQX=K

解得平面方程系数为:The coefficients of the solved plane equation are:

X=Q-1KX=Q -1 K

将求解的a,b,c带入平面方程,在三维空间坐标系画出平面Ax+By+Cz=0。Bring the solved a, b, c into the plane equation, and draw the plane Ax+By+Cz=0 in the three-dimensional space coordinate system.

平面度表示方法用距离标准差表示,如下所示:The flatness representation method is expressed in distance standard deviation as follows:

Figure BDA0003685839540000093
Figure BDA0003685839540000093

其中,di为有效像素点到拟合平面的距离;davg为有效像素点到拟合平面的平均距离。Among them, d i is the distance from the effective pixel point to the fitting plane; davg is the average distance from the effective pixel point to the fitting plane.

步骤(9)对步骤(3)中激光线的区域图像内的所有像素点,模拟偏差20秒时的三维图像,计算偏差20秒时的平面度σ20,对比步骤(8)中的平面度σ,Step (9) For all the pixels in the area image of the laser line in step (3), simulate the three-dimensional image when the deviation is 20 seconds, calculate the flatness σ 20 when the deviation is 20 seconds, and compare the flatness in step (8) σ,

若σ<σ20,则为合格;If σ<σ 20 , it is qualified;

若σ≥σ20,则为不合格。If σ≥σ 20 , it is unqualified.

由于激光头投射出的激光线应是一圈水平的平面,其主要误差在于激光投射出来的角度误差,其工厂工艺标准为偏差小于20秒,故使用当前相机获取的有效像素点,模拟激光线在分划板上偏差20秒的位置图像,并用此图像的像素点放置到三维空间坐标系中,利用最小二乘法进行平面拟合,计算平面离散程度σ2。将其中一个方向上的图像,带入误差20秒的坐标,这里我们更改的是相机1中的坐标,对应步骤F中的(u01,v01)改为(u01+Dmin,v01),再计算σ2Since the laser line projected by the laser head should be a horizontal plane, the main error lies in the angle error projected by the laser. The factory process standard is that the deviation is less than 20 seconds. Therefore, the effective pixel points obtained by the current camera are used to simulate the laser line. A position image with a deviation of 20 seconds on the reticle, and the pixels of this image are placed in the three-dimensional space coordinate system, and the least square method is used to perform plane fitting to calculate the plane dispersion degree σ 2 . Bring the image in one of the directions into the coordinates with an error of 20 seconds. Here we change the coordinates in camera 1, corresponding to (u 01 , v 01 ) in step F to (u 01 +D min , v 01 ) ), and then calculate σ 2 :

Figure BDA0003685839540000094
Figure BDA0003685839540000094

对比如果σ<σ2,表示水平度合格;若σ>σ2,表示水平度不合格。In contrast, if σ<σ 2 , it means that the levelness is qualified; if σ>σ 2 , it means that the levelness is unqualified.

以上所述仅为说明本发明的实施方式,并不用于限制本发明,对于本领域的技术人员来说,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The above description is only to illustrate the embodiments of the present invention, and is not intended to limit the present invention. For those skilled in the art, any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention, All should be included within the protection scope of the present invention.

Claims (6)

1. A method for identifying a laser projection line and judging flatness is characterized by comprising the following steps:
the method comprises the following steps that (1) a light source point of a laser head to be detected is taken as a central original point, 4 parallel light tubes are arranged on the periphery of the laser head to be detected at equal intervals, an included angle between every two adjacent parallel light tubes is 90 degrees, cross scale marks are respectively arranged on a dividing plate in each parallel light tube, and the central points of the 4 cross scale marks and the light source point are located on the same horizontal plane; acquiring a cross scale line image on each reticle to obtain 4 scale line images;
respectively extracting skeletons from the 4 scale mark images in the step (1), and recording the pixel coordinate of the central point of the cross scale mark and the pixel distance of two minimum scales;
step (3) lightening the laser head to be detected, and acquiring laser line images projected on 4 reticles to obtain 4 laser line sample images;
step (4) extracting characteristic colors of the 4 laser line sample images in the step (3) to obtain laser line images;
step (5) segmenting the laser line image in the step (4), extracting skeleton characteristics and obtaining an image of a laser line central line;
step (6) extracting sub-pixel precise contours from the outer contour region of the laser line in the laser line image in the step (4), and performing approximate calculation linear calculation on XLD of the line segment to obtain a complete laser contour line;
step (7) dividing the laser line target area image of the laser line sample image in the step (3), calculating the average gray value of the characteristic area, and expressing the brightness of each section of laser line;
step (8) performing threshold segmentation on the 4 laser line sample images in the step (3), performing coordinate conversion on effective pixel points, changing pixel coordinates into position coordinates by taking the pixel coordinates of the central point of the cross scale line in the step (2) as an original point, and bringing the position coordinates into a three-dimensional space coordinate system, and calculating the dispersion degree of a plane by using a fitting space plane algorithm to express the flatness sigma;
step (9) simulating a three-dimensional image when the deviation is 20 seconds for all pixel points in the area image of the laser line in the step (3), and calculating the planeness sigma when the deviation is 20 seconds 20 Comparing the flatness sigma in the step (8),
if σ < σ 20 If so, the product is qualified;
if sigma is larger than or equal to sigma 20 And then the test is not qualified.
2. The method as claimed in claim 1, wherein the step (2) is specifically as follows:
assuming that S (X) represents the skeleton of X, the skeleton expression of the scale line image is as follows:
Figure FDA0003685839530000011
Figure FDA0003685839530000012
wherein S is n (X) is the nth backbone subset of X, N is
Figure FDA0003685839530000013
The operation erodes X the number of last iterations before empty set,
Figure FDA0003685839530000021
Figure FDA0003685839530000026
denotes etching X with B n successive times, i.e.
Figure FDA0003685839530000022
From S n (X) distinguishing the horizontal scale mark and the vertical scale mark, judging the angle of the horizontal scale mark and the vertical scale mark, and setting deg n Angle of the nth skeleton, if deg n More than 5 degrees, the vertical scale mark is formed, and SV is used n (X) represents, if deg n Less than 5 deg. is horizontal scale mark, using SH n (X) represents; traversing all the scale marks, and taking the longest vertical scale mark SH max (X) and the longest horizontal scale line SV max (X) recording their intersection coordinates (u) 0i ,v 0i ) Wherein (u) 0 ,v 0 ) And i is the pixel coordinate of the center point of the cross scale line, and i is the camera number.
3. The method as claimed in claim 1, wherein the step (4) of extracting the characteristic color from the laser line sample image comprises: r, G, B color components in the RGB color model are not affected each other, and because the laser main line color is green, the G component is separated and converted into a gray scale map, which is beneficial to extracting color features.
4. The method as claimed in claim 1, wherein in the step (7), N is set i Calculating the laser brightness L for the number of pixel points in the region i The formula of (1) is:
Figure FDA0003685839530000023
where i is the camera number and j is the number of each pixel.
5. The method as claimed in claim 1, wherein the step (8) comprises the following steps:
step (8-1) establishing a space coordinate system, taking (0, 0, 0) as a symmetrical midpoint, setting the distance from a camera for taking pictures, which is arranged at the rear end of the collimator tube, to a light source point to be D, and setting the distance from the laser plane distance in each direction to the symmetrical midpoint to be D; with each reticle center intersection (u) 0i ,v 0i ) As the origin of coordinates, to pixel points (u) in the region j ,v j ) Performing coordinate conversion, introducing the position coordinate (x) converted to relative to the origin ij ,y ij ,z ij ) Wherein (u) 0 ,v 0 ) Is the pixel coordinate of the central point of the cross scale mark;
and (8-2) establishing a front view coordinate corresponding to the camera No. 1, wherein the front view coordinate is as follows:
x 1j =D
Figure FDA0003685839530000024
Figure FDA0003685839530000025
establishing the coordinates of the left image corresponding to the camera No. 2 as follows:
Figure FDA0003685839530000031
y 2j =-D
Figure FDA0003685839530000032
establishing the coordinates of the rear graph corresponding to the camera No. 3 as follows:
x 3j =-D
Figure FDA0003685839530000033
Figure FDA0003685839530000034
establishing the coordinates of a right image corresponding to the camera No. 4 as follows:
Figure FDA0003685839530000035
y 4j =-D
Figure FDA0003685839530000036
step (8-3) obtaining effective pixel points (x) obtained by all cameras through three-dimensional point cloud ij ,y ij ,z ij ) And placing the three-dimensional coordinate system in a three-dimensional coordinate system, performing plane fitting by using a least square method, and calculating the plane dispersion degree to obtain the flatness sigma.
6. The method as claimed in claim 4, wherein the step (8-3) is specifically as follows:
obtaining pixel points (x) within a threshold value acquired by all cameras through three-dimensional point cloud ij ,y ij ,z ij ) Put into a three-dimensional coordinate system, the plane equation in the defined space is Ax + By + Cz ═ 0, and here, the standard spatial plane equation is defined as follows:
p(x,y,z)=ax-by+cz+1=0
in the formula, x, y and z are respectively 3 coordinate axes, and a, b and c are respectively 3 coefficients of a space plane equation;
set point set M { (x) 1 ,y 1 ,z 1 ),(x 2 ,y 2 ,z 2 ),...,(x n ,y n ,z n ) The best fit plane of which satisfies the following equation:
Figure FDA0003685839530000037
in the formula: when p (x, y, z) is 0, then:
Figure FDA0003685839530000038
respectively calculating 0 partial derivative for a, b and c, and satisfying
Figure FDA0003685839530000039
If true, the following formula is obtained after finishing:
Figure FDA0003685839530000041
is provided with
Figure FDA0003685839530000042
X=[a,b,c] T
Figure FDA0003685839530000043
Then:
QX=K
the coefficients of the plane equation are solved as follows:
X=Q -1 K
and substituting the solved a, b and c into a plane equation, and drawing a plane Ax + By + Cz as 0 in a three-dimensional space coordinate system.
The flatness representation is expressed in terms of distance standard deviation, as follows:
Figure FDA0003685839530000044
wherein d is i The distance from the effective pixel point to the fitting plane is obtained; d avg And the average distance from the effective pixel point to the fitting plane.
CN202210650078.2A 2022-06-09 2022-06-09 Laser projection line identification and flatness judgment method Active CN115014248B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210650078.2A CN115014248B (en) 2022-06-09 2022-06-09 Laser projection line identification and flatness judgment method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210650078.2A CN115014248B (en) 2022-06-09 2022-06-09 Laser projection line identification and flatness judgment method

Publications (2)

Publication Number Publication Date
CN115014248A true CN115014248A (en) 2022-09-06
CN115014248B CN115014248B (en) 2024-02-06

Family

ID=83073811

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210650078.2A Active CN115014248B (en) 2022-06-09 2022-06-09 Laser projection line identification and flatness judgment method

Country Status (1)

Country Link
CN (1) CN115014248B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116725730A (en) * 2023-08-11 2023-09-12 北京市农林科学院智能装备技术研究中心 A pig vaccine injection method, system and storage medium based on visual guidance
CN117053687A (en) * 2023-08-17 2023-11-14 广州市西克传感器有限公司 Cell height level difference detection method based on laser line scanning 3D camera
CN118049937A (en) * 2024-01-24 2024-05-17 江苏锡滚轴承科技有限公司 Bearing roundness detection device and detection method
CN118376217A (en) * 2024-05-31 2024-07-23 南京景曜智能科技有限公司 Laser positioning projection method

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1571398A (en) * 1976-06-05 1980-07-16 Heidenhain Gmbh Dr Johannes Arrangement for ascertaining deviations from straightness or planeness
JP2013076623A (en) * 2011-09-30 2013-04-25 Avanstrate Inc Flatness measuring method
CN206223112U (en) * 2016-09-29 2017-06-06 七海测量技术(深圳)有限公司 Flatness quick detection equipment based on many laser feelers
CN112958954A (en) * 2021-02-03 2021-06-15 北京博清科技有限公司 Welding seam tracking method and system
CN113137941A (en) * 2021-04-29 2021-07-20 荣旗工业科技(苏州)股份有限公司 Method for testing product flatness based on point laser

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1571398A (en) * 1976-06-05 1980-07-16 Heidenhain Gmbh Dr Johannes Arrangement for ascertaining deviations from straightness or planeness
JP2013076623A (en) * 2011-09-30 2013-04-25 Avanstrate Inc Flatness measuring method
CN206223112U (en) * 2016-09-29 2017-06-06 七海测量技术(深圳)有限公司 Flatness quick detection equipment based on many laser feelers
CN112958954A (en) * 2021-02-03 2021-06-15 北京博清科技有限公司 Welding seam tracking method and system
CN113137941A (en) * 2021-04-29 2021-07-20 荣旗工业科技(苏州)股份有限公司 Method for testing product flatness based on point laser

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘许亮: "基于优化点激光的平面度测量研究", 《应用激光》, vol. 42, no. 5, pages 88 - 94 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116725730A (en) * 2023-08-11 2023-09-12 北京市农林科学院智能装备技术研究中心 A pig vaccine injection method, system and storage medium based on visual guidance
CN116725730B (en) * 2023-08-11 2023-12-05 北京市农林科学院智能装备技术研究中心 Pig vaccine injection method, system and storage medium based on visual guidance
CN117053687A (en) * 2023-08-17 2023-11-14 广州市西克传感器有限公司 Cell height level difference detection method based on laser line scanning 3D camera
CN117053687B (en) * 2023-08-17 2024-06-07 广东西克智能科技有限公司 Cell height level difference detection method based on laser line scanning 3D camera
CN118049937A (en) * 2024-01-24 2024-05-17 江苏锡滚轴承科技有限公司 Bearing roundness detection device and detection method
CN118376217A (en) * 2024-05-31 2024-07-23 南京景曜智能科技有限公司 Laser positioning projection method

Also Published As

Publication number Publication date
CN115014248B (en) 2024-02-06

Similar Documents

Publication Publication Date Title
CN115014248B (en) Laser projection line identification and flatness judgment method
CN108764257B (en) Multi-view pointer instrument identification method
CN106949848B (en) A kind of high-precision laser 3D profile phone structural detection method
CN112818988B (en) Automatic identification reading method and system for pointer instrument
CN109520436B (en) A machine vision-based automatic measuring system for three-dimensional size of butterfly spring and its measuring method
CN109900711A (en) Workpiece defect detection method based on machine vision
CN108335286B (en) Online weld joint forming visual detection method based on double-line structured light
US9052253B2 (en) Method of determining at least one refraction characteristic of an ophthalmic lens
CN106340044A (en) Camera external parameter automatic calibration method and calibration device
CN104913737A (en) Component quality checking device based on line laser three-dimensional measurement and detection method of device
CN109489591A (en) Plane scratch length non-contact measurement method based on machine vision
CN112767364B (en) Image detection system for gate blade surface corrosion and corrosion area rapid measuring and calculating method
CN110261069B (en) Detection method for optical lens
CN110207951A (en) A kind of aircraft cable support assembled state detection method of view-based access control model
US20200388019A1 (en) Method and system for testing field of view
CN111179335A (en) Standing tree measuring method based on binocular vision
CN113008158B (en) Multi-line laser tire pattern depth measuring method
CN112881424A (en) AI + fluorescence permeation small-sized pipe surface defect detection and quality grading method and system
CN112902869B (en) Laser plane adjustment method and device for orbital profile measurement system
CN103884294B (en) Method and device for three-dimensional shape measurement with infrared light with wide field of view
CN112361989A (en) Method for calibrating parameters of measurement system through point cloud uniformity consideration
CN118837379A (en) High-precision defect detection method and system for intelligent glasses
CN113435442A (en) Water level measuring method, device, water gauge and electronic equipment
CN117419663A (en) Rail geometry detecting system
CN117974804A (en) Calibration data verification method and device, acquisition equipment and computer storage medium

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