发明名称 一种冷轧带钢边部在线快速定位方法
摘要 一种冷轧带钢边部在线快速定位方法,将相机采集的带钢实时图像处理后获得带钢同步灰度图像;对灰度图像进行纵向中值滤波,图像上每个像素点的灰度值等于其纵向临近点的灰度中间值;对图像分层,从左右两方向分别确定每层的左右波动最大块;进行图像纵向压缩,分别形成左右压缩图像序列;再分别生成左右边部波动序列;然后,确定每一波动序列的边缘波动阈值,从左右两方向分别搜索波动序列中第一个大于对应边缘波动阈值的波动序列值,确定本层左右边缘位置及跃阶值;最后确定左右边部统计边界。本发明工作可靠,算法简单,运行速度快,处理效果好,不受带钢纹理、色差、缺陷及光照影响,适用于工业生产和各种计算设备应用。
申请公布号 CN105414206A 申请公布日期 2016.03.23
申请号 CN201410483796.0 申请日期 2014.09.19
申请人 鞍钢股份有限公司 发明人 宋宝宇;王军生;王靖震;杨东晓;高冰;李连成;费静;王奎越;柴明亮;赵耕
分类号 B21B38/00(2006.01)I;B21B39/14(2006.01)I 主分类号 B21B38/00(2006.01)I
代理机构 鞍山华惠专利事务所 21213 代理人 赵长芳
主权项 一种冷轧带钢边部在线快速定位方法,其特征在于,定位过程分为左边部的定位和右边部的定位,具体方法和步骤为:(1)边部图像获取:通过相机采集带钢实时同步图像,然后对采集到的图像进行转灰度处理,获得带钢同步灰度图像;(2)图像预处理:对获取的带钢同步灰度图像首先进行纵向中值滤波,即图像上每个像素点的灰度值等于其纵向临近点的灰度中间值;(3)搜索波动最大块:对于左首层波动最大块的搜索,从图像左上部开始,横向M列像素点,纵向N行像素点,以M×N的大小,从左向右分块检验,根据块内灰度初检情况确定左首层波动最大块;对于右首层波动最大块则从图像右上部开始,以M×N的大小,从右向左分块检验,使用同样的方法确定右首层波动最大块;首层以下的每一层最大块纵向范围仍为N行像素点,横向范围则等同于上一层最大块的横向范围;(4)进行图像纵向压缩:对于左边部定位,在步骤(3)确定的左层波动最大块向左右扩展为3M×N区域,左右扩展长度相同,计算每列的纵向平均值,形成3M×1大小的左压缩图像序列FY<sub>L</sub>(x,y),其中x=0,1,2……,y为层数;对于右边部定位,方法相同,形成右压缩图像序列FY<sub>R</sub>(x,y);(5)生成波动序列:对于左边部定位,波动序列FB<sub>L</sub>(x,y)的计算方法为:FB<sub>L</sub>(i,y)=|FY<sub>L</sub>(i,y)‑FY<sub>L</sub>(i‑1,y)|  (i&gt;0)对于右边部定位,波动序列FB<sub>R</sub>(x)的计算方法为:FB<sub>R</sub>(i,y)=|FY<sub>R</sub>(i,y)‑FY<sub>R</sub>(i+1,y)|  (i&lt;压缩图像序列长度)(6)确定层边缘位置:对于左边部定位,首先计算波动序列FB<sub>L</sub>(x,y)中的边缘波动阈值KB<sub>L</sub>(y):KB<sub>L</sub>(y)=AVG<sub>L</sub>(y)+SDV<sub>L</sub>(y)其中,AVG<sub>L</sub>(y)为波动序列FB<sub>L</sub>(x,y)的平均值,SDV<sub>L</sub>(y)为波动序列FB<sub>L</sub>(x,y)的标准差;然后,从左至右搜索波动序列FB<sub>L</sub>(x,y)中第一个大于KB<sub>L</sub>(y)的波动序列值,该值对应的原图像横向像素点位置(x)即为本层左边缘位置LB<sub>L</sub>(Y),该位置所对应的波动序列值记为本层左边缘跃阶值FS<sub>L</sub>(y);对于右边部定位,首先计算波动序列FB<sub>R</sub>(x,y)中的边缘波动阈值KB<sub>R</sub>(y):KB<sub>R</sub>(y)=AVG<sub>R</sub>(y)+SDV<sub>R</sub>(y)其中,AVG<sub>R</sub>(y)为波动序列FB<sub>R</sub>(x,y)的平均值,SDV<sub>R</sub>(y)为波动序列FB<sub>R</sub>(x,y)的标准差;然后,从右至左搜索波动序列FB<sub>R</sub>(x,y)中第一个大于KB<sub>R</sub>(y)的波动序列值,该值对应的原图像横向像素点位置(x)即为本层右边缘位置LB<sub>R</sub>(Y),该位置所对应的波动序列值记为本层右边缘跃阶值FS<sub>R</sub>(y);(7)确定统计边界:对于计算多层图像统计边界的情况,左边部统计边界TGB<sub>L</sub>的计算方法为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>TGB</mi><mi>L</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>Y</mi></munderover><msub><mi>LB</mi><mi>L</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&times;</mo><mfrac><mrow><msub><mi>FS</mi><mi>L</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow><msub><mi>TS</mi><mi>L</mi></msub></mfrac></mrow>]]></math><img file="FDA0000573057970000021.GIF" wi="805" he="148" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>TS</mi><mi>L</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>Y</mi></munderover><msub><mi>FS</mi><mi>L</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000573057970000022.GIF" wi="477" he="141" /></maths>其中,y为层数;右边部统计边界TGB<sub>R</sub>的计算方法为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>TGB</mi><mi>R</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>Y</mi></munderover><msub><mi>LB</mi><mi>R</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&times;</mo><mfrac><mrow><msub><mi>FS</mi><mi>R</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow><msub><mi>TS</mi><mi>R</mi></msub></mfrac></mrow>]]></math><img file="FDA0000573057970000023.GIF" wi="811" he="148" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>TS</mi><mi>R</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>Y</mi></munderover><msub><mi>FS</mi><mi>R</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000573057970000024.GIF" wi="481" he="131" /></maths>所述步骤(3)中,根据块内灰度初检情况确定左首层波动最大块中灰度初检情况判定方法为如果该块首次满足如下检验条件,则确定为左首层波动最大块:检验条件:N<sub>M</sub>&gt;T<sub>K</sub><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>N</mi><mi>M</mi></msub><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mo>(</mo><mo>|</mo><mi>F</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mi>F</mi><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mi>N</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>|</mo><mo>></mo><msub><mi>C</mi><mi>K</mi></msub><mo>)</mo></mrow></munder><mn>1</mn></mrow>]]></math><img file="FDA0000573057970000025.GIF" wi="699" he="144" /></maths>其中,C<sub>K</sub>为带钢基板与非带钢部分的灰度波动预设最小值,T<sub>K</sub>为判断常数。
地址 114021 辽宁省鞍山市铁西区鞍钢厂区内