一种地铁隧道管片渗漏水的远距离测量方法
【技术领域】
[0001] 本发明涉及地铁隧道管片检测领域，具体涉及一种地铁隧道管片渗漏水的远距离 测量方法。
【背景技术】
[0002] 隧道渗漏水现象在运营隧道中非常普遍，是一种常见的病害，若不及时对其进行 修补，将会降低管片结构的强度，并引发钢筋和螺栓腐蚀、管片开裂及混凝土剥落等其他病 害，会严重危害隧道的运营并且可能导致安全事故。
[0003] 目前，国内传统的盾构隧道渗漏水检测方法，是安排受过培训的专业人员对隧道 进行定期检测，根据肉眼观察结果进行判断，其工作量大、耗费人力和物力多且检测精度不 高。随着科技的进步，无损检测的方法应运而生。计算机高性能处理器、大容量存储器以及 图像处理技术的飞速发展，使得基于图像处理技术的隧道衬砌表面病害自动检测技术成为 了可能。但其在隧道病害检测方面的研究还是比较少而且不成熟。黄永杰等提出了基于数 字图像自动识别算法的盾构隧道渗漏水自动检测技术，采用自主研发的检测系统进行现场 渗漏信息的采集、识别、抓取和处理而完成盾构隧道的渗漏水检测。但是处理对象是近距离 拍摄，对于远距离的图像效果并不具有适用性。刘学增等以隧道衬砌渗漏水面积为检测目 标，研究出包括去噪、锐化、分割、修正的一整套数字图像处理算法，但研究成果对照片质量 依赖性较大，针对复杂的图像处理效果并不理想。综上，现有检测技术存在照片依赖性大， 适用范围小等问题，研究一种效率高、检测精确的地铁隧道管片渗漏水远距离检测算法是 十分有必要的。
【发明内容】
[0004] 针对现有技术存在的缺陷，本发明的目的在于提供一种地铁隧道管片渗漏水的远 距离测量方法，从而实现在昏暗的地铁隧道环境下对管片渗漏水的面积的实时、高效和精 确的检测。
[0005] 为了达到以上目的，本发明采用以下技术方案：
[0006] -种地铁隧道管片渗漏水的远距离测量方法，具体包括如下步骤：
[0007] a.将(XD相机置于距离地铁隧道管片6-10米处，正对地铁隧道管片拍摄图像；
[0008] b.根据步骤a采集的彩色图像，将彩色图像进行灰度化处理，灰度化公式如下：
[0009] f = 0.3R+0.59G+0.11B
[0010] 其中，f为图像的灰度值，R，G，B为真彩图像每个像素点的3个参数值；
[0011] c.根据步骤b灰度化后的图像，通过灰度拉伸来提升图像中渗漏水与背景的对比 度，灰度拉伸公式如下： - Ο Ο <./' < a
[0012] /-< a < f < h b-a 255 bSf <25.5
[0013] 上述公式中a，b分别为上下阈值，f'为灰度拉伸后的灰度值；
[0014] d.根据步骤c得到预处理后的图像，采用改进迭代法对地铁隧道管片渗漏水图像 进行二值分割；
[0015] e.根据步骤d得到的二值图像，采用基于连通区域的多级滤波算法处理，具体包括 三个部分，分别为基于连通区域小面积滤波，基于连通区域的圆形度滤波以及特殊噪声滤 波；
[0016] f.对于步骤c得到预处理后的灰度图像，对地铁隧道管片渗漏水图像进行3X2分 块划分，针对每一块区域采用〇stu最大类间方差法图像分割；
[0017] g.经过步骤f分块阈值分割后的图像，能够保留管片的边缘信息，为了提取管片边 缘位置，用Canny算子对其进行边缘检测，将地铁隧道管片渗漏水图像的边缘全部提取出 来，为了有针对的得到管片的边缘在竖直方向，有如下附加条件公式：
[0018] Pl.x-P2.x| <TC
[0019] 其中PI. x和P2.x分别是直线两端的横坐标，将两者做差的绝对值小于阈值TC的直 线提取出来；
[0020] h.根据步骤g提取出来的管片边缘直线，得到边缘直线间的像素点数，同时已知管 片的宽度，对应得到每个像素点的实际宽度，也算出图像上每个像素点的实际面积，进一步 验证测量精度；
[0021] i.根据步骤e通过多级滤波后得到的渗漏水区域，由步骤h得到的换算方法直接换 算成实际面积。
[0022] 所述步骤d中的采用改进迭代法对地铁隧道管片渗漏水图像进行二值分割，具体 步骤如下：
[0023] 1)设定灰度值偏大区域自动转化成背景区域，其阈值为，最小灰度值为0,令 初始阈值为：
[0024] To= {Tk I k = 0}
[0025] To=(0 + T/max)/2
[0026] 2)利用阈值Tk将图像分成两组，Ri和R2，其中
[0027] Ri = {f(x,y) |〇 < f(x,y) < Tk}
[0028] R2 = {f (x,y) I Tk+1 < f (x,y) < T7 max}
[0029] 3)计算区域RjPR2的平均灰度值ZjPZ2,其中 2],/(/"./) x,V(/,./)
[0030] Zj -
[0031] %= - X-Vu,./)
[0032] 式中，f(i，j)是图像上点（i，j)的灰度值，N(i，j)是点（i，j)的权重系数，一般设为 1；
[0033] 4)选择新的阈值：
[0034] Tk+i = (Zi+Z2)/2
[0035] 5)如果|Tk+1_Tk| < 1，则迭代结束，否则k = k+l，跳转步骤2)继续迭代。
[0036] 所述步骤e的采用基于连通区域的多级滤波算法，具体为：
[0037] 1)基于二值图像连通区域小面积滤波
[0038] 渗漏水在二值图像中一般是不规则的团状分布，相比于分布的点状噪声面积更 大，设Pk(x，y)为二值图像的连通区域，N为连通区域的个数，S k为每个连通区域的面积，设定 一个阈值Tz，滤波公式如下： 0, Sk < T,
[0039] Α(Χ,>〇 = < 1 苴他；1=1,2:,3, " ·,况
[0040] 2)基于连通区域圆形度滤波
[0041] 二值图像中存在着大量的块状噪声，其中含有螺栓孔属于圆形结构，与渗漏水的 不规则结构并不一样，因此利用连通域的圆形度进行滤波，圆形度计算公式如下：
[0042] (6=-- cl
[0043] 其中Ck为连通区域的面积，理想情况下，对于圆形，e的值等于l，e越小越不规则， 设定一个阈值T b，滤波公式如下： ? 0, e > Th
[0044] Pk {x；, y) - < . H /, ； k = 1,2,3:, * · ·, ^ ? 1 1，其他
[0045] 3)特殊噪声滤波
[0046] 前两种滤波处理后还存在特殊噪声，由于它的宽度w或者高度h变化不大，逐行扫 描图像计算连通域的宽度和高度，统计其宽度或者高度众数出现的次数，分别为N w和Nh，设 定两个阈值Iw和Th，滤波公式如下所示： ? 0, Nv, > Tw or Nh >Th
[0047] = 1 , 甘仙 ；Κ，2,3，·.?
[1, 其他 ο
[0048] 所述步骤h的验证测量精度是通过根据第一块管片的单个像素实际长度来验证第 二块管片来实现，具体如下：
[0049] 其中管片的实际宽度均为W米，第一块管片经计算像素点数为N，第二块管片像素 点数为Λ，按照提换算公式得到管片2所测得的宽度为： w
[0050] W7 =--N, N
[0051] 第二块管片的实际的宽度W米，计算实际管片的面积误差为： IW72 -W2 I
[0052] σ ^ ~-~~r--U lOO% IV ~ 〇.
[0053] 本发明与现有技术相比较，具有如下显而易见的突出实质性特点和显著优点： [0054]本发明对于地铁隧道管片渗漏水面积测量无需人工参与，提高工作效率，适合在 实时的系统中采用。采用面测量，同时对于渗漏水面积大小可以通过程序换算实时标定，减 小劳动强度。因此，该检测算法效率高、检测精确。
【附图说明】
[0055]图1是本发明的算法的总流程图。
[0056]图2是管片渗漏水图像的灰度图。
[0057]图3是通过改进迭代法处理后的二值图像。
[0058]图4是地铁渗漏水图像中的圆形噪声以及特殊噪声。
[0059]图5是管片边缘特征提取图。
[0060]图6是通过算法最终提取的渗漏水区域。
【具体实施方式】
[0061]下面结合附图和具体实施例对本发明的技术方案作进一步具体说明，
[0062] 如图1所示，一种地铁隧道管片渗漏水的远距离测量方法，具体包括如下步骤： [0063] a.将(XD相机置于距离地铁隧道管片8米处，正对地铁隧道管片拍摄图像；
[0064] b.根据步骤a采集的彩色图像，将彩色图像进行灰度化处理，灰度化公式如下：
[0065] f = 0.3R+0.59G + 0.11B
[0066] 其中，f为图像的灰度值，R，G，B为真彩图像每个像素点的3个参数值，如图2所示；
[0067] c.根据步骤b灰度化后的图像，通过灰度拉伸来提升图像中渗漏水与背景的对比 度，灰度拉伸公式如下： " 0 0</<?： ， 255
[0068] f--(f-a) a < f < h b - a 255 h < f < 255
[0069] 上述公式中a，b分别为上下阈值，f'为灰度拉伸后的灰度值;本实施例中a = 20，b = 200；
[0070] d.根据步骤c得到预处理后的图像，采用改进迭代法对地铁隧道管片渗漏水图像 进行二值分割，如图3所示；
[0071] e.根据步骤d得到的二值图像，采用基于连通区域的多级滤波算法处理，具体包括 三个部分，分别为基于连通区域小面积滤波，基于连通区域的圆形度滤波以及特殊噪声滤 波，如图4所示；
[0072] f.对于步骤c得到预处理后的灰度图像，对地铁隧道管片渗漏水图像进行3X2分 块划分，针对每一块区域采用〇stu最大类间方差法图像分割；
[007