注意：该夜间灯光亮度的单位是 0.1 nWatts/(cm^2 sr) 如果想转换成 nWatts/(cm^2 sr)，需要把数值乘以 0.1 。
今天给大家分享一份 2012年1月19日～2024年12月31日中国各省市区县日度夜间灯光亮度面板数据。
原始数据来源于 NASA 的 VNP46A2 - VIIRS/NPP Gap-Filled Lunar BRDF-Adjusted Nighttime Lights Daily L3 Global 500m Linear Lat Lon Grid 数据：https://ladsweb.modaps.eosdis.nasa.gov/archive/allData/5000/VNP46A2/
该数据把全球根据经纬度划按照 10˚x10˚ 的网格进行了划分，每个分块又包含了 2400x2400 个格点，分辨率大致为 500mx500m，数据均为 HDF5 格式。因此处理该数据的第一步就是要先找到中国所处的区块文件。
这些操作都可以在 R 语言中实现。
该数据的处理大致包含如下步骤：
- 通过地理计算的方法找到组成中国范围的区块（一共是 24 块）；
- 获取所有区块的下载链接，并从中筛选包含中国范围的区块进行下载；
- 读取每个区块文件并设置恰当的范围和坐标参考系；
- 合成栅格数据；
- 对栅格数据进行分省市区县汇总计算得到各省市区县每日的平均夜间灯光亮度面板数据。
前面四步的操作挺复杂，感兴趣的小伙伴可以学习最近的这个课程：
名师讲堂｜使用 R 语言下载和处理日度夜间灯光栅格数据：https://rstata.duanshu.com/#/course/270178a4d6bb4ee2ba8b66ae71de3b8f
至于计算各省市区县的平均夜间灯光亮度，感兴趣的小伙伴可以学习平台上的一些栅格转面板课程，例如：
使用 R 语言处理 Merra2 数据获取各省市区县的比湿、降水量、风速和气压数据：https://rstata.duanshu.com/#/brief/course/3239b87be1384488bdacb679cd467548
计算结果
经过计算一共得到了下面几个文件：
- 2012-01-19~2024-12-31各城市日度平均夜间灯光亮度面板数据.dta
- 2012-01-19~2024-12-31各区县日度平均夜间灯光亮度面板数据.dta
- 2012-01-19~2024-12-31各省份日度平均夜间灯光亮度面板数据.dta
数据预览如下：
绘图展示
面板数据
下图展示了省市区县的夜间灯光数据：
栅格数据
中间合成的日度夜间灯光亮度栅格数据如下所示：
连续性检查
再对连续性进行检查：
看起来还不错!
数据定制：日度夜间灯光亮度提取
由于处理得到的日度夜间灯光栅格数据非常大（数百GB），无法分享，因此特别推出了「日度夜间灯光亮度提取」的定制服务。
该服务仅限 RStata 会员。可以根据研究对象的经纬度提取其所处位置的夜间灯光亮度，也可以提取其周边一定范围的平均夜间灯光亮度。费用为 500 元/次。
数据引用格式
在论文中使用该数据请声明数据来源于RStata 数据中心，或者 RStata Data Center，并使用下述的格式引用：
RStata 数据中心: 2012年1月19日～2024年12月31日中国各省市区县日度夜间灯光亮度面板数据. 2025. https://tidyfriday.cn/rsdb2/
英文文献可以使用下面的格式引用：
RStata Data Center: Daily Nighttime Light Brightness Panel Data for Provinces, Cities, Districts, and Counties in China from January 19, 2012 to December 31, 2024. 2025. https://tidyfriday.cn/rsdb2/
另外也建议在说明数据来源的同时，根据本文的介绍在论文中描述数据的详细处理方法。
点击这里跳转到 RStata 短书平台获取附件：2012年1月19日～2024年12月31日中国各省市区县日度夜间灯光亮度面板数据
评论