一、DMSP/OLS、NPP/VIIRS等夜间灯光数据之GDP空间化的理论介绍

一、前言

自上世纪90年代以来，DMSP/OLS夜间灯光数据问世，由于和其他遥感数据不同，该数据可以探测到地球夜间微弱的灯光。随着数据的积累，可以获取长时间序列的数据之后，随后常用于人类社会经济指标的研究中，例如：人口密度、GDP反演、能源消耗、二氧化碳等指标的研究。

以上社会指标的反演可以统称为XXXX空间化，可能小伙伴有疑问，为什么可以利用夜间灯光数据来做？其他的遥感影像不能用于空间化的研究，主要原因是在于其与社会指标相关性较高，可以进行拟合，进而进行反演，且该数据主要由DN值高低来体现人类社会活动的强弱，可以得到如下一些空间化分布结果：

那么如何我们才能利用夜间灯光数据实现各类社会指标的空间化的可视化并且对空间化的结果进行一些分析呢？

那我们要先理解相关的社会指标的空间化流程，可以仔细阅读以下的硕博论文。

基于夜间灯光数据广东省GDP空间化分析 - 中国知网 (cnki.net)https://kns.cnki.net/kcms/detail/detail.aspx?dbcode=CMFD&dbname=CMFD202101&filename=1020326578.nh&uniplatform=NZKPT&v=JXlCiC9owA7raZivlrkl41hqzHvhJ_DdXcEq1aKLbKL2-wNpQVHdZYyKl9a7ZoXp基于多源数据的人口空间化方法研究 - 中国知网 (cnki.net)https://kns.cnki.net/kcms/detail/detail.aspx?dbcode=CMFD&dbname=CMFD202201&filename=1021634491.nh&uniplatform=NZKPT&v=bwv1ZzB7kP7c1vuUaULhjoG43FX2v3hfVNxiCtLe_0NAPSzyOUcwLVTlBFFriCnB基于多源夜光遥感数据的1993-2018年京津冀城市群发展研究 - 中国知网 (cnki.net)https://kns.cnki.net/kcms/detail/detail.aspx?dbcode=CMFD&dbname=CMFD202102&filename=1021587801.nh&uniplatform=NZKPT&v=FSFTK35IaNEOJYHmgSK0McDTLJhsFg2fyeIMUpLCsOkrFvfFOlo_TSKHji2QZrDZ各位小伙伴可以先看一下上述的文章，在进行这方面的相关研究，主要学习理解一下方法技术这一章。

回归主题，本专栏主要从GDP空间化为例，那么给大家分析一下当前的基于DMSP/OLS、NPP/VIIRS等夜间灯光数据GDP空间化的相关理论介绍。

二、具体操作

（1）研究维度的选择

由于夜间灯光数据目前有三种：DMSP/OLS、NPP/VIIRS、珞珈一号夜间灯光数据，DMSP/OLS研究年份有1992-2013年，NPP/VIIRS研究年份有2012-至今，目前来说NPP/VIIRS还在持续更新中，珞珈一号目前最早的数据应该在2018年左右，珞珈一号优势明显就是分辨率高，可以达到130米，但是缺点也明显，数据覆盖不够，不管是年份还是区域。

 所以也就说目前要做基于夜间灯光数据长时间序列的研究比较多的做法是结合DMSP/OLS和NPP/VIIRS夜间灯光数据，也就是需要对两夜间灯光数据分别完成自身的校正之后，还要对校正完的数据再次进行连续性校正，得到长时间序列可用的夜间灯光数据集。

当然珞珈一号目前的研究热度也是很高，虽然目前比较难形成长时间序列的数据，目前来说不适用长时间序列的GDP空间化研究，但是其空间分辨率比起之前的夜间灯光数据都有很大的优势，达到130m，而且数据丰富（有日数据、月度数据），之前的仅仅只有月度数据，可以提高研究的现势性，未来也是应用前景很广。

（2）研究方法

目前基于夜间灯光数据GDP空间化研究都是分产业且分维度进行研究。

分产业其实很好理解，就是将GDP三个产业进行剥离出来，也就是说第一产业主要是以农、林、渔，你可以理解为这些产业并不需要什么灯光，将其剥离出来，研究表明分产业的一个方式可以进一步的提高拟合的相关性系数，也就是我们所说的R方。

分维度其实这个看别人的论文也好理解，主要有两种叫法时间维度和空间维度，一种是横向和纵向。其实时间维度就是对应纵向，空间维度就是对应横向，这一部分内容可以详细看硕博论文中如何介绍，看一下就能明白。

（3）分析角度

分析角度目前都是结合区域经济学的知识进行分析，包含经济联系强度、城市类型等等进行分析，最后提出相关的建议。

数据准备：对于GDP空间化，数据准备主要就是校正好夜间灯光数据、GDP统计数据、研究区域行政矢量数据。

第二就是从统计年鉴得到分产业的GDP统计数据，这里就要准备两种，一种全省的分产业数据，还有就是各个市的分产业数据，这是为分不同维度研究所准备的统计数据。

如果研究区域是全国的的话，需要准备全国的分产业数据，还有就是分省一级的分产业数据。通常来说做全国的GDP空间化会结合多源数据也就是结合多种数据按照权重实现GDP空间化。

步骤：其实这里主要是对空间化过程，对于初学者来说在研究过程中会碰到很多问题，例如，怎么反演，什么叫模型拟合等，甚至在反演过程中，出现很多不一样的问题。

对于空间化的流程基本都是一样：

准备好统计数据和灯光数据之后，提取夜间灯光指数，那么夜间灯光指数提取有两种方式，一种自动或者手动计算。将夜间灯光指数与GDP数据进行拟合，拟合完成之后得到方程，接下来的重点步骤就是：

第一步：建立渔网（包含一个点，一个面渔网）渔网像元大小，也就是实现像元提取的第一步；

第二步：将夜间灯光影像的DN值提取到第一步建立的渔网（点状渔网）；

第三步：将提取的DN值表导出到Excel；

第四步：将预先模拟获得的模型，用于计算想要的值（GDP、人口、能源等），然后导入到GIS中，连接到渔网中的面状渔网，然后进行渲染即可；

第五步：对反演结果进行校正。

有人会问模型哪里来，很简单，利用夜间灯光数据的平均灯光强度、综合灯光指数等与GDP、人口、能源利用SPSS或者excel直接拟合得到一个方程，这里要注意自变量和因变量的选择。

可能又有小伙伴要问，我怎么知道夜间灯光数据平均灯光强度、综合灯光指数这些东西怎么求？后续文章会一一介绍。