今天给大家分享一份 1992～2024 年中国各省市区县乡镇类 DMSP-OLS 夜间灯光面板数据。
原始数据是中国区域的栅格数据，来源于 Harvard Dataverse（https://dataverse.harvard.edu/dataset.xhtml?persistentId=doi:10.7910/DVN/GIYGJU ），为了方便大家使用，我使用 R 语言把栅格数据转换成了各省、市、区县、乡镇的面板数据。
引用该数据可以使用：
Wu, Yizhen; Shi, Kaifang; Chen, Zuoqi ; Liu, Shirao; Chang, Zhijian, 2021, “An improved time-series DMSP-OLS-like data (1992-2024) in China by integrating DMSP-OLS and SNPP-VIIRS”, https://doi.org/10.7910/DVN/GIYGJU, Harvard Dataverse, V7
数据概览
经过栅格数据到面板数据的转换，得到如下文件：
- 1992～2024年中国各城市夜间灯光面板数据.dta
- 1992～2024年中国各区县夜间灯光面板数据.dta
- 1992～2024年中国各省市夜间灯光面板数据.dta
- 1992～2024年中国各乡镇夜间灯光面板数据.dta
每份数据都包含相应的区域名称和代码以及 mean、min、max、sum、sd 指标。这些统计量可以这么理解：栅格数据是类似一个个格点的数据，每个格点都包含一个值，可以把区域内的这些值求平均得到 mean 指标，求总和得到 sum 指标…… 因此如果想要使用区域平均夜间灯光亮度应该使用 mean 指标。
例如城市面板是这样的：
为了便于大家直观的判断这份数据的可靠性，我绘制了下面这些图表：
栅格数据
原始数据即为中国范围的栅格数据。2024 年各地平均夜光强度绘图效果如下：
面板数据
下图展示了 2024 年各城市平均夜光强度：
连续性检查
由于这份数据是根据 DMSP-OLS 和 SNPP-VIIRS 整合得到的，所以我们再检查下数据的连续性，1992～2024 年中国各省市夜光强度变化趋势如下（可以看到该夜间灯光数据是以 DN 值表示，无单位）：
处理方法
栅格数据转面板的方法都差不多，感兴趣的小伙伴可以学习下面两个课程：
使用 R 语言处理 Merra2 数据获取各省市区县的比湿、降水量、风速和气压数据：https://rstata.duanshu.com/#/brief/course/3239b87be1384488bdacb679cd467548
中国各省市碳排放量是如何计算的？R 语言栅格数据转面板数据：https://rstata.duanshu.com/#/brief/course/75de598dcfcc4ad0b9ff28bff27f6b83
点击这里跳转到 RStata 短书平台获取附件：1992～2024 年中国各省市区县乡镇类 DMSP-OLS 夜间灯光面板数据（施开放版本）
评论