新版本在这里：https://rstata.duanshu.com/#/brief/course/6de2b481c6a44bee8e98d991087b161d
之前给大家分享过 2012～2024 年的日度夜间灯光数据的爬取与方法：
名师讲堂｜使用 R 语言下载和处理日度夜间灯光栅格数据：https://rstata.duanshu.com/#/brief/course/270178a4d6bb4ee2ba8b66ae71de3b8f
也给大家分享过基于日度夜间灯光栅格数据处理的省市区县面板数据：
2012年1月19日～2024年5月20日中国各省市区县日度夜间灯光亮度面板数据：https://rstata.duanshu.com/#/brief/course/1854d495604b47d8b700c3009f23fae6
今天给大家分享一份基于这个数据计算的 2012～2024 年各省市区县政府周边的夜间灯光亮度及是否加班面板数据。
计算方法
该数据的计算方法参考了附件中的文献：「超时加班与劳动收入份额：基于卫星夜间灯光的经验证据」
论文中提供了两个判断政府是否加班的标准：
- 时间维度标准：该政府处的夜间灯光亮度大于该年节假日夜间灯光亮度的中位数；
- 空间维度标准：该政府处的夜间灯光亮度大于周边8个网格的中位数；
- 时空双维度：同时满足上面的时间和空间标准的日子才被视为加班日。
法定节假日数据之前也有给大家分享过：
2007~2024 年法定节假日及星期数据：https://rstata.duanshu.com/#/brief/course/63147c1f06924ddb89b152d1df84b78c
栅格数据的本质是一个个的格子：
高分辨率的栅格数据，这个夜间灯光数据的分辨率是 500mx500m的，很难看出来这种本质，但是低分辨率的栅格数据大家一看就明白了，例如这个栅格数据：
来源于课程「使用 R 语言绘制历年中国省市区县地图（小地图版本+长版）」：https://rstata.duanshu.com/#/course/379d19770956478ebc4d919910fca74c
每个政府都会落在其中的某个格子里面，例如 2024 年 12 月 31 日这天的安庆市政府：
数据概览
经过计算得到了如下数据：
- 2012～2024年各区县政府年加班天数
- 2012～2024年各区县政府周边的夜间灯光亮度及是否加班数据
- 2012～2024年各省政府年加班天数
- 2012～2024年各省政府周边的夜间灯光亮度及是否加班数据
- 2012～2024年各市政府年加班天数
- 2012～2024年各市政府周边的夜间灯光亮度及是否加班数据
注意：该夜间灯光亮度的单位是 0.1 nWatts/(cm^2 sr) 如果想转换成 nWatts/(cm^2 sr)，需要把数值乘以 0.1 。
其中，各省市区县政府周边的夜间灯光亮度及是否加班数据分别根据时间和空间维度标准判断了加班情况，变量名称的含义一目了然。
各省市区县政府加班天数数据就是根据各省市区县政府周边的夜间灯光亮度及是否加班数据统计的每个政府每年加班的天数，如果想计算比例，除以该年的总天数即可，需要注意的是 2012 年的数据是从 1 月 19 日开始的。
下图展示了各年区县政府的平均加班天数：
除此之外，大家还可以计算一些拓展指标，例如亮度差值：
点击这里跳转到 RStata 短书平台获取附件：旧版本｜2012～2024 年各省市区县政府周边的夜间灯光亮度及是否加班数据
评论