新版本在这里：https://rstata.duanshu.com/#/course/49a7c949a7b04ef2b2006b46841f9139
注意：该夜间灯光亮度的单位是 0.1 nWatts/(cm^2 sr) 如果想转换成 nWatts/(cm^2 sr)，需要把数值乘以 0.1 。
今天给大家分享一份 2012年1月19日～2024年5月20日中国各省市区县日度夜间灯光亮度面板数据。
原始数据来源于 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 语言处理 Merra2 数据获取各省市区县的比湿、降水量、风速和气压数据：https://rstata.duanshu.com/#/brief/course/3239b87be1384488bdacb679cd467548
计算结果
经过计算一共得到了下面几个文件：
- 省级：2012-01-19~2024-05-20各省份日度平均夜间灯光亮度面板数据.dta
- 市级：2012-01-19~2024-05-20各城市日度平均夜间灯光亮度面板数据.dta
- 县级：2012-01-19~2024-05-20各区县日度平均夜间灯光亮度面板数据.dta
绘图展示
下图展示了省市区县的夜间灯光数据：
中间合成的日度夜间灯光亮度栅格数据如下所示：
再对连续性进行检查：
看起来还不错!
数据定制：日度夜间灯光亮度提取
由于处理得到的日度夜间灯光栅格数据非常大（数百GB），无法分享，因此特别推出了「日度夜间灯光亮度提取」的定制服务。
该服务仅限 RStata 会员。可以根据研究对象的经纬度提取其所处位置的夜间灯光亮度，也可以提取其周边一定范围的平均夜间灯光亮度。费用为 500 元/次。工企、上市公司、税调的就不用了，后面很快就会分享～
点击这里跳转到 RStata 短书平台获取附件：旧版本｜2012年1月19日～2024年5月20日中国各省市区县日度夜间灯光亮度面板数据
评论