发明名称 一种塔式起重机结构的非接触振动检测方法
摘要 本发明涉及一种塔式起重机结构的非接触振动检测方法,属于检测技术领域。本方法为非接触式测量,通过配置多个摄像头实现多点和在线测试,无需在被测对象上安装传感器,没有任何负载效应,容易实施。图像信息直接数字化,构成全自动全数字化振动测试处理、记录系统,信息量比传统方法大。其中使用ROI(Regions of Interest)技术,缩小图像跟踪区域和计算量大小,并使用SURF方法加速特征点的提取速度,使用拉伊达准则剔除“离群点”改进SURF方法的精度;提高了现有技术中起重机安全监控系统的可靠性。
申请公布号 CN103033253B 申请公布日期 2015.07.15
申请号 CN201210524875.2 申请日期 2012.12.10
申请人 太原科技大学 发明人 姚峰林;赵婕;宁少慧;高崇仁
分类号 G01H9/00(2006.01)I 主分类号 G01H9/00(2006.01)I
代理机构 代理人
主权项 一种塔式起重机结构的非接触振动检测方法,其特征在于:包括以下步骤:步骤1,采用工业级CCD摄像头对塔式起重机结构进行拍摄,得到一幅目标图像和一路塔机运行过程中的视频,并对视频进行分解,得到多幅连续的视频分解图像;步骤2,对目标图像和视频分解图像分别使用ROI方法获取感兴趣区域;并将目标图像的感兴趣区域作为对比图像,将视频分解图像的感兴趣区域作为图像跟踪区域;步骤3,根据步骤2得到的图像跟踪区域,在步骤1得到的每一幅视频分解图像中,选取对应区域的图像作为待匹配图像;步骤4,对步骤2得到的对比图像以及步骤3得到的多幅待匹配图像,分别进行Hessian变换,并从中提取兴趣点;步骤5,对步骤4得到的两组兴趣点分别进行Harr小波变换,并提取特征点,得到两个特征点集;步骤6,对步骤5得到的两个特征点集进行粗匹配,得到各个待匹配图像与对比图像之间的匹配点集;粗匹配采用最邻近向量匹配方法,具体实现过程如下:采用特征点的特征向量的欧式距离作为相似性测度,设N<sub>1</sub>和N<sub>2</sub>分别是对比图像P和任意一幅待匹配图像Q的特征点集合,对于N<sub>1</sub>中的第m个特征点p<sub>m</sub>,在N<sub>2</sub>中与p<sub>m</sub>的特征向量的欧氏距离最小的特征点为q<sub>mj</sub>和q′<sub>mj</sub>,对应的欧氏距离分别为d<sub>mj</sub>和d′<sub>mj</sub>,如果d<sub>mj</sub>/d′<sub>mj</sub>≤T,T为阈值,则认为p<sub>m</sub>和q<sub>mj</sub>为对应匹配点;步骤7,采用拉伊达准则,对步骤6得到的多个匹配点集分别进行精确匹配,剔除“离群”点;对其中一个匹配点集的精确匹配方法为:求得该匹配点集中所有粗匹配点的欧氏距离平均值<img file="FSB0000136379800000011.GIF" wi="53" he="38" />令常数<img file="FSB0000136379800000012.GIF" wi="140" he="48" />设y<sub>i</sub>=x<sub>i</sub>‑a<sub>0</sub>,其中x<sub>i</sub>为第i对粗匹配点的欧氏距离;则x<sub>i</sub>=y<sub>i</sub>+a<sub>0</sub>;得到用常数a<sub>0</sub>以及所有粗匹配对的欧氏距离与a<sub>0</sub>的差表示的平均值表达式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>x</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>+</mo><msub><mi>a</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>a</mi><mn>0</mn></msub><mo>+</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>y</mi><mi>i</mi></msub><mo>;</mo></mrow>]]></math><img file="FSB0000136379800000013.GIF" wi="798" he="177" /></maths>带入到贝赛尔公式<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&sigma;</mi><mo>=</mo><msqrt><mfrac><mn>1</mn><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow>]]></math><img file="FSB0000136379800000021.GIF" wi="423" he="135" /></maths>得到标准差为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>&sigma;</mi><mo>=</mo><msqrt><mfrac><mn>1</mn><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>a</mi><mn>0</mn></msub><mo>-</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>y</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>=</mo><msqrt><mfrac><mn>1</mn><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>y</mi><mi>i</mi><mn>2</mn></msubsup><mo>-</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><msup><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>y</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo></msqrt></mrow>]]></math><img file="FSB0000136379800000022.GIF" wi="1167" he="140" /></maths>计算欧氏距离值x<sub>i</sub>的偏差<img file="FSB0000136379800000023.GIF" wi="208" he="49" />设定一个动态参数k=k<sub>0</sub>+k<sub>w</sub>,k<sub>0</sub>为一个恒定的初始值,k<sub>w</sub>为k的变化梯度;滤波过程中实时调整k<sub>w</sub>,直到滤波处理后数据的有效匹配点个数小于滤波前个数且大于滤波前个数的3/4;k<sub>0</sub>=1.3,k∈[1.3,5];若|v<sub>i</sub>|>kσ,则判定此匹配点为离群点,剔除其对应的匹配对,得到精确匹配点集;步骤8,对步骤7得到的多个精确匹配点集进行特征点亚像素内插;步骤9,用最近领域法,在步骤8得到的亚像素内插点找到最匹配的亚像素内插点;步骤10,将待匹配图像中亚像素内插点与目标图像中的亚像素内插点相减,得到图像跟踪区域内塔式起重机感兴趣区域的振动曲线。
地址 030024 山西省太原市万柏林区瓦流路66号