发明名称 一种地铁隧道裂缝宽度检测及校正算法
摘要 本发明涉及轨道交通中地铁隧道裂缝宽度的检测和估计。在地铁隧道表面裂缝检测过程中,由于地铁隧道内环境的复杂性,如灰尘掩盖、阴影影响等,对裂缝宽度的检测产生干扰,以及图像处理过程中对裂缝宽度测量的不准确的限制,简单通过一次图像处理得到的裂缝宽度数据并不能完全准确代表裂缝真实宽度,可能存在较大的误差,不能满足安全性要求。因此,本算法提出了一种有效的地铁隧道裂缝宽度检测及校正算法,本发明算法流程分为三个阶段:数据采集阶段、数据处理阶段、校正阶段。本发明对图像处理得到的结果和历史数据进行融合,提高了裂缝宽度预测的真实性。
申请公布号 CN103486971B 申请公布日期 2016.07.13
申请号 CN201310353919.4 申请日期 2013.08.14
申请人 北京交通大学 发明人 张振江;刘云;付俊松;亓大鹏;吴晓阳;赵阳阳;付慧蕾
分类号 G01B11/02(2006.01)I;G06F19/00(2011.01)I 主分类号 G01B11/02(2006.01)I
代理机构 北京市商泰律师事务所 11255 代理人 毛燕生
主权项 一种地铁隧道裂缝宽度检测及校正算法,其特征在于,包括依次进行的三个阶段:数据采集阶段、数据处理阶段、数据校正阶段;通过挖掘历史数据变化规律来估计裂缝最佳值,并且进一步通过一种校验机制来对估计结果进行修正,以确保对裂缝宽度估计的精确性;在所述数据采集阶段中的图像获取技术为采用图像采集设备,在车辆移动时逐行扫描地铁隧道表面,得到隧道表面图像;通过图像处理技术计算裂缝宽度,对裂缝进行识别、定位并且要对每条裂缝的宽度信息进行记录,建立裂缝数据库;在所述数据处理阶段中,根据裂缝数据库中的历史数据和上述步骤获得的最新裂纹数据,挖掘裂缝变化规律,考虑到环境复杂各种噪声的影响,利用裂缝宽度估计算法,通过裂缝宽度的预测步以及随后的更新步,再通过迭代步骤,来实时获得裂缝宽度的最佳估计值;其中,所述的裂缝宽度预测步包括裂缝宽度预测以及裂缝宽度协方差预测;所述的更新步同样包括对裂缝宽度及裂缝宽度协方差的更新;所述的数据校正阶段中,样本库中各个裂缝宽度的人工实际测量值分别为{C<sub>m1</sub>,C<sub>m2</sub>,C<sub>m3</sub>,C<sub>m4</sub>,C<sub>m5</sub>......},所述的估计值修正具体步骤为:(Ⅰ)对选取的每一个裂缝进行裂缝宽度的估计,得到对应的估计值{C<sub>e1</sub>,C<sub>e2</sub>,C<sub>e3</sub>,C<sub>e4</sub>,C<sub>e5</sub>......};(Ⅱ)计算每条裂缝的估计值与测量值的偏差,得到偏差集{C<sub>d1</sub>,C<sub>d2</sub>,C<sub>d3</sub>,C<sub>d4</sub>,C<sub>d5</sub>......},其中C<sub>di</sub>=C<sub>ei</sub>‑C<sub>mi</sub>;(Ⅲ)根据裂缝宽度的不同进行分类处理:宽度为W<sub>i</sub>的裂缝的估计的偏差值为<img file="FDA0000955399080000011.GIF" wi="669" he="85" />若<img file="FDA00009553990800000110.GIF" wi="235" he="78" />说明样本波动小,同一宽度的裂缝估计值偏差小;若<img file="FDA00009553990800000111.GIF" wi="234" he="70" />说明此裂缝的估计值波动较大,有些裂缝在机器识别时出现了错误,需要剔除;查找出现错误的估计值,依次剔除k<sub>ii</sub>,计算剩余样本的方差<img file="FDA0000955399080000012.GIF" wi="79" he="81" />若出现<img file="FDA0000955399080000013.GIF" wi="236" he="80" />则将剩余样本作为最终的样本,如果<img file="FDA0000955399080000014.GIF" wi="59" he="82" />始终大于<img file="FDA0000955399080000015.GIF" wi="103" he="73" />则选取最小<img file="FDA0000955399080000016.GIF" wi="57" he="82" />对应的样本值,再次进行剔除;重复之前的计算,直至方差小于<img file="FDA0000955399080000017.GIF" wi="103" he="74" />(Ⅳ)经过上述计算后,得到宽度为W<sub>i</sub>的裂缝估计的偏差值的最终样本集<img file="FDA0000955399080000018.GIF" wi="660" he="85" />计算该样本的平均值<img file="FDA0000955399080000019.GIF" wi="495" he="128" />由此得到判断值<img file="FDA0000955399080000021.GIF" wi="303" he="86" /><img file="FDA0000955399080000022.GIF" wi="75" he="69" />与<img file="FDA0000955399080000023.GIF" wi="66" he="86" />相对应;(Ⅴ)随着时间的变化,隧道中裂缝的宽度也会有所改变,定时对选取的裂缝再次进行测量,得到一组新的样本集,按照之前的处理方式,得到新的<img file="FDA0000955399080000024.GIF" wi="68" he="68" />和<img file="FDA0000955399080000025.GIF" wi="86" he="85" />需要再利用多次的测量结果得到一个准确的数据,作为最终的判定值和修订值<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><msubsup><mi>w</mi><mi>f</mi><mi>i</mi></msubsup><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><mover><msub><mi>w</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover><mo>*</mo><msub><mi>log</mi><mn>2</mn></msub><mfrac><mi>j</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac></mrow><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>log</mi><mn>2</mn></msub><mfrac><mi>j</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000955399080000026.GIF" wi="710" he="270" /></maths>其中<img file="FDA0000955399080000027.GIF" wi="74" he="94" />为最终的判定值,<img file="FDA0000955399080000028.GIF" wi="85" he="79" />表示宽度为W<sub>i</sub>的裂缝的第j次测量的判定值,n表示宽度为W<sub>i</sub>的裂缝的样本集中样本总数;由于时间的不同,裂缝会出现变化,所以对每次的判定值都添加了相应的权重,时间越久的数据所占的权重就越小,最近的数据所占的权重越大,第j次的数据权重用<img file="FDA0000955399080000029.GIF" wi="206" he="117" />来表示,随着j的增大,权重也相应的增大;同理可以得到最终的修订值:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><msubsup><mi>c</mi><mrow><mi>d</mi><mi>f</mi></mrow><mi>i</mi></msubsup><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><mover><msubsup><mi>c</mi><msub><mi>d</mi><mi>J</mi></msub><mi>i</mi></msubsup><mo>&OverBar;</mo></mover><mo>*</mo><msub><mi>log</mi><mn>2</mn></msub><mfrac><mi>j</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac></mrow><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>log</mi><mn>2</mn></msub><mfrac><mi>j</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA00009553990800000210.GIF" wi="687" he="269" /></maths>(Ⅵ)根据所得到的<img file="FDA00009553990800000211.GIF" wi="73" he="94" />和<img file="FDA00009553990800000212.GIF" wi="88" he="94" />来对估计值进行修订,设估计值为C,找到与C最相近的判定值<img file="FDA00009553990800000213.GIF" wi="93" he="94" />得到裂缝最终的宽度为:<img file="FDA00009553990800000214.GIF" wi="411" he="96" />
地址 100044 北京市海淀区西直门外上园村3号