今天给大家分享一份2012～2024年各乡镇年度、月度、日度平均夜间灯光亮度。
之前给大家分享过各省市区县的日度夜间灯光面板数据：
2012年1月19日～2024年5月20日中国各省市区县日度夜间灯光亮度面板数据：https://rstata.duanshu.com/#/brief/course/1854d495604b47d8b700c3009f23fae6
由于乡镇的计算量会大非常多，所以当时并没有计算。最近有小伙伴表示还需要乡镇的，所以花了大概一周的时间又把乡镇的计算了下。
注意：该夜间灯光亮度的单位是 0.1 nWatts/(cm^2 sr) 如果想转换成 nWatts/(cm^2 sr)，需要把数值乘以 0.1 。
日度夜间灯光原始数据来源于 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 语言下载和处理日度夜间灯光栅格数据：https://rstata.duanshu.com/#/brief/course/270178a4d6bb4ee2ba8b66ae71de3b8f
中有讲解过。
按照上面的课程中讲解的方法可以处理得到日度夜间灯光的栅格数据。栅格数据分乡镇平均就可以得到乡镇的面板数据了，该方法可以学习这个课程：
使用 R 语言处理 Merra2 数据获取各省市区县的比湿、降水量、风速和气压数据：https://rstata.duanshu.com/#/brief/course/3239b87be1384488bdacb679cd467548
这里也不再赘述了。
年度和月度的数据是根据日度数据分组平均得到，为了方便大家使用，我也放在了附件中。附件中也提供了计算使用的乡镇矢量数据，该矢量数据来源于各省的天地图网站，大致是 2021 年的，不过部分省份不完全是。
附件中包含了如下三个文件：
- 2012～2024年各乡镇平均月度夜间灯光亮度.csv
- 2012～2024年各乡镇平均日度夜间灯光亮度.csv
- 2012～2024年各乡镇平均年度夜间灯光亮度.csv
注意 csv 格式的文件直接使用 Excel 打开会乱码，建议使用 Stata 或者 R 语言读取使用。
下图展示了 2024年1月各乡镇平均夜间灯光亮度：
2024年2月11日各乡镇平均夜间灯光亮度：
这两幅图均是使用 Stata 绘制，感兴趣的小伙伴可以学习如下课程：
使用 Stata 绘制中国乡镇地图及根据经纬度判断所处的乡镇：https://rstata.duanshu.com/#/brief/course/aaf11610c38a488c84b8e1d574561950
使用 R 语言的小伙伴可以学习这个课程：
使用 R 语言绘制中国乡镇地图及根据经纬度判断所处的乡镇：https://rstata.duanshu.com/#/brief/course/eea9db6b6cd3420e85ab4d8dbc9b4eda
点击这里跳转到 RStata 短书平台获取附件：2012～2024 年各乡镇年度、月度、日度平均夜间灯光亮度
评论