一、NPP/VIIRS夜间灯光数据校正一理论介绍

一、前言

我们在进行校正之前，先把NPP/VIIRS夜间灯光数据给大家介绍清楚，由于DMSP/OLS夜间灯光数据在2013年之后不再提供数据，美国为了能够持续获得该类型的数据，在2011年顺利的发射NPP卫星，其实NPP就是为替代DMSP卫星，其搭载了全新一代的传感器——VIIRS，相比OLS传感器具有高清晰度、高敏感性的优势，而且不存在饱和的情况，整体上来说具有较大优势。

但是该数据同样需要完成校正，由于该数据并未消除极光、火灾、渔火等，存在负值和极个别异常值和突变值的情况。也需要对该数据进行相关校正，当前主流的校正流程为阈值法、掩膜法，主要就是为了完成去负值、去异常值，完成NPP/VIIRS夜间灯光数据的连续性校正，得到可用连续性强的2012-2020年NPP/VIIRS夜间灯光数据。

数据下载链接：VIIRS Nighttime Light (mines.edu)https://eogdata.mines.edu/products/vnl/#monthly

校正流程包含：

（1）年度数据合成

NPP/VIIRS夜间灯光数据目前有日度、月度、年度数据，但是NOAA官网年度数据只有2015、2016年，为了得到长时间序列的可用的灯光数据，通常是下载月度数据合成年度数据。年度数据合成方法主要有平均值，说白了就是求平均值，将多个月份加在一起，除以月份数量。

（2）投影、裁剪、重采样

年度数据合成之后，将年度数据转换成Lambert等角圆锥投影坐标系，并重采样为1Km，方便后续和DMSP/OLS夜间灯光数据拟合校正。

（3）去负值、异常值

 也就是这些数量少且Value值比较大的就是我们说的异常值，其实为什么会这样？这个NPP/VIIRS夜间灯光数据不存在上限有关，导致其DN值波动范围较大，当然这和NPP/VIRRS夜间灯光数据传感器也有关系。

通常的做法就是对其进行设置阈值，也就是阈值法，将灯光DN值大于63设置为63，当然这里考虑将其设置为63的原因就是为了后续更好的和DMSP/OLS夜间灯光数据拟合，这只是一种做法。

其实也有一种做法就是将改年经济最发达的区域的最高DN值作为全年的DN值上限，其实这也很好的理解，全国通常经济经济最发达就是北上广深，例如北京市某一年DN值最高也就是259，那么全国应该不可能有太多地方的DN值高于259，也就是此时就是将259设置为最高DN值，当然也可以参考上海、广州、深圳市的DN值。

（4）连续性校正

NPP/VIIRS夜间灯光数据除了异常值的问题，DN值还存在突变问题，由于经济发展情况与夜间灯光数据关联性比较强，所以中国改革开放40年来，GDP是不断在增加，由此我们可以知道夜间灯光数据应该也是这样增加的，也就是说每一年的NPP/VIIRS夜间灯光数据DN值理论上来说保持前一年原有的一个水平要增加。也可以理解为假设2012年的像元亮度不会在2013年消失，且2013年的所有像元DN值不应该小于2012年，这里只是以这两年为例，其他年份就是以此类推。那么也就是需要对NPP/VIIRS夜间灯光数据进行连续性校正。