夜间灯光数据dn值_中国区域DMSP/OLS夜间灯光影像的校正

DMSP/OLS夜间灯光影像由于传感器未经过星上定标，导致长时间序列的各期影像之间缺乏连续性和可比性，且由于光谱分辨率的限制，使得大量的城市中心像元出现DN值为63的集聚饱和现象等问题，因此学者们对以上存在的问题提出相应的解决办法。2009年，文献[12]提出了不变目标区域的方法对全球的夜间灯光影像的不连续和不可比问题进行了校正，该方法以Sicily为不变动的参考区域，以F121999为参考数据集，将待校正影像与F121999参考数据集进行比较建立回归方程。该方法较适用于全球夜间灯光数据的校正，对于某一国家尺度的适用性不强。文献[3]在Elvidge的基础上，提出了适合中国区域的夜间灯光影像的校正方法。首先统计分析了中国主要城市1992—2008年的GDP统计数据和建成区面积数据，确定黑龙江省的鸡西市为校正中国区域的不变动参考区域，选择F16的2007年数据作为参考数据集，实现长时间序列DMSP/OLS夜间灯光数据的相互校正、年内融合和年际间校正，却忽略了影像的饱和问题。文献[13]基于地理探测器模型验证分析，提出利用单元路网长度去饱和算法。结合2006年的辐射定标产品数据对2013年DMSP/OLS稳定灯光数据进行饱和校正。文献[14]结合DMSP/OLS数据和Landsat数据，建立去除水体并且缓解灯光强度饱和的指数RwNTLI，结果表明该方法在一定程度上降低了影像的饱和问题。以上几位学者提出的方法能够在一定程度上缓解影像的饱和问题，但是仅局限在单期影像的校正，无法在长时间序列上进行研究。因此，本文旨在对1992—2013年长时间序列的DMSP/OLS夜间灯光影像进行饱和校正、相互校正、年内融合及年际间校正等一系列校正内容，从而解决夜间灯光影像的饱和、影像间不连续及不可比等的问题。

1 校正方法

具体校正思路为：①确定黑龙江省的鸡西市为参考区域；②选取已经经过饱和校正的辐射定标产品数据中的2006年F16传感器的数据作为对连续稳定灯光影像进行饱和校正、相互校正的参考数据集；③对1992年至2013年34期夜间灯光影像进行相同年份的年内融合及不同年份的年际间校正，从而解决灯光影像的饱和、影像间不连续及不可比等问题。

1.1 影像预处理

原始稳定灯光影像的坐标为WGS-84，空间分辨率为30″。研究区域为中国区域，因此根据原始影像裁剪出中国区域影像。为了使投影面积变形最小，对坐标系进行变换，椭球体为D_Krassovsky_1940，投影坐标为Krassovsky_1940_Albers。统一影像分辨率为1000 m。

1.2 中国区域的DMSP/OLS夜间灯光数据校正1.2.1 影像数据分析

图 1所示为未经校正的中国区域的所有亮值像元(像元DN值大于0)的数量，图 2所示为未经校正的中国区域的所有亮值像元的DN值总和(TDV)。可以看出：①同一年份的数据可由不同