本发明属于图像处理方法技术领域，涉及一种基于图像局部对比度的能见度测量方法。
背景技术
近年来，雾霾天气的频发给人们的日常生活带来诸多不便，大气能见度及时准确的预报对航空，河运和高速公路等交通运输有重要意义。现有的能见度测量方法主要有目测法，器测法和数字摄像法。目测法规范性、客观性相对较差；器测法成本高、调试复杂、稳定性易受环境限制；数字摄像法利用图像的视觉特征测量大气能见度，克服了目测法和器测法的局限性，成本低应用广，受到国内外学者的广泛关注。
陶善昌等通过计算图像中地平线附近两个不同距离的目标物与背景天空亮度差的比值求取白天能见度，但该方法对测量环境与目标物的选择要求较高。许茜等通过提取基于像素对比度的图像特征及向量，用支持向量机对能见度模型进行训练，算法复杂，计算量大。baumer等提出根据图像局部对比度强度来估计视觉范围，但需对大量的目标物进行标定，且系统对采集设备的参数设置很敏感，鲁棒性较差。
技术实现要素：
本发明的目的是提供一种基于图像局部对比度的能见度测量方法，能够更加精确的测量雾霾天气的能见度。
本发明采用的技术方案是，基于图像局部对比度的能见度测量方法，其特征在于，具体按照下述步骤进行：
步骤1，选定目标物，用图像采集设备采集选定的目标物的无雾霾图像，备用；
在采集目标物的无雾霾图像的相同位置，使用相同的图像采集设备采集目标物的有雾霾图像；
步骤2，将步骤1得到的无雾霾图像进行预处理得到有雾霾基准图像，将步骤1得到的有雾霾图像进行预处理得到无雾霾基准图像；
步骤3，计算步骤2得到的无雾霾基准图像的图像梯度cj(x,y)和有霾基准图像的图像梯度ji(x,y)；
步骤4，步骤3完成后，建立大气能见度对比特征模型计算雾霾天气的能见度：
其中，(x,y)为无雾霾基准图像的图像坐标，(x,y)为有霾基准图像的图像坐标，r为景深，v为大气能见度。
本发明的特点还在于：
步骤2中通过下述步骤对无雾霾图像进行预处理：
将步骤1得到的无雾霾图像进行灰度处理后得到无雾霾灰度图像，将无雾霾灰度图像进行感兴趣区域提取后得到无雾霾感兴趣图像，将无雾霾感兴趣图像进行特征点提取并进行特征点匹配后得到无雾霾基准图像；
所述步骤2中通过下述步骤对有雾霾图像进行预处理：
将步骤1得到的有雾霾图像进行灰度处理后得到无有雾霾灰度图像，将无有雾霾灰度图像进行感兴趣区域提取后得到有雾霾感兴趣图像，将有雾霾感兴趣图像进行特征点提取并进行特征点匹配后得到有雾霾基准图像。
无雾霾感兴趣图像和有雾霾感兴趣图像进行特征点提取时均使用sift算法进行角点、边缘点、暗区的亮点及亮区的暗点进行提取，提取后均通过计算特征点的方向进行图像匹配。
步骤3中计算有雾霾基准图像的图像梯度ji(x,y)按照下述步骤进行：
步骤3.1，将步骤3得到的有雾霾基准图像进行分区，得到多个有雾霾基准分区图像；
步骤3.2，在每个有雾霾基准分区图像中均任选一像素点δi(xi,yi),将每个像素点δi(xi,yi)在0°、45°、90°和135°方向进行卷积，得到点δi(xi,yi)在0°方向的卷积结果si1,点δi(xi,yi)在45°方向的卷积结果si2,点δi(xi,yi)在90°方向的卷积结果si3和将点δi(xi,yi)在135°方向的卷积结果si4,将每个象素点的在0°、45°、90°和135°方向的卷积结果求取平均值,即为该点梯度值gi(xi,yi)：
si1＝g(xi,yi-2)+g(xi,yi-1)+g(xi,yi)+g(xi,yi+1)+g(xi,yi+2)(2)；
si2＝g(xi-2,yi)+g(xi-1,yi)+g(xi,yi)+g(xi+1,yi)+g(xi+2,yi)(3)；
si3＝g(xi-2,yi-2)+g(xi-1,yi-1)+g(xi,yi)+g(xi+1,yi+1)+g(xi+2,yi+2)
(4)；
si4＝g(xi+2,yi-2)+g(xi+1,yi-1)+g(xi,yi)+g(xi-1,yi+1)+g(xi-2,yi+2)
(5)；
gi(xi,yi)＝mean(si1，si2，si3，si4)(6)；
其中，i＝1,2,....,n，n为有雾霾基准分区图像的个数；
步骤3.3，将每个有雾霾基准分区图像中选取的所有像素点取平均值即可得到无雾霾基准图像的图像梯度cl(x，y)：
ji(x,y)＝mean(g1(x1,y1),g2(x2,y2),...,gn(xn,yn)(7)。
本发明的有益效果在于：
本发明的基于图像局部对比度的能见度测量方法，有效地提高了能见度检测精度，能够更加准确地测量雾霾天气的能见度。
附图说明
图1为本发明基于图像局部对比度的能见度测量方法对雾霾图像进行roi提取的示意图；
图2为本发明图像局部对比度的能见度测量方法中未匹配到特征点的示意图；
图3为本发明图像局部对比度的能见度测量方法对雾霾头像进行分区的示意图；
图4为本发明像局部对比度的能见度测量方法实施例中计算结果与cjy-1g能见度测量结果对比图。
具体实施方式
下面根据附图和具体实施方式对本发明进行详细说明。
基于图像局部对比度的能见度测量方法，具体按照下述步骤进行：
步骤1，选定目标物，在空气质量指数(aoi)为优时，基本无空气污染，即为无雾霾天气时，使用图像采集设备采集选定的目标物的无雾霾图像，备用；
在采集目标物的无雾霾图像的相同位置，使用相同的图像采集设备采集目标物的有雾霾图像；
步骤2，将步骤1得到的无雾霾图像进行灰度处理后得到无雾霾灰度图像，如图1所示，将无雾霾灰度图像进行感兴趣区域提取后得到无雾霾感兴趣图像，将无雾霾感兴趣图像，使用sift算法进行角点、边缘点、暗区的亮点及亮区的暗点进行提取，提取后均通过计算特征点的方向进行图像匹配后得到无雾霾基准图像，如图2所示，若是匹配搭配的特征点较少时，将无雾霾感兴趣图像展宽后再进行角点、边缘点、暗区的亮点及亮区的暗点，提取后再进行匹配得到无雾霾基准图像；
将步骤1得到的有雾霾图像进行灰度处理后得到有雾霾灰度图像，将有雾霾灰度图像进行感兴趣区域提取后得到有雾霾感兴趣图像，将有雾霾感兴趣图像使用sift算法进行角点、边缘点、暗区的亮点及亮区的暗点进行提取，提取后均通过计算特征点的方向进行图像匹配后得到有雾霾基准图像；如图2所示，若是匹配搭配的特征点较少时，将有雾霾感兴趣图像展宽后再进行角点、边缘点、暗区的亮点及亮区的暗点，提取后再进行匹配得到有雾霾基准图像
步骤3，按照下述方法计算步骤2得到的无雾霾基准图像的图像梯度cj(x,y)：
步骤a，如图3所示，将步骤2得到的无雾霾基准图像使用3×16的栅格进行分区，得到48个无雾霾基准分区图像；
步骤b，在每个无雾霾基准分区图像中均任选一像素点δj(xj,yj),将每个像素点δj(xj,yj)在0°、45°、90°和135°方向进行卷积，得到点δj(xj,yj)在0°方向的卷积结果sj1,点δj(xj,yj)在45°方向的卷积结果sj2,点δj(xj,yj)在90°方向的卷积结果sj3和将点δj(xj,yj)在135°方向的卷积结果sj4,将每个象素点的在0°、45°、90°和135°方向的卷积结果求取平均值,即为该点梯度值δj(xj,yj)：
sj1＝g(xj,yj-2)+g(xj,yj-1)+g(xj,yj)+g(xj,yj+1)+g(xj,yj+2)(8)；
sj2＝g(xj-2,yj)+g(xj-1,yj)+g(xj,yj)+g(xj+1,yj)+g(xj+2,yj)(9)；
sj3＝g(xj-2,yj-2)+g(xj-1,yj-1)+g(xj,yj)+g(xj+1,yj+1)+g(xj+2,yj+2)
(10)；
sj4＝g(xj+2,yj-2)+g(xj+1,yj-1)+g(xj,yj)+g(xj-1,yj+1)+g(xj-2,yj+2)
(11)；
gj(xj,yj)＝mean(sj1，sj2，sj3，si4)(12)；
其中，j＝1,2,....,48；
步骤c，将每个有雾霾基准分区图像中选取的所有像素点取平均值即可得到无雾霾基准图像的图像梯度cj(x,y)：
cj(x,y)＝mean(g1(x1,y1),g2(x2,y2),...,g48(x48,y48)(13)；
按照下述方法计算步骤2得到的有雾霾基准图像的图像梯度ji(x,y)：
步骤3.1，如图3所示将步骤2得到的有雾霾基准图像使用3×16的栅格进行分区，得到48个有雾霾基准分区图像；
步骤3.2，在每个有雾霾基准分区图像中均任选一像素点δi(xi,yi),将每个像素点δi(xi,yi)在0°、45°、90°和135°方向进行卷积，得到点δi(xi,yi)在0°方向的卷积结果si1,点δi(xi,yi)在45°方向的卷积结果si2,点δi(xi,yi)在90°方向的卷积结果si3和将点δi(xi,yi)在135°方向的卷积结果si4,将每个象素点的在0°、45°、90°和135°方向的卷积结果求取平均值,即为该点梯度值gi(xi,yi)：
si1＝g(xi,yi-2)+g(xi,yi-1)+g(xi,yi)+g(xi,yi+1)+g(xi,yi+2)(2)；
si2＝g(xi-2,yi)+g(xi-1,yi)+g(xi,yi)+g(xi+1,yi)+g(xi+2,yi)(3)；
si3＝g(xi-2,yi-2)+g(xi-1,yi-1)+g(xi,yi)+g(xi+1,yi+1)+g(xi+2,yi+2)
(4)；
si4＝g(xi+2,yi-2)+g(xi+1,yi-1)+g(xi,yi)+g(xi-1,yi+1)+g(xi-2,yi+2)
(5)；
gi(xi,yi)＝mean(si1，si2，si3，si4)(6)；
其中，i＝1,2,....,48；
步骤3.3，将每个有雾霾基准分区图像中选取的所有像素点取平均值即可得到无雾霾基准图像的图像梯度cl(x，y)：
ji(x,y)＝mean(g1(x1,y1),g2(x2,y2),...,g48(x48,y48)(7)；
步骤4，建立大气能见度对比特征模型计算雾霾天气的能见度：
其中，r为景深，v为大气能见度。
实施例
使用m5’回归树对气能见度对比特征模型进行训练，使用10-foldercross-validation对实验数据进行交叉验证。
随机选取1500组数据样本，能见度最低为625m，最高4761m。以安装在西安理工大学教6楼楼顶的前向散射仪cjy-1g测量的能见度数据为真值，用本发明的方法与m5’模型树结合的测量值与真值进行拟合，如图4所示，实验样本中有1373个点落在±20％的误差线内，占实验样本的91.53％，其中，实验值与真值的相对误差最大为35.17％，最小为0。实验结果表明，本文方法的测量值与前向散射仪测量的能见度值基本保持一致，能够较为准确地对能见度进行测量，简单方便。