今天给大家分享一份各省市区县乡镇 1986～2024 年夜间灯光面板数据。
原始数据是全球范围的栅格数据。为了方便大家使用，我提取了中国范围的栅格数据并汇总成了省市区县、乡镇面板。
引用该数据可以使用：
Tian, Y., Cheng, K. M., Zhang, Z., Zhang, T., Li, S., Yan, D., & Xu, B. (2025). A novel modeling framework and data product for extended VIIRS-like artificial nighttime light image reconstruction (1986-2024). arXiv preprint arXiv: 2508.00590.
数据概览
经过栅格数据到面板数据的转换，得到如下文件：
- 1986～2024年中国各城市夜间灯光面板数据.dta
- 1986～2024年中国各区县夜间灯光面板数据.dta
- 1986～2024年中国各省市夜间灯光面板数据.dta
- 1986～2024年中国各乡镇夜间灯光面板数据.dta
每份数据都包含相应的区域名称和代码以及 mean、min、max、sum、sd 指标。这些统计量可以这么理解：栅格数据是类似一个个格点的数据，每个格点都包含一个值，可以把区域内的这些值求平均得到 mean 指标，求总和得到 sum 指标…… 因此如果想要使用区域平均夜间灯光亮度应该使用 mean 指标。
各文件预览如下：
为了便于大家直观的判断这份数据的可靠性，我绘制了下面这些图表：
栅格数据
首先从原始数据中提取中国范围的栅格数据，绘制 2024 年各地平均夜光强度图示：
面板数据
下图展示了 2024 年各城市平均夜光强度：
连续性检查
我们再来检查下数据的连续性，历年各省市夜光强度变化趋势如下：
连续性还不错！
处理方法
栅格数据转面板的方法都差不多，感兴趣的小伙伴可以学习下面两个课程：
- 使用 R 语言处理 Merra2 数据获取各省市区县的比湿、降水量、风速和气压数据：https://rstata.duanshu.com/#/brief/course/3239b87be1384488bdacb679cd467548
- 中国各省市碳排放量是如何计算的？R 语言栅格数据转面板数据：https://rstata.duanshu.com/#/brief/course/75de598dcfcc4ad0b9ff28bff27f6b83
代码
为方便大家学习，附件中还提供了数据处理和绘图代码以供参考：
数据引用格式
由于该数据包含较多 RStata 处理的内容，在研究中使用该数据请使用清晰的方式注明数据来源于 RStata 或者 RStata 数据中心，并使用如下格式引用：
RStata 数据中心: 1986～2024 年中国各省市区县、乡镇夜间灯光面板数据（田一禾版本）. 2026. https://tidyfriday.cn/rsdb2/
英文文献可以使用下面的格式引用：
RStata Data Center: Panel Data on Nighttime Lights by Province, City, District, County, and Township in China, 1986–2024 (Version by Yihe Tian). 2026. https://tidyfriday.cn/rsdb2/
点击这里跳转到 RStata 短书平台获取附件：1986～2024 年中国各省市区县、乡镇夜间灯光面板数据（田一禾版本）
评论