CN101422787A - 基于单步相移法的带钢平坦度测量方法 - Google Patents
基于单步相移法的带钢平坦度测量方法 Download PDFInfo
- Publication number
- CN101422787A CN101422787A CNA2008102391574A CN200810239157A CN101422787A CN 101422787 A CN101422787 A CN 101422787A CN A2008102391574 A CNA2008102391574 A CN A2008102391574A CN 200810239157 A CN200810239157 A CN 200810239157A CN 101422787 A CN101422787 A CN 101422787A
- Authority
- CN
- China
- Prior art keywords
- band steel
- flatness
- steel
- phase
- measurement
- 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
Links
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明将相位测量法引入平坦度检测领域,提供了一种基于单步相移法的带钢平坦度测量方法。首先利用计算机程序编制正弦条纹图像,灰度值按正弦规律变化,将条纹图像通过DLP投影仪倾斜投影到带钢表面。若带钢有浪形,则带钢纵向纤维相对于参考平面会发生偏移和扭曲,带钢表面的条纹也因而产生变形,所以变形条纹图像中包含了所有的板形特征信息。采用CCD摄像头垂直拍摄变形条纹,由条纹图相位与钢板表面高度的对应关系求得钢板上各点的高度值,最后由平坦度计算模型即可计算出平坦度参数。此测量方法属于全场式三维测量方法,能克服带钢运动对测量造成的不良影响。
Description
技术领域
本发明属于轧钢技术领域,特别是提供了一种带钢平坦度在线测量方法,以及具体的实施方式。
背景技术
在带钢轧制生产过程中,带钢的平坦度是重要的控制目标,直接影响带钢的最终质量,故在板带轧机出口设置平坦度检测装置十分重要。传统的非接触式平坦度测量以三角法、光切法和莫尔法为代表,但它们却存在着不容忽视的缺点:激光三角法原理简单,响应速度快,在线数据处理容易实现,但带钢振动和飘摆给带钢纤维长度的测量带来了不确定误差。光切法的优点是包含有一定的带钢三维形状信息,在某种程度上避免了带钢垂直运动的影响,简化了板形计算模型,使在线测量成为可能。但是激光扫描需要大功率的激光器作光源,造价昂贵,安装调试不便,并且强激光束对人眼有危害。莫尔法可以实时地测量带钢形状,通过适当的数据处理可以克服带钢跳动和摇摆对板形测量的影响,但在实际测量中需要一个较大的矩形耐热格栅,格栅的加工、耐热、变形以及靠近被测热轧带钢的安装问题都妨碍了系统的可靠性,而且成本极高。
相位测量法以测量投影到物体上的变形光栅像的相位为基础,通过相位与高度的映射关系得到被测物体的三维轮廓,具有较高的灵敏度。相位测量法的基本思想是将正弦光栅条纹投射到被测物表面,从另一角度可以观察到由于受物体高度的调制而变形的条纹,这种变形可解释为相位和振幅均被调制的空间载波信号。采集变形条纹并对其进行解调,恢复出代表物体高度的相位信息,再经过相位展开和几何计算就可以获得被测物体形面的三维几何信息。因为相位测量法是直接采集被测物体表面上变形光栅的图像,充分利用了物体表面的信息,变形光栅图像中的每个像素都对应着确定的相位值和高度值等数据信息,因此具有测量速度快、精度高等特点。
求解相位是通过相位测量法获取物体三维轮廓的关键,目前常用的求解相位方法有三步相移、四步相移、五步相移和单步相移等,其中三步相移、四步相移、五步相移法等都需要两幅以上的图像,数据采集量大,采集时间长,对那些具有速度要求的场合不可能采集到足够的数据,不适合于运动物体的测量。
单步相移算法是将投影光栅线移动0°,也就是不作相移(即单步),但提取相位时利用传统的相移法提取技术。该方法只需采集一幅条纹图,对于该条纹图中的任意一个像素点,把与它相邻的像素的灰度值看成是该点的相位分别移动某一固定值的结果,也就是将一幅图拆分成几幅相移图,然后利用传统的三步相移、四步相移、五步相移等相位提取技术进行相位的求解。
发明内容
本发明将相位测量法引入带钢平坦度的测量,提供了一种基于单步相移法的非接触式平坦度检测方法,它属于全场式三维测量方法,能有效克服带钢振动和飘摆对测量造成的不良影响,与其它测量平坦度方法相比,能更准确地测量出带钢平坦度。
首先利用计算机程序编制正弦条纹图像(灰度值按正弦规律变化),将条纹图像通过DLP投影仪倾斜投影到带钢表面。若带钢有浪形,则带钢纵向纤维相对于参考平面会发生偏移和扭曲,带钢表面的条纹也因而产生变形,所以变形条纹图像中包含了所有的板形特征信息。由垂直安装在带钢中心侧正上方的面阵CCD摄像机拍摄变形条纹图像,并将图像保存于计算机中进行相位提取及去包裹处理。然后由变形条纹的相位与钢板表面高度的对应关系求得钢板上各点的高度值,最后由平坦度计算模型即可计算出板形参数。
一种利用连续3个像素计算相位的方法,其计算原理如下。
利用软件生成周期不同的正弦条纹,这样在斜向投影时可通过调节条纹周期和CCD成像透镜使得CCD拍摄到的未变形的条纹周期大约为8个像素的大小。因此,在垂直于条纹方向的直线上相邻两像素可认为有的相移量。条纹的图像光强分布可表示为:
式(3)中,i,j分别为水平和垂直方向像素数;B(i,j)为背景光强,C(i,j)为被测物表面的光强;为与像素位置和物体高度分布相关的变形条纹的相位。
假设B(i,j)和C(i,j)在6像素范围内基本不变。则由式(1)可推导出提取相位的公式:
式(2)提取相位的3个灰度数据取自连续6个像素范围中的3个间隔位置。为了提高空间分辨率,本发明推导出利用连续3个像素计算相位的公式:
式(5)是假设B(i,j)和C(i,j)在3像素范围内基本不变的情况下得到的,此假设比6像素法更为合理,空间分辨率比式(3)高出一倍。利用式(5)计算出的相位包括物体相位、交叉光轴相位和部分频率相位,其中后两项并不表征物体表面的高度起伏信息,对于给定的光学系统而言,它们是恒定的常量,因此,可以先测量一个标准平面,将其作为系统的零分布,存储其相位作为系统基本数据在获取被测带钢的相位数据后,将其减去基本数据即可得到单纯的被测带钢的相位信息。
一种针对相位测量法的平坦度计算模型,其计算原理如下:
带钢形成浪形的实质是带钢横向各点沿纵向产生不均匀延伸,所以平坦度可用带钢纵向纤维的相对延伸来表示:
式(4)中:
ε(χj):χj点的相对延伸
Lj:沿带钢宽度不同位置的纵向纤维长度
L0:平直纤维的水平长度
考虑到实际测量中所得高度值z是离散采样值z0,z1,z2,z3L zn,则被测纤维长度可按下式计算:
得到纤维长度后,再用式(4)计算带钢纤维的相对延伸,作为平坦度评价参数。
由式(5)可知,计算过程无需考虑时间因素和带钢瞬时速度值,只与测量点的三维坐标值有关,因而可以有效避免由带钢速度测量偏差以及带钢运动所引起的测量误差。
本发明的有益效果是:摄像机一次拍摄一组条纹图像,纤维条上的测量点的数据都是在同一时刻测得,不含有时间和带钢速度的因素,避免了带钢振动和飘摆以及带钢速度对测量的干扰。由于投影正弦条纹的级数多,测量点数量远远多于其他平坦度检测方法,故在测量精度和灵敏度方面都有明显提高,可以实现全场测量,具有精度高、响应速度快、结构简单等优点。
附图说明
下面结合附图对本发明进行进一步说明。
图1为基于单步相移法的带钢平坦度测量方法结构图。图1中,1.DLP(Digital LightProjector)投影仪,2.面阵CCD摄像机。
具体实施方式
在图1中,将DLP(1)倾斜安装于带钢中心侧上方,投影的角度和范围可根据实际情况进行适当调整。通常情况下,DLP(1)光轴及投影光切面的入射角度在20°~50°范围内较为适宜。面阵CCD摄像机(2)的安装高度与DLP(1)相同,光轴垂直于带钢表面,在面阵CCD摄像机(2)前安装滤光片,以滤除杂质光源(如热轧板带的近红外光)。
DLP投影仪将软件编制的正弦光栅条纹(灰度值按正弦规律分布)图像倾斜投射到带钢表面上,由垂直安装在带钢中心侧正上方的面阵CCD摄像机拍摄变形条纹图像,并将图像保存于计算机中。根据式(5)进行相位提取及去包裹处理,从而得到包含带钢三维轮廓信息的相位差值。
为了通过相位差必须进一步转化为被测带钢的实际高度值,从而计算平坦度,这就需要对系统进行标定。为计算方便,本发明采用参考平面作为世界坐标系的XW-YW平面,摄像机光轴作为世界坐标系的ZW轴。被测带钢高度与相位差满足关系式:
由式(6)可以看出,只要确定出a(x,y)、b(x,y)这两个参数,就能根据相位差得到带钢任意位置的实际高度Z值,因此只要给出两组不同的高度h及其对应的相位差即可求出这两个参数。具体标定时,利用相对于参考平面平行移动距离为h1、h2的两个标定平面来获得标定数据,从而求得这两个参数,建立起高度和相位差的映射关系。
通过传统的摄像机标定可以求得摄像机的投影矩阵 这样就可以得到关于XW,YW,ZW的两个线性方程组:
通过已经可以得到的任意点Z值,带入方程组(7)即可求出该点对应的实际坐标(X,Y)。根据式(4)和式(5),即可得到带钢平坦度的评价指标。
Claims (1)
1、一种基于单步相移法的带钢平坦度测量方法,其特征在于利用计算机产生正弦条纹图像,灰度值按正弦规律变化,安装于带钢中心侧上方的DLP投影仪将条纹图像倾斜投影到带钢表面,由垂直安装在带钢中心侧正上方的面阵CCD摄像机拍摄变形条纹图像,并将图像保存于计算机中进行相位提取,然后由变形条纹的相位与钢板表面高度的对应关系求得钢板上各点的高度值,最后由平坦度计算模型即可计算出板形参数,所述的平坦度计算模型计算原理如下:
带钢形成浪形的实质是带钢横向各点沿纵向产生不均匀延伸,所以平坦度可用带钢纵向纤维的相对延伸来表示:
上式中:
ε(χj):χj点的相对延伸
Lj:沿带钢宽度不同位置的纵向纤维长度
L0:平直纤维的水平长度
考虑到实际测量中所得高度值z是离散采样值z0,z1,z2,z3L zn,则被测纤维长度
得到纤维长度后,再用上式计算带钢纤维的相对延伸,作为平坦度评价参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102391574A CN101422787B (zh) | 2008-12-10 | 2008-12-10 | 基于单步相移法的带钢平坦度测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102391574A CN101422787B (zh) | 2008-12-10 | 2008-12-10 | 基于单步相移法的带钢平坦度测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101422787A true CN101422787A (zh) | 2009-05-06 |
CN101422787B CN101422787B (zh) | 2011-04-20 |
Family
ID=40613796
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008102391574A Expired - Fee Related CN101422787B (zh) | 2008-12-10 | 2008-12-10 | 基于单步相移法的带钢平坦度测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101422787B (zh) |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102169095A (zh) * | 2009-09-24 | 2011-08-31 | 凯德易株式会社 | 检查系统以及检查方法 |
CN102261896A (zh) * | 2011-04-19 | 2011-11-30 | 长春东瑞科技发展有限公司 | 一种基于相位测量的物体三维形貌测量方法及系统 |
CN102319743A (zh) * | 2011-05-24 | 2012-01-18 | 重庆大学 | 带钢跑偏与漂浮量激光扫描检测方法及纠偏系统 |
CN102620663A (zh) * | 2012-03-21 | 2012-08-01 | 莱芜钢铁集团有限公司 | 热轧板带坯料在线测长装置及其使用方法 |
CN102667400A (zh) * | 2009-10-19 | 2012-09-12 | 住友金属工业株式会社 | 板材的平坦度测量方法以及使用该方法的钢板的制造方法 |
CN104807406A (zh) * | 2014-01-27 | 2015-07-29 | 康耐视公司 | 确定物体上3d表面特征和不平度的系统和方法 |
CN105178149A (zh) * | 2015-08-25 | 2015-12-23 | 吉林大学 | 车载路面裂纹机器视觉检测系统的投影式标定系统 |
CN105229413A (zh) * | 2013-05-14 | 2016-01-06 | 新日铁住金株式会社 | 板材的平面度测定方法、板材的平面度测定装置以及钢板的制造方法 |
CN105598186A (zh) * | 2015-12-30 | 2016-05-25 | 北京首钢自动化信息技术有限公司 | 一种基于带钢横向平坦度分布的平坦度评价系统及其方法 |
CN107282648A (zh) * | 2017-06-21 | 2017-10-24 | 北京科技大学 | 一种热轧带钢全幅宽平坦度的控制方法 |
CN107466356A (zh) * | 2016-04-06 | 2017-12-12 | 4维传感器有限公司 | 测量方法、测量装置、测量程序以及记录了测量程序的计算机可读取记录介质 |
TWI619560B (zh) * | 2017-06-01 | 2018-04-01 | 中國鋼鐵股份有限公司 | 鋼板量測系統及其方法 |
CN108053392A (zh) * | 2017-12-05 | 2018-05-18 | 南京航空航天大学 | 结合深度测量的车轮踏面擦伤视觉检测方法 |
CN108398088A (zh) * | 2017-02-05 | 2018-08-14 | 鞍钢股份有限公司 | 一种带钢宽度测量中的振动误差校正方法 |
CN108474644A (zh) * | 2015-08-10 | 2018-08-31 | 慧咨环球有限公司 | 体积估计方法、设备和系统 |
CN110503048A (zh) * | 2019-08-26 | 2019-11-26 | 中铁电气化局集团有限公司 | 刚性接触网悬挂装置的识别系统和方法 |
CN110864650A (zh) * | 2019-11-25 | 2020-03-06 | 天津大学 | 基于条纹投影的平面度测量方法 |
CN111521129A (zh) * | 2020-04-20 | 2020-08-11 | 北京科技大学 | 基于机器视觉的板坯翘曲检测装置及方法 |
CN112595261A (zh) * | 2020-12-01 | 2021-04-02 | 哈尔滨石油学院 | 一种基于结构光的锈斑检测系统 |
CN113066164A (zh) * | 2021-03-19 | 2021-07-02 | 河南工业大学 | 一种基于相移轮廓术的单条纹运动物体三维重建方法 |
CN113310438A (zh) * | 2020-02-26 | 2021-08-27 | 广东博智林机器人有限公司 | 平整度测量方法、装置、计算机设备和存储介质 |
CN115798625A (zh) * | 2023-01-17 | 2023-03-14 | 中国科学技术大学 | 一种耦合单步化学机理的反应流数值求解方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2695049B1 (fr) * | 1992-09-03 | 1994-11-10 | Lorraine Laminage | Procédé et dispositif d'évaluation de la planéité d'une bande et ligne de fabrication d'une bande en continu comportant un dispositif d'évaluation de la planéité de la bande. |
FR2725512B1 (fr) * | 1994-10-07 | 1996-10-31 | Lorraine Laminage | Dispositif d'evaluation de la planeite d'un produit en defilement |
DE19758466B4 (de) * | 1997-03-11 | 2007-10-04 | Betriebsforschungsinstitut VDEh - Institut für angewandte Forschung GmbH | Planheits-Regelungssystem für Metallband |
DE10205132A1 (de) * | 2002-02-07 | 2003-08-28 | Bfi Vdeh Inst Angewandte Forschung Gmbh | Verfahren und Vorrichtung zum optischen Messen der Oberflächenform und zur optischen Oberflächeninspektion von bewegten Bändern in Walz- und Weiterbearbeitungsanlagen |
-
2008
- 2008-12-10 CN CN2008102391574A patent/CN101422787B/zh not_active Expired - Fee Related
Cited By (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102169095A (zh) * | 2009-09-24 | 2011-08-31 | 凯德易株式会社 | 检查系统以及检查方法 |
CN102169095B (zh) * | 2009-09-24 | 2014-09-03 | 凯德易株式会社 | 检查系统以及检查方法 |
CN102667400B (zh) * | 2009-10-19 | 2014-11-05 | 新日铁住金株式会社 | 板材的平坦度测量方法以及使用该方法的钢板的制造方法 |
CN102667400A (zh) * | 2009-10-19 | 2012-09-12 | 住友金属工业株式会社 | 板材的平坦度测量方法以及使用该方法的钢板的制造方法 |
CN102261896A (zh) * | 2011-04-19 | 2011-11-30 | 长春东瑞科技发展有限公司 | 一种基于相位测量的物体三维形貌测量方法及系统 |
CN102319743A (zh) * | 2011-05-24 | 2012-01-18 | 重庆大学 | 带钢跑偏与漂浮量激光扫描检测方法及纠偏系统 |
CN102319743B (zh) * | 2011-05-24 | 2013-07-10 | 重庆大学 | 带钢跑偏与漂浮量激光扫描检测方法及纠偏系统 |
CN102620663A (zh) * | 2012-03-21 | 2012-08-01 | 莱芜钢铁集团有限公司 | 热轧板带坯料在线测长装置及其使用方法 |
CN105229413A (zh) * | 2013-05-14 | 2016-01-06 | 新日铁住金株式会社 | 板材的平面度测定方法、板材的平面度测定装置以及钢板的制造方法 |
CN105229413B (zh) * | 2013-05-14 | 2017-12-01 | 新日铁住金株式会社 | 板材的平面度测定方法、板材的平面度测定装置以及钢板的制造方法 |
CN104807406A (zh) * | 2014-01-27 | 2015-07-29 | 康耐视公司 | 确定物体上3d表面特征和不平度的系统和方法 |
CN108474644A (zh) * | 2015-08-10 | 2018-08-31 | 慧咨环球有限公司 | 体积估计方法、设备和系统 |
CN105178149A (zh) * | 2015-08-25 | 2015-12-23 | 吉林大学 | 车载路面裂纹机器视觉检测系统的投影式标定系统 |
CN105598186A (zh) * | 2015-12-30 | 2016-05-25 | 北京首钢自动化信息技术有限公司 | 一种基于带钢横向平坦度分布的平坦度评价系统及其方法 |
CN107466356A (zh) * | 2016-04-06 | 2017-12-12 | 4维传感器有限公司 | 测量方法、测量装置、测量程序以及记录了测量程序的计算机可读取记录介质 |
CN108398088A (zh) * | 2017-02-05 | 2018-08-14 | 鞍钢股份有限公司 | 一种带钢宽度测量中的振动误差校正方法 |
TWI619560B (zh) * | 2017-06-01 | 2018-04-01 | 中國鋼鐵股份有限公司 | 鋼板量測系統及其方法 |
CN107282648A (zh) * | 2017-06-21 | 2017-10-24 | 北京科技大学 | 一种热轧带钢全幅宽平坦度的控制方法 |
CN107282648B (zh) * | 2017-06-21 | 2019-01-22 | 北京科技大学 | 一种热轧带钢全幅宽平坦度的控制方法 |
CN108053392A (zh) * | 2017-12-05 | 2018-05-18 | 南京航空航天大学 | 结合深度测量的车轮踏面擦伤视觉检测方法 |
CN110503048A (zh) * | 2019-08-26 | 2019-11-26 | 中铁电气化局集团有限公司 | 刚性接触网悬挂装置的识别系统和方法 |
CN110864650A (zh) * | 2019-11-25 | 2020-03-06 | 天津大学 | 基于条纹投影的平面度测量方法 |
CN113310438A (zh) * | 2020-02-26 | 2021-08-27 | 广东博智林机器人有限公司 | 平整度测量方法、装置、计算机设备和存储介质 |
CN111521129A (zh) * | 2020-04-20 | 2020-08-11 | 北京科技大学 | 基于机器视觉的板坯翘曲检测装置及方法 |
CN111521129B (zh) * | 2020-04-20 | 2021-08-31 | 北京科技大学 | 基于机器视觉的板坯翘曲检测装置及方法 |
CN112595261A (zh) * | 2020-12-01 | 2021-04-02 | 哈尔滨石油学院 | 一种基于结构光的锈斑检测系统 |
CN113066164A (zh) * | 2021-03-19 | 2021-07-02 | 河南工业大学 | 一种基于相移轮廓术的单条纹运动物体三维重建方法 |
CN115798625A (zh) * | 2023-01-17 | 2023-03-14 | 中国科学技术大学 | 一种耦合单步化学机理的反应流数值求解方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101422787B (zh) | 2011-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101422787B (zh) | 基于单步相移法的带钢平坦度测量方法 | |
CN102620685B (zh) | 一种基于史托克维尔变换的改进窗口傅里叶三维测量法 | |
CN100489446C (zh) | 基于相位法的三维轮廓测量方法 | |
CN103383360B (zh) | 一种薄带连铸坯表面缺陷正弦光栅相移检测装置及检测方法 | |
CN103575227B (zh) | 一种基于数字散斑的视觉引伸计实现方法 | |
CN203231736U (zh) | 一种基于双目视觉的镜面物体测量装置 | |
Su et al. | Dynamic 3-D shape measurement method based on FTP | |
CN101986098B (zh) | 基于三色光投影的傅里叶变换三维测量法 | |
CN102032877B (zh) | 基于小波变换的三维测量方法 | |
CN101901501B (zh) | 一种生成激光彩色云图的方法 | |
CN106978774B (zh) | 一种路面坑槽自动检测方法 | |
CN104160241A (zh) | 利用高维亮度信息的条纹图像的相位分布分析方法、装置及其程序 | |
CN110514143A (zh) | 一种基于反射镜的条纹投影系统标定方法 | |
CN104111036A (zh) | 一种基于双目视觉的镜面物体测量装置及方法 | |
CN103776390A (zh) | 三维自然纹理数据扫描机及多视场数据拼接方法 | |
CN103267496B (zh) | 一种基于小波变换的改进窗口傅里叶三维测量法 | |
CN104061879A (zh) | 一种连续扫描的结构光三维面形垂直测量方法 | |
CN109506592A (zh) | 基于条纹光流的物体三维面形测量方法和装置 | |
CN109751971A (zh) | 高精度短相干光三维形貌快速测量算法 | |
CN1266452C (zh) | 复合编码多分辨三维数字成像方法 | |
CN111765851B (zh) | 基于双屏透射显示的镜面物体三维形貌测量方法及装置 | |
CN102589516A (zh) | 一种基于双目线扫描摄像机的动态距离测量系统 | |
CN109631798A (zh) | 一种基于π相移方法的三维面形垂直测量方法 | |
CN106017358A (zh) | 一种基于预校正光栅投影的三维面形测量方法 | |
CN105043301A (zh) | 一种用于三维测量的光栅条纹相位求解方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110420 Termination date: 20131210 |