一种基于夜光遥感数据的道路夜间照明质量监测方法

1.本发明涉及道路照明质量监测技术领域，具体而言涉及一种基于夜光遥感数据的道路夜间照明质量监测方法。


背景技术：

2.道路夜间照明是现代交通设施的重要组成部分，最重要的目的是为夜间出行车辆以及行人创造良好的光照环境。道路夜间照明的质量对车辆行驶安全有很大影响，过暗的照明条件降低了夜间能见度，存在交通隐患。而亮度变化过于大或者对比过于强烈，也会干扰驾驶员的反应速度和判断力。提高道路照明质量，对于保障夜间交通安全、提高交通效率具有重要意义。根据国际照明委员会(cie) 报告，良好的夜间照明能让城市道路交通事故率降低30%。
3.为了营造良好的道路夜间照明环境，需要对道路照明质量进行有效监测。传统的道路照明质量监测方法是布点法，在选定路段按照一定间隔并利用照度计或者亮度计测量若干样点的照度或者亮度，然后计算平均照度或者平均亮度作为道路照明水平的表征，计算照度或者亮度最大值与最小值之差作为道路照明均匀性的表征。这样的测量方式存在一些缺陷：1）测量地点有限，只能以点代面，不能快速监测大范围下各条道路的照明质量；2）需要耗费大量的人力物力；3）测量作业会影响正常交通。夜间灯光遥卫星遥感的发展为道路照明质量监测提供了新的思路。与布点监测相比，卫星遥感能够在大范围内快速监测道路照明质量，并且不会影响交通，具有独特的优势。目前我国长光公司的jl1-3b遥感卫星能够提供空间分辨率高达0.92m的夜间灯光遥感图像，具备了道路照明监测的潜力。然而目前尚未有利用夜间灯光遥感的道路照明质量监测相关报道。


技术实现要素：

4.本发明针对现有技术中的不足，提供一种基于夜光遥感数据的道路夜间照明质量监测方法，能够基于夜间灯光遥感数据快速有效进行道路夜间照明质量的监测与评估，为改善道路夜间照明环境、保障交通安全提供技术支撑。
5.为实现上述目的，本发明采用以下技术方案：本发明实施例提出了一种基于夜光遥感数据的道路夜间照明质量监测方法，所述监测方法包括以下步骤：s1，对jl1-3b夜间灯光遥感数据进行去噪和辐射定标处理，将灰度值转换为辐射亮度值，得到相应的辐射亮度图；以地形图数据或者具有准确定位的昼间高分辨率遥感图像为基准图像，对辐射定标后的辐射亮度图进行几何精校正，消除空间定位误差，得到预处理后的jl1-3b遥感图像；s2，在卫星过境同时，选取典型区域作为观测点开展地表照度实地观测，记录观测点的经纬度数据和实测照度；提取观测点对应的jl1-3b遥感图像像元的第1、2、3波段辐射亮度值进行求和，再与实测照度进行对数回归，构建地面照度遥感反演模型：
式中，为照度值，单位为lux；、、分别为jl1-3b遥感图像像元的第1、2、3波段的辐射亮度值；和为经验系数，通过回归拟合得到；将构建的地面照度遥感反演模型应用于预处理后的jl1-3b遥感图像，得到研究区的夜间地表照度空间分布；s3，基于高分辨率昼间遥感图像，得到各条道路的无遮挡空间范围；s4，分别提取各条道路无遮挡空间范围内所有像元的地表照度值，计算各条道路的照度平均值、标准差与变异系数，对道路照明质量进行评估。
6.进一步地，步骤s1中，得到预处理后的jl1-3b遥感图像的过程包括以下子步骤：s11，对jl1-3b夜间灯光遥感数据的第1、2、3波段进行rgb真彩色合成，再进行hsv变换，按照下面的判断准则提取并标识噪声像元：式中，h、s、v分别为hsv颜色空间的色调、饱和度和亮度值；s12，针对每个噪声像元，分别以周边非噪声像元第1、2、3波段值对相应波段进行空间插补，以插补值替换噪声像元值，去除噪声；s13，对去噪后的jl1-3b夜间灯光遥感数据进行辐射定标处理，将灰度值转换为辐射亮度值，得到jl1-3b辐射亮度图：式中，为第i波段辐射亮度值，单位为w/(m2
·
sr)，为第i波段灰度值，和为第i波段为定标系数，i=1,2,3；定标系数和的值分别为：，，，，，。
7.s14，以地形图数据或者具有准确定位的昼间高分辨率遥感图像为基准图像对辐射定标后jl1-3b辐射亮度图进行几何精校正，消除空间定位误差。
8.进一步地，步骤s2中，在卫星过境同时，选取典型区域作为观测点开展地表照度实地观测，记录观测点的经纬度数据和实测照度的过程包括以下子步骤：s21，在地面选取若干个具有空旷无遮挡特性的典型区域作为观测点开展照度实地测量，观测时间为卫星过境时间前后1小时；s22，将照度计的感光头水平朝上放置，记录照度值，并利用gnss设备记录观测点经纬度值；s23，结合人工筛选，使用3σ法剔除照度观测数据中的异常值。
9.进一步地，步骤s2中，提取各观测点的照度值以及对应的预处理后jl1-3b夜光遥感图像像元的第1、2、3波段辐射亮度值，通过最小二乘法回归拟合求解系数和的值，构建得到夜间地表照度遥感反演模型，并通过10折交叉验证对模型精度进行验证。
10.进一步地，步骤s3中，基于高分辨率昼间遥感图像，得到各条道路的无遮挡空间范围的过程包括以下子步骤：s31，基于高分辨率昼间遥感图像计算归一化植被指数ndvi：式中，ndvi为归一化植被指数值；nir和r分别为遥感图像近红外波段和红波段像元值；s32，将高分辨率昼间遥感图像的红、绿、蓝波段进行真彩色合成，采用gram-schmidt正交变换融合方法与ndvi图像进行融合，在真彩色图像基础上突出植被信息；s33，基于融合后图像，通过目视解译提取道路空间范围，并提取道路范围内的行道树冠层范围，去除道路中的行道树冠层遮挡范围，得到各条道路的无遮挡空间范围。
11.进一步地，步骤s4中，分别提取各条道路无遮挡空间范围内所有像元的地表照度值，计算各条道路的照度平均值、标准差与变异系数的过程包括以下子步骤：s41，采用下述公式计算得到指定道路的照度平均值：式中，为指定道路的照度平均值，单位为lux；为该指定道路范围内第j个像元的照度值，单位为lux；n为该指定道路范围内像元总数；s42，采用下述公式计算得到指定道路的照度标准差：式中，为指定道路的照度标准差，单位为lux；为该指定道路的照度平均值，单位为lux；s43，采用下述公式计算得到指定道路的照度变异系数：式中，为指定道路的照度变异系数。
12.s44，按照不同道路的不同等级，根据它们的照度平均值、照度标准差和照度变异系数对道路照明质量进行统计分析和评价；其中，照度平均值用于反映道路的平均照明亮度水平，标准差与变异系数用于反映道路照明亮度的均匀度。
13.本发明基于0.92m分辨率的jl1-3b夜间灯光遥感数据开展道路夜间照明质量监测，在反演夜间地面照度空间分布图的基础上，结合道路空间分布信息计算各条道路照度平均值、标准差与变异系数，对道路照明质量进行监测和定量评估。
14.本发明的有益效果是：
本发明提出的基于夜光遥感数据的道路夜间照明质量监测方法，利用jl1-3b夜间灯光遥感数据反演大范围的夜间地面照度，结合道路空间范围信息对道路夜间照明质量进行监测和评估。相比之下，传统的地面布点监测方法只能够对少数道路的有限地点进行监测，缺乏效率并且会妨碍交通。本发明提供的方法能对大范围的道路夜间照明质量进行有效监测与评估，为保障夜间交通安全、提高交通效率提供技术支撑，测量效率明显优于传统监测方法，并且不会对交通造成影响。
附图说明
15.图1为本发明实施例的基于夜光遥感数据的道路夜间照明质量监测方法流程图。
16.图2a为本发明实施例的去噪前的jl1-3b图像示意图。
17.图2b为本发明实施例的去噪后的jl1-3b图像示意图。
18.图3为本发明实施例的经过去噪、辐射定标和几何精校正的jl1-3b夜光遥感图像示意图。
19.图4为本发明实施例的地面照度遥感反演交叉验证散点图。
20.图5为本发明实施例的遥感反演夜间地表照度空间分布图。
21.图6为本发明实施例的道路提取结果示意图。
具体实施方式
22.现在结合附图对本发明作进一步详细的说明。
23.需要注意的是，发明中所引用的如“上”、“下”、“左”、“右”、“前”、“后”等的用语，亦仅为便于叙述的明了，而非用以限定本发明可实施的范围，其相对关系的改变或调整，在无实质变更技术内容下，当亦视为本发明可实施的范畴。
24.图1为本发明实施例的基于夜光遥感数据的道路夜间照明质量监测方法流程图。参见图1，该监测方法包括以下步骤：s1，对jl1-3b夜间灯光遥感数据进行去噪和辐射定标处理，将灰度值转换为辐射亮度值，得到相应的辐射亮度图；以地形图数据或者具有准确定位的昼间高分辨率遥感图像为基准图像，对辐射定标后的辐射亮度图进行几何精校正，消除空间定位误差，得到预处理后的jl1-3b遥感图像。
25.s2，在卫星过境同时，选取典型区域作为观测点开展地表照度实地观测，记录观测点的经纬度数据和实测照度；提取观测点对应的jl1-3b遥感图像像元的第1、2、3波段辐射亮度值进行求和，再与实测照度进行对数回归，构建地面照度遥感反演模型：式中，为照度值，单位为lux；、、分别为jl1-3b遥感图像像元的第1、2、3波段的辐射亮度值；和为经验系数，通过回归拟合得到；将构建的地面照度遥感反演模型应用于预处理后的jl1-3b遥感图像，得到研究区的夜间地表照度空间分布。
26.s3，基于高分辨率昼间遥感图像，得到各条道路的无遮挡空间范围。
27.s4，分别提取各条道路无遮挡空间范围内所有像元的地表照度值，计算各条道路
的照度平均值、标准差与变异系数，对道路照明质量进行评估。
28.本实施例的监测方法在对jl1-3b夜间灯光遥感进行去噪、辐射定标、几何精校正处理的基础上，结合实测数据构建夜间地表照度遥感模型反演地面照度，结合道路空间分布图计算各条道路夜间照度的平均值、标准差与变异系数，评估道路的照明质量。具体包括以下步骤：1）jl1-3b夜间灯光遥感数据处理jl1-3b夜间灯光遥感数据存在较多的点状纯色噪声，需要进行去噪处理。对jl1-3b数据的第1、2、3波段进行rgb真彩色合成，然后进行hsv变换，按照下面的判断准则判断是否是噪声像元：式中，h、s、v分别为hsv颜色空间的色调、饱和度和亮度值。
29.按照上述标准提取并标识噪声像元，分别以周边非噪声像元第1、2、3波段值对相应波段进行空间插补，以插补值替换噪声像元值，去除噪声。
30.对去噪后的jl1-3b夜间灯光遥感数据进行辐射定标处理，将灰度值转换为辐射亮度值：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（1）式中，为第i波段辐射亮度值，单位为w/(m2
·
sr)，为第i波段灰度值，和为第i波段为定标系数（数值见表1）。
31.表1 jl1-3b夜间灯光遥感数据定标系数值以地形图数据或者具有准确定位的昼间高分辨率遥感图像为基准图像对辐射定标后jl1-3b辐射亮度图进行几何精校正，消除空间定位误差。
32.2）地表照度遥感反演在卫星过境同时，选取典型区域开展地表照度实地观测。将照度计的感光头水平放置，记录照度值，并记录观测点经纬度数据。提取实测点对应jl1-3b遥感图像像元的第1、2、3波段（即r、g、b波段）辐射亮度值求和，与实测照度进行对数回归，构建地面照度遥感反演模型：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（2）式中，e为照度值，单位为lux；、、分别为jl1-3b第1、2、3波段的辐射亮度值；和为经验系数，需要通过回归拟合得到。
33.将构建的地面照度遥感反演模型应用于预处理后的jl1-3b遥感图像，得到研究区的夜间地表照度空间分布。
34.3）道路照明质量评估基于高分辨率昼间遥感图像计算归一化植被指数ndvi
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（3）式中，ndvi为归一化植被指数值；nir和r分别为遥感图像近红外波段和红波段像元值。
35.将高分辨率昼间遥感图像的红、绿、蓝波段进行真彩色合成，然后采用gram-schmidt正交变换融合方法与ndvi图像进行融合，在真彩色图像基础上突出植被信息。
36.基于融合后图像，通过目视解译提取道路空间范围，并提取道路范围内的行道树冠层范围，去除道路中的行道树冠层遮挡范围，得到各条道路的无遮挡空间范围。
37.4）道路照明质量评估分别提取各条道路无遮挡空间范围内所有像元的地表照度值，计算各条道路的照度平均值、标准差与变异系数：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（4）式中，为某条道路的照度值平均，单位为lux；为该道路范围内第j个像元的照度值，单位为lux；n为该道路范围内像元总数。
38.ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（5）式中，为某条道路的照度标准差，单位为lux；为该道路的照度值平均，单位为lux；为该道路范围内第j个像元的照度值，单位为lux；n为该道路范围内像元总数。
39.ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（6）式中，为某条道路的照度变异系数；为该道路的照度标准差，单位为lux；为该道路的照度值平均，单位为lux。
40.平均值反映了道路的平均照明亮度水平，标准差与变异系数反映了道路照明亮度的均匀度。按照不同道路的不同等级，根据它们的照度平均值、标准差和变异系数对道路照明质量进行统计分析和评价。
41.实例本实例以南京主城区为研究区，利用jl1-3b夜光遥感数据监测南京主城区各道路夜间照明质量，以下为示例的具体实施步骤。
42.步骤1，对南京主城区范围内的jl1-3b夜光遥感图像进行去噪处理。jl1-3b夜光遥感图像有3个波段：第1波段（红波段）、第2波段（绿波段）和第3波段（蓝波段）。各个波段都存在一定比例的亮点噪声，进行rgb真彩色合成之后点状噪声体现为纯色（红色、绿色或者蓝
色）像元，详见图2a。将rgb图像进行hsv变换，转换到hsv颜色空间。然后利用下面的判断准则标识噪声像元：式中，h、s、v分别为hsv颜色空间的色调、饱和度和亮度值。
43.按规则识别出噪声像元后，将噪声像元掩膜为无效值，然后运用三角网空间插值通过非噪声像元进行空间插值，以噪声像元空间插值结果作为这些噪声像元的插补值，完成图像的去噪处理，详见图2b。
44.步骤2，基于jl1-3b的辐射定标参数对经过去噪处理的夜光遥感图像进行辐射定标，将dn值转化为辐射亮度值。
45.式中，li分别为波段i（i=1,2,3）的辐射亮度值，dni分别为波段i（i=1,2,3）的灰度值。
46.步骤3，jl1-3b夜间灯光遥感数据存在几何畸变。以具有准确空间定位信息的昼间高分辨率遥感图像（如gf2、worldview等）为基准图像，以jl1-3b夜光图像为待校正图像，满幅选取一定数量的控制点进行几何精校正，控制点的数量大于30。由于jl1-3b遥感图像为面阵成像，单景图像为中心投影，研究区图像由多景图像镶嵌而成，整个图像的几何畸变存在局部差异，无法使用统一的坐标转换方程拟合，需要采用局部三角网模型进行坐标转换，得到具有准确空间定位的jl1-3b夜光遥感图像，详见图3。
47.步骤4，在jl1-3b卫星过境当晚，在地面选取若干典型地点开展照度实地测量，观测时间为卫星过境时间前后1小时，观测仪器为tes-1399r照度计，其测量范围为0.01 lux至999900 lux。选取空旷的、没有树木等遮挡的路面作为观测点，将照度计的感光头水平朝上放置，记录照度值，并利用gnss设备记录观测点经纬度值。结合人工筛选，使用3σ法剔除照度观测数据中的异常值。
48.步骤5，提取各观测点的照度值以及对应的预处理后jl1-3b夜光遥感图像像元的第1、2、3波段辐射亮度值，基于公式（2）通过最小二乘回归拟合求解系数β0和β1的值，构建夜间地表照度遥感反演模型，并通过10折交叉验证对模型精度进行验证。图4给出了交叉验证结果散点图。
49.步骤6，将建立的夜间地表照度模型应用于经过去噪、辐射定标和几何精校正等处理的南京市区jl1-3b辐射亮度数据，得到南京市区夜间地表照度空间分布图，如图5所示。
50.步骤7，根据昼间高分辨率遥感图像运用公式（3）计算归一化植被指数ndvi，然后将图像红、绿、蓝进行真彩色合成，采用gram-schmidt融合方法与ndvi图像进行融合。参考南京市区路网矢量gis文件，通过人机交互目视解译从融合图像中提取各条道路的空间范围。考虑到街道两边行道树树冠对地面光辐射有遮蔽效应，通过人机交互目视解译从融合图像中提取行道树冠层空间范围，在提取出的道路空间范围中去除行道树冠层空间范围，得到道路无遮蔽空间范围，如图6所示。
51.步骤8，将南京市区道路无遮蔽空间范围与夜间地表照度空间分布图进行空间叠加分析，提取各条道路空间范围内的地表照度值。运用公式（4）、公式（5）、公式（6）计算各条道路空间范围内的地表照度平均值、标准差和变异系数，如表2所示。根据各道路的照度平均值、标准差和变异系数，结合《国家道路照明设计标准》和各道路等级对道路夜间照明质量进行监测和评价。
52.表2 南京主城区各道路夜间照明质量指标统计表本实例的实验结果说明，基于jl1-3b夜光遥感数据的道路夜间照明质量监测方法能够便捷高效地完成大范围道路夜间照明质量监测的任务，反映道路夜间照度分布情况，测量效率明显优于传统监测方法，并且不会对交通造成影响，能够为道路照明设计规划和道路安全性评估提供技术支撑。
53.以上仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，应视为本发明的保护范围。