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