夜间灯光数据dn值_DMSP/OLS和VIIRS/DNB夜间灯光影像的校正及拟合

原标题：DMSP/OLS和VIIRS/DNB夜间灯光影像的校正及拟合

夜间灯光数据能够表征人类活动强度和城市化进程，应用最广泛的夜间灯光数据为美国国防气象卫星(DMSP)搭载的可见红外成像线性扫描业务系统(OLS)数据及美国新一代国家极轨卫星(Suomi-NPP)搭载的可见光近红外成像辐射(VIIRS)传感器数据。夜间灯光影像将城镇夜光及其他发光体表征为亮值斑块即灯光区，使城镇等发光体明显区别于黑暗的无灯光背景区，避免了传统遥感影像城镇区域与非城镇区域混淆问题，为监测人类活动提供了数据源，在宏观尺度的城市研究中具有巨大的潜力和应用前景。

由于OLS传感器设计局限导致DMSP/OLS影像之间不连续、像元DN值过饱和等问题，数据于2013年截止，现被VIIRS白天/夜间波段(DNB)夜间灯光数据取代。为获得长时间序列稳定的夜间灯光数据，需要对两者夜间灯光数据进行校正与拟合。目前，部分学者针对DMSP/OLS数据提出了多种不变目标区域的相互校正方法，这些研究选择亮值年际变化小的区域作为标准区域，确立参考影像数据，对其他年份影像数据建立各类回归模型进行相互校正，形成了适用于我国区域完整的夜间灯光影像校正的技术方法。后又有学者提出了各种基于植被指数构建城市灯光指数为DMSP/OLS数据去饱和的方法，该类方法侧重单一影像校正，不同年际影像校正需获取相应年植被指数数据。而VIIRS/DNB数据在城市研究中需要去除火点噪声，且该数据存在时间短，不足以支持长时间序列的研究。大多数学者在长时间序列研究中对两种夜间灯光数据设置不同的处理方法，降低了结果的可比性，对两者数据拟合的研究鲜有 道。

为研究校正后获得的长时间序列夜间灯光数据的效果，对其模拟 会经济参量的能力进行探索，选择新疆连续年份的夜间灯光数据模拟 会经济参量。以往基于两种夜间灯光数据的 会经济参量研究在省级和县级两个尺度上与区域 会经济参量进行回归分析，大多数研究是单一年份数据的线性模拟或两种数据的对比模拟，且在部分模拟研究中没有对数据进行校正，少有研究长时间序列下校正后的模拟。利用校正后的夜间灯光数据对 会经济参量模拟的研究有助于提升数据可靠性与解释力，对数据广泛应用具有重要意义。

1 数据源及数据概况1.1 数据源

DMSP/OLS稳定夜间灯光影像数据集由美国国家海洋与大气管理局(NOAA)下属的美国国家地球物理数据中心(NGDC)发布，选择全部34期影像，影像数据集由6个不同的DMSP卫星F10(1992—1994年)、F12(1994—1999年)、F14(1997—2003年)、F15(2000—2007年)、F16(2004—2009年)、F18(2010—2013年)获取。全部影像均在NGDC的 站下载。

北半球的VIIRS/DNB合成的月平均灯光辐射数据集，冬季像元DN值受积雪影响，夏季像元DN值受植被影响，因此，选择2012—2017年10月影像数据，全部影像在NGDC的 站下载。

1.2 数据概况分析

DMSP/OLS稳定夜间灯光辐射数据产品包括由多个DMSP卫星传感器获取的1992—2013年共34期影像，影像去除了极光、野火等不稳定光源，以及月光、云的干扰，最后的数据值为无云图幅稳定光的年均灰度(DN)值，其范围是0~63。该影像数据集是由DMSP不同传感器获得的，而不同DMSP卫星搭载的OLS传感器具有不同的辐射探测性能、飞行前辐射标定，探测光谱波段也存在细微差别，且随着时间的推移，各传感器的辐射探测能力也逐渐衰退。影像数据未进行星上辐射校正和相互校正，造成同一个卫星传感器获取的连续不同年度的影像间相同位置的亮值像元DN值之间的异常。且由不同的传感器获取的同一年度的影像数据之间存在差异，主要表现为影像中的亮值像元的DN值总和不相等，以及影像间相同位置的亮值像元的DN值不同。长时间序列的DMSP/OLS稳定夜间灯光影像数据集存在的这些问题导致不同年份、不同传感器间的长时间序列数据不具有连续性和可比性。因此，利用该数据集进行长时间尺度的应用研究时，必须进行相互校正。但影像中像元最大值DN值为63，且许多城市核心区域出现了大量DN值为63的像元聚集情况，由此可推断影像中的像元DN值出现了饱和现象。像元DN值饱和会导致影像所表示的某一区域的灯光强度值低于其真实值，进而在灯光强度指标的相关研究中会产生误差，影响研究的准确性，因此影像需进行饱和校正。

VIIRS传感器的DNB探测通道继承了DMSP上OLS传感器的微光探测能力，增强了探测敏感度，提高了空间分辨率，两者传感器参数对比见表 1。VIIRS/DNB合成的月平均灯光辐射数据产品自2012年4月至今每月一期影像，影像过滤了杂散光、闪电、月光、云覆盖的影响，同时保留了极光、火、船及其他暂时性灯光。且该数据不存在DMSP/OLS稳定夜间灯光影像数据的问题，解决了像元DN值过饱和现象，其在轨检验校正程序也进一步提高了影像的清晰度。

表 1 OLS和VIIRS传感器的参数比较

传感器

OLS

VIIRS

通道

—

DNB

中心波长/μm

0.600

0.700

星下点分辨率/km

2.2×2.2

0.74×0.74

边缘分辨率/km

5.4×5.4

0.74×0.74

采样合并方向

扫描和轨道

无

可侦测信 最小值/

(W·m-2·sr-1)

3×10-5

4×10-5

应用ArcGIS软件对两者影像进行数据预处理，以我国县级行政边界面状数据为掩膜，提取所有数据的我国行政区域影像。所有数据的地理坐标系均为WGS-84坐标系，投影坐标系为兰伯特等面积投影坐标系，同时将影像数据空间重采样为1 km。应用ArcGIS软件数据管理模块的栅格属性工具集为影像创建栅格属性数据表，根据属性表提取影像DN值大于0的像元，对像元DN值进行统计分析，统计出影像的亮值像元DN值总和(total DN value，TDN)和亮值像元的总数量(total lit pixel，TLP)随时间变化的趋势，分析结果如图 1所示。其影像亮值像DN值总和(TDN)的计算公式为

图 2 影像校正流程

2.1 DMSP/OLS稳定夜间灯光影像数据的校正

在空间分辨率为1 km的夜间灯光数据中我国区域城市不存在衰退现象，前一年存在的亮值区域应在后一年继续存在；但分析数据发现前一年存在的亮值区域在后一年消失，且不同传感器同年影像的亮值区域也不一致，导致该问题的原因是存在不稳定像元与不一致像元。这些像元不能反映夜间灯光的空间变化，需剔除保留稳定亮值像元。在ArcGIS软件中提取DMSP/OLS 34期影像像元DN值大于0的区域，对不同传感器同年份影像进行相交分析，保留公共亮值像元；以2013年影像与下一年影像进行相交分析保留公共亮值像元为下一年稳定像元，以此类推，得到1992—2013年的稳定亮值像元影像，再用Extraction工具得到每个像元的DN值，按行列导入Excel中获得影像像元DN值矩阵。

2.1.1 传感器间相互校正

多传感器影像像元DN值异常波动是造成影像数据不连续的主要原因，为提高影像连续性，对影像进行传感器间的相互校正。参考文献[7-8]选取F16 2007年鹤岗市作为标定数据。选取一元二次方程式建立校正模型式

式中，DN、DNc分别为校正前、后像元灰度值；a、b、c为回归参数。

将每一期待校正影像的像元DN值矩阵与标定数据的像元DN值矩阵进行拟合，其回归参数见表 2。

表 2 DMSP/OLS数据集影像相互校正模型参数

传感器

年份

a

b

c

R2

F10

1992

-0.003 1

1.285 4

0.516 4

0.797 2

1993

-0.001 4

1.163 4

1.303 2

0.804 5

1994

0.002 2

0.953 1

0.857 7

0.819 8

F12

1994

0.002 1

0.978 1

1.992 9

0.798 2

1995

0.006 5

0.699 1

3.982 6

0.801 6

1996

0.009 3

0.562 7

4.208 5

0.815 6

1997

0.010 4

0.420 9

6.931 8

0.810 4

1998

0.011 8

0.476 4

6.289 5

0.813 9

1999

0.008 2

0.605 7

3.927 6

0.834 8

F14

1997

0.005 4

0.788 6

2.598 5

0.826 8

1998

0.006 1

0.754 3

2.981 2

0.798 6

1999

0.000 9

1.024 6

2.939 0

0.842 0

2000

0.005 4

0.786 1

9.033 7

0.854 1

2001

0.000 3

1.214 3

1.139 7

0.848 9

2002

0.000 0

1.237 7

0.793 0

0.810 3

2003

-0.003 2

1.415 6

0.954 5

0.880 2

F15

2000

0.007 5

0.645 2

3.707 0

0.812 5

2001

0.009 8

0.501 9

1.456 0

0.852 0

2002

0.000 0

0.990 1

2.453 9

0.837 9

2003

-0.003 9

1.391 7

1.816 4

0.835 7

2004

-0.008 4

1.806 1

0.904 8

0.829 8

2005

-0.005 1

1.532 8

0.946 2

0.836 8

2006

-0.004 2

1.607 2

1.035 9

0.803 9

2007

-0.008 9

1.830 7

0.628 9

0.796 1

F16

2004

-0.001 4

1.166 8

1.070 3

0.867 9

2005

-0.002 8

1.404 2

0.165 0

0.910 2

2006

-0.005 1

1.511 3

0.019 5

0.934 8

2007

0.000 0

1.000 0

0.000 0

1.000 0

2008

0.005 6

0.675 7

1.830 9

0.924 2

2009

0.007 1

0.585 9

2.959 7

0.925 1

F18

2010

0.008 5

0.439 6

3.620 7

0.896 4

2011

0.006 9

0.499 6

3.989 7

0.794 5

2012

0.008 5

0.398 7

4.059 8

0.822 3

2013

0.007 4

0.483 3

2.875 5

0.803 1

2.1.2 影像间连续性校正

传感器间相互校正后的影像数据削弱了像元DN值饱和程度，但影像不连续问题仍然存在，主要表现为多传感器同年影像相同位置像元DN值之间存在差异, 以及多传感器获取的多年影像的数据像元DN值异常波动，故需要进行影像间连续性校正。利用式(3)对同年不同传感器相互校正后的影像数据进行校正，形成1992—2013年每年唯一的数据集。

式中，DN(n-1, i)a、DN(n, i)a、DN(n+1, i)a分别表示数据集影像上第i个像元在n-1、n、n+1年度的DN值; n=1992, 1993, …, 2013。

2.2 VIIRS/DNB合成的月平均灯光辐射数据的校正

3 结果分析与讨论3.1 影像校正结果分析3.1.1 影像定性评价

图 4 校正后影像数据DN值统计分析

图 4显示校正整理后影像的TLP和TDN呈递增趋势，TLP在整个时间序列中缓慢增长，其2009—2017年的年间统计值差值逐渐减小；TDN在1992—2005年间平稳增长，2006—2013年间急剧增长，2014—2017年间增长速度又趋于平稳。这可能是由于城市发展到中后期阶段城区中心的灯光强度趋于稳定，灯光影像表现为稳定像元，城市通过向周边地区的扩张和合并继续发展，这一点通过观察VIIRS/DNB影像浮点型数据也可得出。

以上研究结果表明校正后的灯光影像解决了DMSP/OLS影像数据本身的问题和DMSP/OLS数据与VIIRS/DNB数据不可比的问题。VIIRS/DNB影像数据在数据处理时将空间分辨率重采样降低为1000 m，而用DMSP/OLS影像数据为参考重分类和回归校正后使VIIRS/DNB影像数据像元的DN值表达更像DMSP/OLS影像数据，使VIIRS/DNB影像数据损失了大量细节信息，丧失了数据本身在大气气溶胶、夜间云等应用研究中的优势[13-14]。但不可否认的是，校正后的数据更适合与DMSP/OLS影像数据组成长时间序列的夜间灯光数据，从而应用到城市化研究中。

城市化发展直接体现的是土地利用方式的改变，但夜间灯光影像亮值区域范围和数值的变化反映出城市土地利用的变化，且影像容易获取，校正的过程具有可重复性，校正后的影像连续稳定，像元DN值的饱和程度也有所降低，建立的灯光指数适合动态监测大尺度城市化强度与时空分析。

3.2 连续影像数据与 会经济参量的关系

图 5 夜间灯光影像拟合 会经济参量

4 结语

李雪萍, 贡璐. DMSP/OLS和VIIRS/DNB夜间灯光影像的校正及拟合[J]. 测绘通 ，2019(7)：138-146.

DOI:10.13474/j.cnki.11-2246.2019.0236.

相关资源：软件标书范本(技术部分)_软件技术标书-项目管理文档类资源-CSDN文库