将 DMSP/OLS 夜间灯光时间序列扩展到 2013 年以后

1.综述

        国防气象卫星计划的业务线扫描系统 (DMSP-OLS) 传感器收集的数据已由美国国家海洋和大气管理局 (NOAA) 的地球观测组 (EOG) 存档和处理，以制作自 1994 年以来的全球夜间图像地图。多年来，EOG 开发了自动算法，通过消除火灾和渔船中的瞬态光，从 OLS 可见波段数据制作稳定光复合材料。短暂的灯光因其高亮度和短持续时间而被移除。但是，原来收集DMSP数据的6颗卫星逐渐从昼/夜轨道转移到了黎明/黄昏轨道，也就是更早的穿越时间。2014年初，F18卫星不再收集可用的夜间数据，并且重点已转移到处理来自可见红外成像辐射计套件 (VIIRS) 日/夜波段 (DNB) 数据的全球夜间图像。尽管如此，人们很快发现 F15 和 F16 卫星从 2012 年起开始收集黎明前的夜间数据。因此，从 2013 年开始，前几年建立的算法被扩展到处理 OLS 数据。此外，2013-DMSP 卫星 F18 和 F15 分别来自傍晚和黎明前的三个立交时间的夜间数据的存在，以及午夜之后的 VIIRS 的夜间数据的存在，使得对三个不同立交时间的图像进行相互校准成为可能并研究夜间灯光的昼夜模式。人们很快发现 F15 和 F16 卫星从 2012 年开始收集黎明前的夜间数据。因此，从 2013 年开始，前几年建立的算法被扩展到处理 OLS 数据。此外，2013-DMSP 卫星 F18 和 F15 分别来自傍晚和黎明前的三个立交时间的夜间数据的存在，以及午夜之后的 VIIRS 的夜间数据的存在，使得对三个不同立交时间的图像进行相互校准成为可能并研究夜间灯光的昼夜模式。人们很快发现 F15 和 F16 卫星从 2012 年开始收集黎明前的夜间数据。因此，从 2013 年开始，前几年建立的算法被扩展到处理 OLS 数据。此外，2013-DMSP 卫星 F18 和 F15 分别来自傍晚和黎明前的三个立交时间的夜间数据的存在，以及午夜之后的 VIIRS 的夜间数据的存在，使得对三个不同立交时间的图像进行相互校准成为可能并研究夜间灯光的昼夜模式。

2.材料和方法

2.1材料

       自 1970 年代中期以来，美国空军、国防部 (DoD) 一直在使用极地轨道国防气象卫星计划的操作线扫描系统 (DMSP-OLS) 传感器。自 1992 年以来，DMSP 数据已在 NOAA 的国家环境信息中心 (NCEI) 存档。夜间灯光处理由 NOAA/NCEI 的地球观测组于 1994 年建立，并转移到科罗拉多矿业学院佩恩公共政策研究所在 2020 年。

OLS 数据在全球范围内可用，空间分辨率为 2.7 公里，由两个光谱带组成，可见带 (0.4–1.1 µm) 和热带 (10.5–12.6 µm)。

        OLS 设计用于检测月光云，光电倍增管 (PMT) 用于增强可见波段信号。然而，由于这种强化，来自城市的灯光、火灾、渔船、煤气灯等也会被检测到。其中一些灯光的来源是短暂的，例如火灾和渔船，而城市灯光则更为持久。一年一度的“稳定”灯复合材料是通过检查进入复合产品的可见光带的分布，将瞬态源与更持久的源隔离开来的。虽然 OLS 在收集微光成像数据的能力上是独一无二的，但它也存在一些缺陷，例如空间分辨率粗糙、6 位量化以及有限的动态范围导致城市中心饱和。在以下部分中。

2.2创建 OLS 稳定灯组合的步骤

  2.21预处理     

     包含几个预处理步骤，以确保在稳定灯产品中只包含高质量的无云夜间数据。这是通过为 OLS 轨道创建伴旗图像来完成的。标志图像用于将 OLS 像素放入类别中。标志图像是 16 位的，它们是按位处理的。因此，通过打开特定位，每个 OLS 像素可以同时属于多个类。Stable Lights 处理中使用的标志类别是 DAYTIME、NIGHTTIME MARGINAL、ZERO LUNAR ILLUMINANCE、CLOUDS PRESENT 和 NO DATA。

太阳高度角有助于确定 DAYTIME 和 NIGHTTIME MARGINAL 标志。DAYTIME标志位为大于-6的太阳高度角设置，NIGHTTIME MARGINAL标志位为-15和-6之间的太阳高度角设置。NIGHTTIME MARGINAL 标志位涵盖了从夜间到白天的过渡区，与较暗的夜间数据相比，该终止区中的数据质量有所下降。

图示：亚洲上空轨道 F15201801162050 的夜间部分 ( a ) OLS 可见波段 (0.4–1.1 µm)；( b ) 伴随旗带，DAYTTIME 显示为红色，NIGHTTIME MARGINAL 显示为绿色。 

        零月球照度标志是使用月球照度值设置的，该值通过从美国海军天文台 [ 25 ]获取的算法计算得出。该算法根据月相、方位角和仰角来估计地球表面的月球照度，而月相、方位角和仰角又是根据每个 OLS 像素的纬度和经度以及每个扫描线的最低点像素处的时间计算得出的。当返回的月球照度小于 0.0005 lux 时，设置零月球照度标志。有时，OLS 夜间子轨道内可能存在没有实际数据但作为占位符存在的扫描线。为这些区域设置了 NO DATA 标志。

       轨道处理中的这一步骤包括与之前的方式相比最重要的变化之一。具有零月球照度标志位的子轨道进入稳定灯产品的处理。然而，这些亚轨道可能有极光、太阳眩光和突然的增益变化。通过 F15 在 OLS 轨道上看到的突然增益变化是由阈值引起的，阈值逐渐改变以跟上 OLS 变化的最大增益。当轨道从夜间转向白天时，OLS 的最大增益被移动以防止饱和。尽管每周调整增益阈值，但它们并未记录为数据流的一部分。以前，分析师必须通过每个子轨道并目视筛选增益的突变（通过 F15 卫星）以及极光和太阳眩光的存在，选择要包含在 Stable Lights 产品中的数据的开始和结束扫描线。一名分析师花费了几个小时来遍历一个月的所有子轨道并完成这项任务。然而，在处理扩展系列时，神经网络 (NN) 被用于分割轨道并找到极限并移除被极光和太阳眩光照亮的 DMSP 轨道部分。