1998～2014 年工企周边的污染物浓度年度面板数据

1017 字 约 3 分钟 数据资料 面板数据空气质量上市公司

之前给大家分享过很多卫星反演污染物数据，例如 PM2.5、SO2等指标：

也给大家分享了不少关于工企周边的数据，例如工企周边的碳排放量、夜间灯光等指标：

最近有小伙伴提出是不是可以使用类似的方法计算工企附近的卫星反演污染物数据，考虑到计算量比较大，因此我就直接帮忙处理好了。

数据概览

为了方便大家使用，我计算了工企周边 1～5 km 范围各种污染物指标的年度数据，时间范围为 1998～2014 年，经过计算之后一共得到了 471.8 万条结果：

这里我提供的数据格式为 dta 格式，数据包含的变量如下：

group、年份、gqid、企业匹配唯一标识码、组织机构代码、企业名称、县代码、县、省代码、省、市代码、市、经度、纬度、PM1_1km、PM1_2km、PM1_3km、PM1_4km、PM1_5km、PM10_1km、PM10_2km、PM10_3km、PM10_4km、PM10_5km、PM2_5_1km、PM2_5_2km、PM2_5_3km、PM2_5_4km、PM2_5_5km、SO2_1km、SO2_2km、SO2_3km、SO2_4km、SO2_5km、SO4_1km、SO4_2km、SO4_3km、SO4_4km、SO4_5km、二甲硫醚_1km、二甲硫醚_2km、二甲硫醚_3km、二甲硫醚_4km、二甲硫醚_5km、粉尘_1km、粉尘_2km、粉尘_3km、粉尘_4km、粉尘_5km、黑碳_1km、黑碳_2km、黑碳_3km、黑碳_4km、黑碳_5km

绘图展示

为了更直观地感受这份数据，我以常州电梯厂有限公司为例，绘制了该公司周边 1~5km 范围内年均 PM10 和 SO2 浓度。PM10 浓度：

SO2 浓度：

处理方法

处理方法类似之前的课程「使用 R 语言计算工企周边一定范围内的平均 PM2.5 浓度」，感兴趣的小伙伴可以学习：

使用 R 语言计算工企周边一定范围内的平均 PM2.5 浓度：https://rstata.duanshu.com/#/course/016330182fb64ecd9fab8ada0d816130

工企的经纬度数据可以使用之前分享的：

1998～2014 年工业企业数据库地理位置数据（含经纬度、所处省市区县、南北方属性、东西部属性及相关距离）：https://rstata.duanshu.com/#/course/1c0d065e272c46b49cc528a445c2d3d9

具体处理方法如下：

  1. 首先分别生成工企地址周边 1km、2km、3km、4km、5km 的缓冲区（因为计算量非常大，所以这次没有计算 10km、15km 和 20km 的）。

  2. 读取各种污染物指标的栅格数据。

  3. 使用上面的缓冲区数据和读取到的栅格数据进行地理计算，就可以得到每个缓冲区内的各种污染物指标的均值了。为了便于大家理解这个计算过程，我绘制了如下的图，其中方块表示的是污染物栅格数据，圆圈表示的是工企周边 5km 范围。计算过程实际上就是提取圆圈范围内的区域求平均值。圆圈边缘的按照比例计算进入圈内的部分。

  4. 按照这个思路循环各个年份的即可（不过运算量非常大，建议大家不要尝试，非常耗时）。

  5. 处理之后再整理即可得到分享给大家的工企周边的卫星反演年度污染物数据.dta了。

附件中也提供了匹配使用的代码供大家参考。

点击这里跳转到 RStata 短书平台获取附件：1998～2014 年工企周边的污染物浓度年度面板数据

相关文章

评论