[go: up one dir, main page]

CN105783782A - Surface curvature abrupt change optical contour measurement method - Google Patents

Surface curvature abrupt change optical contour measurement method Download PDF

Info

Publication number
CN105783782A
CN105783782A CN201610311407.5A CN201610311407A CN105783782A CN 105783782 A CN105783782 A CN 105783782A CN 201610311407 A CN201610311407 A CN 201610311407A CN 105783782 A CN105783782 A CN 105783782A
Authority
CN
China
Prior art keywords
light band
band image
image
curvature
average gray
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
CN201610311407.5A
Other languages
Chinese (zh)
Other versions
CN105783782B (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201610311407.5A priority Critical patent/CN105783782B/en
Publication of CN105783782A publication Critical patent/CN105783782A/en
Application granted granted Critical
Publication of CN105783782B publication Critical patent/CN105783782B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Image Analysis (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Image Processing (AREA)

Abstract

The present invention discloses a surface curvature abrupt change optical contour measurement method. According to the method, an optical band image is obtained; mean filtering and binaryzation processing are performed on the optical band image; average gray calculation is performed on the binarized optical band image; an average gray algorithm is adopted to determine the special position of the curvature abrupt change of a step; in a 3D contour reconstruction process, corresponding recognition processing can be performed on a curvature abrupt change position, and the contour line of the curvature abrupt change position can be filtered out; and the contour of an object of which the surface is subjected to curvature abrupt change can be restored accurately.

Description

表面曲率突变光学轮廓测量方法Optical profilometry method for surface curvature abrupt change

技术领域technical field

本发明属于图像处理技术领域,具体涉及一种表面曲率突变光学轮廓测量方法。The invention belongs to the technical field of image processing, and in particular relates to an optical profile measurement method for sudden changes in surface curvature.

背景技术Background technique

随着机械、汽车、航天等制造业和服务、制鞋、玩具等民用工业的发展,在模具制造业中对复杂表面进行测量的需要越来越多。迄今为止,光学三维轮廓测量技术已经有了很长时间的发展,已有许多成熟的测量方法。但现有的三维轮廓系统的测量精度、自动化、自适应程度普遍偏低,尤其对复杂结构轮廓的测量,现尚存在许多没有解决的难题。With the development of manufacturing industries such as machinery, automobiles, and aerospace, and civil industries such as service, shoemaking, and toys, there are more and more needs to measure complex surfaces in the mold manufacturing industry. So far, optical three-dimensional profilometry technology has been developed for a long time, and there are many mature measurement methods. However, the measurement accuracy, automation, and self-adaptability of the existing 3D contour systems are generally low, especially for the measurement of complex structure contours, there are still many unsolved problems.

三维测量是各个应用领域,如制造业、检验、文档分析、工程设计、刑事侦查现场痕迹分析、自动在线检测、质量控制、机器人、医疗诊断和军事等领域中各种智能自主系统中不可分割的一部分。随着现代检测技术的进步,三维测量技术逐步成为人们的研究重点,特别是随着激光技术、计算机技术以及图像处理技术等高新技术的发展,使得光学式三维测量技术得到广泛的应用。Three-dimensional measurement is an integral part of various intelligent autonomous systems in various application fields, such as manufacturing, inspection, document analysis, engineering design, criminal investigation scene trace analysis, automatic online inspection, quality control, robotics, medical diagnosis and military. part. With the advancement of modern detection technology, three-dimensional measurement technology has gradually become the focus of people's research, especially with the development of high and new technologies such as laser technology, computer technology and image processing technology, optical three-dimensional measurement technology has been widely used.

而随着计算机技术和光电技术的发展,基于光学原理、以计算机图像处理为主要手段的三维复杂曲面非接触式快速测量技术也应运而生。逆向工程,是指从实物上采集大量的三维坐标点,并由此建立该物体的几何模型,进而开发出同类产品的先进技术,复杂曲面的逆向工程也是逆向工程的难点之一,它是根据曲面的测量规划,对曲面进行数字化扫描检测,再将数字点云进行曲面重构生成CAD模型,随着机械、汽车、航天等制造业和服务、制鞋、玩具等民用工业的发展,在模具制造业中对复杂表面进行测量的需要越来越多。With the development of computer technology and optoelectronic technology, the non-contact rapid measurement technology of three-dimensional complex curved surface based on optical principle and computer image processing as the main method has also emerged as the times require. Reverse engineering refers to the collection of a large number of three-dimensional coordinate points from the real object, and the establishment of the geometric model of the object, and then the development of advanced technology for similar products. The reverse engineering of complex surfaces is also one of the difficulties of reverse engineering. It is based on Surface measurement planning, digital scanning and detection of the surface, and then surface reconstruction of the digital point cloud to generate a CAD model. The need to measure complex surfaces in manufacturing is increasing.

光学三维轮廓测量由于其非接触性、高精度与高分辨率,在CAD/CAE、反求工程、在线检测与质量保证、多媒体技术、医疗诊断、机器视觉等领域得到日益广泛的应用,被公认是最有前途的三维轮廓测量方法。Due to its non-contact, high precision and high resolution, optical three-dimensional profile measurement has been widely used in CAD/CAE, reverse engineering, online inspection and quality assurance, multimedia technology, medical diagnosis, machine vision and other fields, and is recognized is the most promising method for 3D profile measurement.

对于曲率突变的几何物体,一般的三维测量方法都会失效。例如线结构激光法,当线形激光器光刀照射于曲率突变的表面时,受表面突变影响,激光带会变成无大面积无规则形状,从而使轮廓测量失败。For geometric objects with sudden changes in curvature, general three-dimensional measurement methods will fail. For example, in the line-structured laser method, when the line-shaped laser light knife is irradiated on the surface with a sudden change in curvature, affected by the sudden change in the surface, the laser band will become irregular without a large area, thus making the contour measurement fail.

发明内容Contents of the invention

本发明的目的是解决上述问题,提供一种表面曲率突变光学轮廓测量方法,可准确的判断几何物体曲率突变的位置并作出相应的处理,使得轮廓测量不受影响。The purpose of the present invention is to solve the above problems and provide an optical profile measurement method for sudden changes in surface curvature, which can accurately determine the position of sudden changes in the curvature of geometric objects and perform corresponding processing, so that the profile measurement will not be affected.

为解决上述技术问题,本发明的技术方案是:一种表面曲率突变光学轮廓测量方法,包括以下步骤:In order to solve the above technical problems, the technical solution of the present invention is: a method for measuring optical profilometry with abrupt changes in surface curvature, comprising the following steps:

S1:待测物处于位置j时,利用计算机软件发出脉冲调制信号对激光器进行调制,激光器输出线型光束对待测物进行照明;同时利用计算机软件发出触发信号,触发CCD图像传感器同步采集待测物表面的光带图像,并将采集的图像记为A1,然后关闭激光器,CCD摄像机再次采集相同位置的图像,记为A11S1: When the object under test is at position j, use computer software to send pulse modulation signals to modulate the laser, and the laser outputs a linear beam to illuminate the object under test; at the same time, use computer software to send a trigger signal to trigger the CCD image sensor to synchronously collect the object under test The light band image on the surface, and record the collected image as A 1 , then turn off the laser, and the CCD camera collects the image at the same position again, and record it as A 11 ;

S2:对由锁定成像所得的图像A1、A11进行处理,在此位置获得一幅去除背景光的光带图像IjS2: Process the images A 1 and A 11 obtained by locking imaging, and obtain a light band image I j at this position with the background light removed;

S3:对上述S2中获得的光带图像Ij进行二值化处理;S3: Binarize the light band image I j obtained in the above S2;

S4:对二值化处理后的光带图像Ij进行平均灰度计算,设该图片中有N个像素点,第i个像素点的灰度值为Ni,则光带图像Ij的平均灰度为: S4: Calculate the average gray level of the light band image I j after binarization, assuming that there are N pixels in the picture, and the gray value of the i-th pixel is N i , then the light band image I j The average grayscale is:

S5:将待测物从位置j沿导轨平移距离L到达位置j+1处,重复步骤S1至S4,得到与光带图像I位置连续的光带图像Ij+1,并计算出光带图像Ij+1的平均灰度 S5: Translate the object to be measured from position j along the guide rail to the position j+1, repeat steps S1 to S4, obtain the light band image I j+1 continuous with the position of the light band image I, and calculate the light band image I Average gray level of j+1

S6:将光带图像Ij+1的平均灰度除以光带图像Ij的平均灰度,得到比值进行判断,若比值则光带图像Ij+1为曲率突变图像,该位置为曲率突变位置,在待测物的每个位置,计算机软件均采用三维表面轮廓重建技术进行同步的轮廓重建,若该位置为曲率突变位置,则滤出该位置轮廓线后返回S5;若比值则该位置没有曲率突变,返回S5;直到完成待测物整体的轮廓重建。S6: Divide the average gray level of the light band image I j +1 by the average gray level of the light band image I j to obtain the ratio Make a judgment, if the ratio Then the light band image I j+1 is a curvature mutation image, and this position is the curvature mutation position. At each position of the object to be measured, the computer software uses the three-dimensional surface contour reconstruction technology to perform simultaneous contour reconstruction. If the position is a curvature mutation position, then filter out the contour line of the position and return to S5; if the ratio Then there is no abrupt curvature change at this position, and return to S5 until the overall contour reconstruction of the object to be measured is completed.

优选地,步骤S3中光带图像在二值化处理处理前进行均值滤波。Preferably, in step S3, the light band image is subjected to mean value filtering before binarization processing.

优选地,步骤S6中k的大小与待测物的表面曲率相关。Preferably, the size of k in step S6 is related to the surface curvature of the object to be tested.

本发明的有益效果是:本发明所提供的表面曲率突变光学轮廓测量方法,采用平均灰度算法判断台阶曲率突变的特殊位置,在三维轮廓重建的过程中,能对曲率突变位置作出相应的识别处理,准确的还原出该表面曲率突变物体的轮廓。The beneficial effects of the present invention are: the optical profile measurement method for sudden changes in surface curvature provided by the present invention uses the average gray scale algorithm to judge the special position of sudden changes in step curvature, and can identify the positions of sudden changes in curvature during the process of three-dimensional contour reconstruction processing, and accurately restore the outline of the object with a sudden change in surface curvature.

附图说明Description of drawings

图1是本发明表面曲率突变光学轮廓测量方法流程图;Fig. 1 is a flow chart of the optical profilometry method for sudden changes in surface curvature of the present invention;

图2是本发明表面曲率突变光学轮廓测量装置的示意图;Fig. 2 is the schematic diagram of the surface curvature abrupt change optical profilometry device of the present invention;

图3是本发明被测物台阶的示意图;Fig. 3 is the schematic diagram of the measured object step of the present invention;

图4是本发明第一台阶位置的灰度图像;Fig. 4 is the grayscale image of the first step position of the present invention;

图5是本发明台阶临界位置的灰度图像;Fig. 5 is the gray scale image of step critical position of the present invention;

图6是本发明第二台阶位置的灰度图像;Fig. 6 is the grayscale image of the second step position of the present invention;

图7是本发明被测物台阶平均灰度变化图;Fig. 7 is a step average gray scale change diagram of the measured object in the present invention;

图8是本发明被测物台阶的轮廓重建效果图;Fig. 8 is an effect diagram of contour reconstruction of the steps of the measured object according to the present invention;

附图标记说明:1、固定架;2、CCD摄像机;3、线结构激光器;4、位移平台;5、待测物。Explanation of reference signs: 1. Fixing frame; 2. CCD camera; 3. Line structure laser; 4. Displacement platform; 5. Object to be measured.

具体实施方式detailed description

下面结合附图和具体实施例对本发明做进一步的说明:The present invention will be further described below in conjunction with accompanying drawing and specific embodiment:

如图2所示,本发明的表面曲率突变光学轮廓测量装置,包括固定架1、两台轴对称的CCD摄像机2、三只共面的线结构激光器3、位移平台4和一台计算机。两台CCD摄像机2通过三维调节平台固定到固定架1的两个相邻的支柱上,并通过调节三维调节平台保证两台CCD摄像机的光轴的交点在测量区域的中心。三只线结构激光器3通过二维调节平台固定在固定架1的另外两个支柱及下方位置,使三只线结构激光器3的夹角成120°,以确保每一只激光器到被测物中心轴等距,这样保证照射到物体表面的激光强度相同。As shown in Fig. 2, the optical profile measuring device with sudden change in surface curvature of the present invention includes a fixed frame 1, two axisymmetric CCD cameras 2, three coplanar line-structured lasers 3, a displacement platform 4 and a computer. Two CCD cameras 2 are fixed to two adjacent pillars of the fixed frame 1 through a three-dimensional adjustment platform, and the intersection of the optical axes of the two CCD cameras is ensured at the center of the measurement area by adjusting the three-dimensional adjustment platform. The three line-structure lasers 3 are fixed on the other two pillars and the lower position of the fixed frame 1 through a two-dimensional adjustment platform, so that the angle between the three line-structure lasers 3 is 120° to ensure that each laser reaches the center of the measured object The axes are equidistant, so that the intensity of the laser light on the surface of the object is guaranteed to be the same.

由于CCD摄像机和线结构激光器均固定在固定架上,要实现激光对被测物的纵向扫描来获取系列的轮廓线,如果移动CCD摄像机和线结构激光器来扫描,那么系统的稳定性和可操作性将会下降,并且很难控制扫描速度。为了解决上述问题,本测量系统采用固定CCD摄像机和线结构激光器的位置,移动待测物相对于激光的位置来扫描物体,为方便被测物的移动速度和距离的控制,将位移平台4放置于固定架1中部下方的导轨上并确保固定在位移平台4上的待测物5于测量区域的中心点上,位移平台4通过步进电机进行驱动可以沿导轨移动。Since the CCD camera and the line-structure laser are fixed on the fixed frame, it is necessary to realize the longitudinal scanning of the laser to the measured object to obtain a series of contour lines. If the CCD camera and the line-structure laser are moved to scan, then the stability and operability of the system Performance will drop and it will be difficult to control the scanning speed. In order to solve the above problems, this measurement system adopts the position of fixed CCD camera and line structure laser, and moves the position of the object to be measured relative to the laser to scan the object. In order to facilitate the control of the moving speed and distance of the object to be measured, the displacement platform 4 is placed On the guide rail below the middle part of the fixed frame 1 and ensure that the object to be measured 5 fixed on the displacement platform 4 is at the center point of the measurement area, the displacement platform 4 is driven by a stepping motor to move along the guide rail.

计算机连接CCD摄像机来捕获图像,连接步进电机来驱动控制位移平台移动,连接线结构激光器以控制线结构激光器的亮或灭。The computer is connected to a CCD camera to capture images, a stepper motor is connected to drive and control the movement of the displacement platform, and a line-structure laser is connected to control the light on or off of the line-structure laser.

本发明还提供一种基于上述装置的表面曲率突变光学轮廓测量方法,其流程图如图1所示,调节二维调整平台使线结构激光器的激光面共面,三个线结构激光器发射的激光在待测物表面形成一个待测物横截面轮廓的激光光带,再使用CCD摄像机获取光带图像,将获取的图像传送到计算机中计算被测物的轮廓信息,具体包括以下步骤:The present invention also provides a method for measuring the optical profile of a sudden change in surface curvature based on the above device, the flow chart of which is shown in Fig. Form a laser light band of the cross-sectional profile of the test object on the surface of the test object, and then use the CCD camera to obtain the image of the light band, and transmit the acquired image to the computer to calculate the profile information of the test object, specifically including the following steps:

S1:待测物处于位置j时,利用计算机软件发出脉冲调制信号对激光器进行调制,激光器输出线型光束对待测物进行照明;同时利用计算机软件发出触发信号,触发CCD图像传感器同步采集待测物表面的光带图像,并将采集的图像记为A1,然后关闭激光器,CCD摄像机再次采集相同位置的图像,记为A11S1: When the object under test is at position j, use computer software to send pulse modulation signals to modulate the laser, and the laser outputs a linear beam to illuminate the object under test; at the same time, use computer software to send a trigger signal to trigger the CCD image sensor to synchronously collect the object under test The light band image on the surface, and record the collected image as A 1 , then turn off the laser, and the CCD camera collects the image at the same position again, and record it as A 11 ;

S2:对由锁定成像所得的图像A1、A11进行处理,在此位置获得一幅去除背景光的光带图像IjS2: Process the images A 1 and A 11 obtained by locking imaging, and obtain a light band image I j at this position with the background light removed;

S3:对上述S2中获得的光带图像Ij进行均值滤波和二值化处理;S3: Perform mean value filtering and binarization processing on the light band image I j obtained in the above S2;

S4:对二值化处理后的光带图像Ij进行平均灰度计算,设该图片中有N个像素点,第i个像素点的灰度值为Ni,则光带图像Ij的平均灰度为: S4: Calculate the average gray level of the light band image I j after binarization, assuming that there are N pixels in the picture, and the gray value of the i-th pixel is N i , then the light band image I j The average grayscale is:

S5:将待测物从位置j沿导轨平移距离L到达位置j+1处,重复步骤S1至S4,得到与光带图像I位置连续的光带图像Ij+1,并计算出光带图像Ij+1的平均灰度 S5: Translate the object to be measured from position j along the guide rail to the position j+1, repeat steps S1 to S4, obtain the light band image I j+1 continuous with the position of the light band image I, and calculate the light band image I Average gray level of j+1

S6:将光带图像Ij+1的平均灰度除以光带图像Ij的平均灰度,得到比值进行判断,若比值则光带图像Ij+1为曲率突变图像,该位置为曲率突变位置,在待测物的每个位置,计算机软件均采用三维表面轮廓重建技术进行同步的轮廓重建,若该位置为曲率突变位置,则滤出该位置轮廓线后返回S5;若比值则该位置没有曲率突变,返回S5;直到完成待测物整体的轮廓重建。S6: Divide the average gray level of the light band image I j +1 by the average gray level of the light band image I j to obtain the ratio Make a judgment, if the ratio Then the light band image I j+1 is a curvature mutation image, and this position is the curvature mutation position. At each position of the object to be measured, the computer software uses the three-dimensional surface contour reconstruction technology to perform simultaneous contour reconstruction. If the position is a curvature mutation position, then filter out the contour line of the position and return to S5; if the ratio Then there is no abrupt curvature change at this position, and return to S5 until the overall contour reconstruction of the object to be measured is completed.

本发明的表面曲率突变光学轮廓测量装置经过激光锁定成像技术取得了复杂形面轮廓的光带图像,消除了环境光对光带图像的影响,然而,由于被测物体本身的轮廓不匀称、表面反射率变化剧烈不均匀等原因,获取的光带图像会产生缺损,或光带边缘毛刺等,严重影响了光带图像的质量,所以在一般的图像处理中我们对复杂结构的轮廓光带进行了一些预处理,主要有图像均值滤波算法和图像二值化。The surface curvature abrupt change optical profile measurement device of the present invention obtains the light band image of the complex surface profile through the laser locking imaging technology, and eliminates the influence of ambient light on the light band image. However, due to the uneven profile of the measured object itself, the surface For reasons such as sharp and uneven changes in reflectivity, the acquired light band image will have defects, or light band edge burrs, etc., which seriously affect the quality of the light band image. Therefore, in general image processing, we carry out the Some preprocessing is done, mainly image mean filtering algorithm and image binarization.

在对光带图像进行二值化的同时,我们同步采用了平均灰度的算法从图片中分辨出台阶曲率突变的位置。设计算机采集的一副经过处理的灰度图像中有N个像素点,第i个像素点的灰度值为Ni,则这幅图像的平均灰度为:当计算机采集的位置到达物体曲率突变的位置时,显示出来的平均灰度这一数据将于其他非突变位置有明显的差异,从而轻松的判断出来这些特殊位置。While binarizing the light strip image, we simultaneously use the average gray level algorithm to distinguish the position of the step curvature mutation from the image. Suppose there are N pixels in a processed grayscale image collected by the computer, and the grayscale value of the i-th pixel is N i , then the average grayscale of this image is: When the position collected by the computer reaches the position where the curvature of the object changes suddenly, the displayed average gray level data will be significantly different from other non-mutated positions, so that these special positions can be easily judged.

本实施例的待测物以台阶为例进行说明,如图3所示,该待测物具有四个台阶,从下往上的两个台阶分别为第一台阶和第二台阶,它们之间用过一个垂直面过渡,该垂直面为台阶临界位置。The object to be tested in this embodiment is described by taking steps as an example. As shown in Figure 3, the object to be tested has four steps, and the two steps from bottom to top are respectively the first step and the second step. A vertical plane transition is used, which is the critical position of the step.

图4至图6所示为计算机采集到的并经过预处理的三幅图像,此时都为灰度图像。它们分别为激光扫到被测物的第一台阶位置所得的图像,扫到台阶临界位置(此位置为该物体曲率突变的位置)所得的图像,扫到第二台阶所得的图像。Figures 4 to 6 show the three images collected by the computer and preprocessed, all of which are grayscale images at this time. They are the image obtained by scanning the laser to the first step position of the measured object, the image obtained by scanning to the critical position of the step (this position is the position where the curvature of the object changes suddenly), and the image obtained by scanning to the second step.

经过反复的实验测试,我们得出第一台阶位置和第二台阶位置所测得图片平均灰度在0.650~0.750之前浮动,而对于曲率突变的台阶临界位置,此位置的平均灰度在3.50~3.90之间浮动。可见台阶临界位置的平均灰度要比第一台阶位置和第二台阶位置的平均灰度大的多,前者为后者的五到六倍,从而本实施例中待测物台阶的k=5。After repeated experiments and tests, we found that the average gray level of the images measured at the first step position and the second step position fluctuates before 0.650-0.750, while for the critical position of the step with a sudden curvature change, the average gray level at this position is between 3.50 and 0.750. Floating between 3.90. It can be seen that the average gray level of the critical position of the step is much larger than the average gray level of the first step position and the second step position, and the former is five to six times that of the latter, so k=5 of the step of the object to be tested in this embodiment .

因为测量系统经历的是连续扫描的过程,电脑采集的图像也是具有一定的规律,先采集一部分第一台阶位置的图片,当第一台阶被激光扫描完毕,接着采集到一两张台阶临界位置的图片,来到第二台阶的位置,有采集到大量的第二台阶的图片,灰度值在此过程中的变化如图7所示,可以很直观的发现待测物位于位置坐标6时,激光恰好扫到曲率突变的位置。Because the measurement system undergoes a continuous scanning process, the images collected by the computer also have certain rules. First, a part of the pictures of the position of the first step is collected. After the first step is scanned by the laser, one or two images of the critical position of the step are collected. Picture, when we came to the position of the second step, a large number of pictures of the second step were collected. The change of the gray value during this process is shown in Figure 7. It can be intuitively found that when the object under test is located at position coordinate 6, The laser scans exactly where the curvature changes.

用上述表面曲率突变光学轮廓测量方法对该台阶进行测量的具体过程如下:The specific process of measuring the step with the above-mentioned surface curvature abrupt change optical profilometry method is as follows:

台阶固定在位移平台上处于位置1处,利用计算机软件发出脉冲调制信号对激光器进行调制,激光器输出线型光束对待测物进行照明;同时利用计算机软件发出触发信号,触发CCD图像传感器同步采集待测物表面的光带图像,并将采集的图像记为A1,然后关闭激光器,CCD摄像机再次采集相同位置的图像,记为A11;对由锁定成像所得的图像A1、A11进行处理,在此位置获得一幅去除背景光的光带图像I1;对光带图像I1进行均值滤波和二值化处理后计算出光带图像I1的平均灰度 The steps are fixed on the displacement platform at position 1, and the laser is modulated by sending a pulse modulation signal using computer software, and the laser outputs a linear beam to illuminate the object to be measured; at the same time, the computer software is used to send a trigger signal to trigger the CCD image sensor to synchronously collect the object to be tested The light band image on the surface of the object, and record the collected image as A 1 , then turn off the laser, and the CCD camera collects the image at the same position again, which is recorded as A 11 ; process the images A 1 and A 11 obtained by locking imaging, Obtain a light strip image I 1 with background light removed at this position; perform mean filtering and binarization on the light strip image I 1 to calculate the average gray level of the light strip image I 1

步进电机驱动控制位移平台沿导轨平移距离40mm到达位置2,重复上述步骤获得光带图像I2,计算出光带图像I2的平均灰度小于台阶的k值5,位置2没有曲率突变。The stepper motor drives and controls the displacement platform to move along the guide rail at a distance of 40mm to reach position 2, repeat the above steps to obtain the light band image I 2 , and calculate the average gray level of the light band image I 2 but The k value of 5 is less than the step, and there is no abrupt curvature change at position 2.

步进电机驱动控制位移平台沿导轨平移距离40mm到达位置3,重复上述步骤获得光带图像I3,计算出光带图像I3的平均灰度为0.698;则小于台阶的k值5,位置3也没有曲率突变。The stepper motor drives and controls the displacement platform to move 40mm along the guide rail to reach position 3, repeat the above steps to obtain the light band image I 3 , and calculate the average gray level of the light band image I 3 is 0.698; then The k value of 5 is smaller than that of the step, and there is no abrupt curvature change at position 3.

继续重复上述步骤,依次测得 Continue to repeat the above steps to measure

当台阶移动到位置6时,测得其光带图像I6的平均灰度为3.792,其与光带图像I5的平均灰度的比值则位置6为曲率突变位置,同步轮廓重建滤出该位置轮廓线。When the step moves to position 6, the average gray level of its light band image I 6 is measured is 3.792, which is the same as the average gray level of the light band image I 5 ratio of Then the position 6 is the abrupt change position of the curvature, and the synchronous contour reconstruction filters out the contour line of this position.

重复上述步骤,直到完成待测物整体的轮廓重建。The above steps are repeated until the overall contour reconstruction of the object under test is completed.

如图8所示,为台阶从位置1到位置11的轮廓重建效果图,其中由于曲率突变位置6的轮廓线被滤出,故会留下缝隙,但是由于该缝隙的宽度极小,可以忽略不计。As shown in Figure 8, it is the outline reconstruction effect diagram of the steps from position 1 to position 11, in which the contour line at position 6 with sudden curvature change is filtered out, so a gap will be left, but because the width of the gap is extremely small, it can be ignored Excluding.

本发明的表面曲率突变光学轮廓测量方法,得到满足被测物体目标面形数据,所需的激光光带的灰度图像后,需要对图像进行处理,本发明采用细化处理和和均值滤波消除噪声,之后对图像进行二值化处理,本发明采用平均灰度算法可以判断出不利于轮廓重建的曲率突变的位置。获得轮廓数据后本发明采用最短对角线的轮廓线拼接来重构曲率突变物体的三维表面轮廓,最后借助VTK(visualizationtoolkit)进行显示三维轮廓。The surface curvature abrupt change optical profile measurement method of the present invention needs to process the image after obtaining the grayscale image of the laser light band that satisfies the target surface data of the measured object, and the present invention adopts thinning processing and mean value filtering to eliminate noise, and then binarize the image, and the present invention uses the average gray level algorithm to determine the position of the sudden curvature that is not conducive to contour reconstruction. After obtaining the contour data, the present invention adopts the shortest diagonal contour splicing to reconstruct the three-dimensional surface contour of the object with abrupt curvature, and finally displays the three-dimensional contour by means of VTK (visualization toolkit).

本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。Those skilled in the art will appreciate that the embodiments described here are to help readers understand the principles of the present invention, and it should be understood that the protection scope of the present invention is not limited to such specific statements and embodiments. Those skilled in the art can make various other specific modifications and combinations based on the technical revelations disclosed in the present invention without departing from the essence of the present invention, and these modifications and combinations are still within the protection scope of the present invention.

Claims (3)

1.一种表面曲率突变光学轮廓测量方法,其特征在于,包括以下步骤:1. A surface curvature abrupt change optical profilometry method, is characterized in that, comprises the following steps: S1:待测物处于位置j时,利用计算机软件发出脉冲调制信号对激光器进行调制,激光器输出线型光束对待测物进行照明;同时利用计算机软件发出触发信号,触发CCD图像传感器同步采集待测物表面的光带图像,并将采集的图像记为A1,然后关闭激光器,CCD摄像机再次采集相同位置的图像,记为A11S1: When the object under test is at position j, use computer software to send pulse modulation signals to modulate the laser, and the laser outputs a linear beam to illuminate the object under test; at the same time, use computer software to send a trigger signal to trigger the CCD image sensor to synchronously collect the object under test The light band image on the surface, and record the collected image as A 1 , then turn off the laser, and the CCD camera collects the image at the same position again, and record it as A 11 ; S2:对由锁定成像所得的图像A1、A11进行处理,在此位置获得一幅去除背景光的光带图像IjS2: Process the images A 1 and A 11 obtained by locking imaging, and obtain a light band image I j at this position with the background light removed; S3:对上述S2中获得的光带图像Ij进行二值化处理;S3: Binarize the light band image I j obtained in the above S2; S4:对二值化处理后的光带图像Ij进行平均灰度计算,设该图片中有N个像素点,第i个像素点的灰度值为Ni,则光带图像Ij的平均灰度为: S4: Calculate the average gray level of the light band image I j after binarization, assuming that there are N pixels in the picture, and the gray value of the i-th pixel is N i , then the light band image I j The average grayscale is: S5:将待测物从位置j沿导轨平移距离L到达位置j+1处,重复步骤S1至S4,得到与光带图像I位置连续的光带图像Ij+1,并计算出光带图像Ij+1的平均灰度Navej+1S5: Translate the object to be measured from position j along the guide rail to the position j+1, repeat steps S1 to S4, obtain the light band image I j+1 continuous with the position of the light band image I, and calculate the light band image I The average gray level Nave j+ 1 of j+1 ; S6:将光带图像Ij+1的平均灰度除以光带图像Ij的平均灰度,得到比值进行判断,若比值则光带图像Ij+1为曲率突变图像,该位置为曲率突变位置,在待测物的每个位置,计算机软件均采用三维表面轮廓重建技术进行同步的轮廓重建,若该位置为曲率突变位置,则滤出该位置轮廓线后返回S5;若比值则该位置没有曲率突变,返回S5;直到完成待测物整体的轮廓重建。S6: Divide the average gray level of the light band image I j +1 by the average gray level of the light band image I j to obtain the ratio Make a judgment, if the ratio Then the light band image I j+1 is a curvature mutation image, and this position is the curvature mutation position. At each position of the object to be measured, the computer software uses the three-dimensional surface contour reconstruction technology to perform simultaneous contour reconstruction. If the position is a curvature mutation position, then filter out the contour line of the position and return to S5; if the ratio Then there is no abrupt curvature change at this position, and return to S5 until the overall contour reconstruction of the object to be measured is completed. 2.根据权利要求1所述的表面曲率突变光学轮廓测量方法,其特征在于:所述步骤S3中光带图像在二值化处理处理前进行均值滤波。2. The method for measuring optical profilometry with sudden changes in surface curvature according to claim 1, characterized in that: in the step S3, the light band image is subjected to mean value filtering before binarization processing. 3.根据权利要求1所述的表面曲率突变光学轮廓测量方法,其特征在于:所述步骤S6中k的大小与待测物的表面曲率相关。3. The method for measuring optical profilometry with sudden changes in surface curvature according to claim 1, characterized in that: the size of k in the step S6 is related to the surface curvature of the object to be measured.
CN201610311407.5A 2016-05-12 2016-05-12 Surface curvature is mutated optical profilometry methodology Expired - Fee Related CN105783782B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610311407.5A CN105783782B (en) 2016-05-12 2016-05-12 Surface curvature is mutated optical profilometry methodology

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610311407.5A CN105783782B (en) 2016-05-12 2016-05-12 Surface curvature is mutated optical profilometry methodology

Publications (2)

Publication Number Publication Date
CN105783782A true CN105783782A (en) 2016-07-20
CN105783782B CN105783782B (en) 2018-06-15

Family

ID=56402100

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610311407.5A Expired - Fee Related CN105783782B (en) 2016-05-12 2016-05-12 Surface curvature is mutated optical profilometry methodology

Country Status (1)

Country Link
CN (1) CN105783782B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106645168A (en) * 2016-12-19 2017-05-10 湘潭大学 Detection method for surface concave-convex defect of boom cylinder of crane
CN113029014A (en) * 2019-12-25 2021-06-25 隆基绿能科技股份有限公司 Method and device for detecting diameter of single crystal silicon rod

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6246480B1 (en) * 1999-09-01 2001-06-12 Lucent Technologies Inc. Stepped etalon
CN103080940A (en) * 2010-08-11 2013-05-01 金属板材技术解决方案股份有限公司 Method for generating optimized, constant curvature 2D or 3D roller profile curves and corresponding device
CN103278090A (en) * 2013-05-14 2013-09-04 陕西科技大学 Visual measurement method for volume of irregular object
CN103499585A (en) * 2013-10-22 2014-01-08 常州工学院 Non-continuity lithium battery thin film defect detection method and device based on machine vision
CN104574502A (en) * 2014-12-22 2015-04-29 博迈科海洋工程股份有限公司 Laser cross section feature identification method based on steel structure model
CN104764413A (en) * 2014-12-22 2015-07-08 博迈科海洋工程股份有限公司 Novel method for measuring welding deformation of marine structure deck plate
CN104808688A (en) * 2015-04-28 2015-07-29 武汉大学 Unmanned aerial vehicle curvature continuous adjustable path planning method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6246480B1 (en) * 1999-09-01 2001-06-12 Lucent Technologies Inc. Stepped etalon
CN103080940A (en) * 2010-08-11 2013-05-01 金属板材技术解决方案股份有限公司 Method for generating optimized, constant curvature 2D or 3D roller profile curves and corresponding device
CN103278090A (en) * 2013-05-14 2013-09-04 陕西科技大学 Visual measurement method for volume of irregular object
CN103499585A (en) * 2013-10-22 2014-01-08 常州工学院 Non-continuity lithium battery thin film defect detection method and device based on machine vision
CN104574502A (en) * 2014-12-22 2015-04-29 博迈科海洋工程股份有限公司 Laser cross section feature identification method based on steel structure model
CN104764413A (en) * 2014-12-22 2015-07-08 博迈科海洋工程股份有限公司 Novel method for measuring welding deformation of marine structure deck plate
CN104808688A (en) * 2015-04-28 2015-07-29 武汉大学 Unmanned aerial vehicle curvature continuous adjustable path planning method

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106645168A (en) * 2016-12-19 2017-05-10 湘潭大学 Detection method for surface concave-convex defect of boom cylinder of crane
CN106645168B (en) * 2016-12-19 2019-05-03 湘潭大学 A method for detecting concave and convex defects on the surface of a crane boom
CN113029014A (en) * 2019-12-25 2021-06-25 隆基绿能科技股份有限公司 Method and device for detecting diameter of single crystal silicon rod
CN113029014B (en) * 2019-12-25 2023-03-31 隆基绿能科技股份有限公司 Method and device for detecting diameter of single crystal silicon rod

Also Published As

Publication number Publication date
CN105783782B (en) 2018-06-15

Similar Documents

Publication Publication Date Title
CN114049294B (en) Road surface technical condition detection method and equipment based on three-dimensional contour
CN204730814U (en) A kind of parts passer based on line laser three-dimensional measurement
CN105526882B (en) Track switch Abrasion detecting system and its detection method based on structural light measurement
CN109612390B (en) Large-size workpiece automatic measuring system based on machine vision
US6909514B2 (en) Wheel profile inspection apparatus and method
CN104913737A (en) Component quality checking device based on line laser three-dimensional measurement and detection method of device
CN107764205A (en) High-frequency resistance welding (HFRW) seam center three-dimensional detection device and detection method are scanned based on line-structured light
Tao et al. Weak scratch detection and defect classification methods for a large-aperture optical element
CN101004389A (en) Method for detecting 3D defects on surface of belt material
CN102589516B (en) Dynamic distance measuring system based on binocular line scan cameras
Deng et al. Binocular video-based 3D reconstruction and length quantification of cracks in concrete structures
CN104930985A (en) Binocular vision three-dimensional morphology measurement method based on time and space constraints
CN105203552A (en) 360-degree tread image detecting system and method
CN101639452A (en) 3D detection method for rail surface defects
CN101458072A (en) Three-dimensional contour outline measuring set based on multi sensors and measuring method thereof
CN118602979B (en) Weld joint morphology and surface defect online detection and evaluation method and system
CN102175692A (en) System and method for detecting defects of fabric gray cloth quickly
CN105133471A (en) Linear structured light pavement surface detection system-based pavement depth image production method
CN112284258B (en) System for measuring cable structure size parameters based on machine vision algorithm
CN102798349A (en) Three-dimensional surface extraction method based on equal-gray line search
CN115791801A (en) 3D glass online monitoring platform based on machine vision
CN111102928A (en) Calibration system and method for detecting hard spots of contact net
CN116402792A (en) A 3D Point Cloud-Based Spatial Hole Location Docking Method
TWI490481B (en) On - line Inspection Method for Panel 3D Defects
CN104515479A (en) Medium plate plane shape measuring system and method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180615

Termination date: 20200512