之前给大家分享过一些关于税调周边的数据，例如税调周边的碳排放量、污染物等指标：
- 2007～2016 年税调周边 1～5km 范围内的碳排放总量：https://rstata.duanshu.com/#/course/0fade2f56be342ea81092c7db266b603
- 2007～2016 年税调周边的卫星反演年度污染物数据：https://rstata.duanshu.com/#/course/99edfc3dbf8b44bc864c6b679eeaa489
- …（更多数据可以在平台搜索）
也给大家分享了很多气象数据，例如平均气压、相对湿度等指标：
- 1980～2020年中国各省市区县分年、分月、逐日平均气压面板数据：https://rstata.duanshu.com/#/course/9289a403138141aeae22150b280cbe36
- 1980～2020年中国各省市区县分年、分月、逐日相对湿度面板数据：https://rstata.duanshu.com/#/course/9f63d4f663dc4d8084af8f4b68751eaa
- …（更多气象数据可以在平台搜索）
最近有小伙伴提出是不是可以使用类似的方法计算税调附近的气象数据，考虑到计算量比较大，因此我就直接帮忙处理好了。
数据概览
这里我提供的是税调周边 1～5km 范围内的年度气象面板数据，时间范围为 2007～2016 年。经过计算之后一共得到了 622.4 万条结果：
- 税调周边的年度气象数据.dta
包含的变量如下：
企业名称、法人代码、年份、sdgroup、sdid、平均地表气温_dist1km、平均地表气温_dist2km、平均地表气温_dist3km、平均地表气温_dist4km、平均地表气温_dist5km、平均气压_dist1km、平均气压_dist2km、平均气压_dist3km、平均气压_dist4km、平均气压_dist5km、平均气温_dist1km、平均气温_dist2km、平均气温_dist3km、平均气温_dist4km、平均气温_dist5km、平均风速_dist1km、平均风速_dist2km、平均风速_dist3km、平均风速_dist4km、平均风速_dist5km、日照时数_dist1km、日照时数_dist2km、日照时数_dist3km、日照时数_dist4km、日照时数_dist5km、相对湿度_dist1km、相对湿度_dist2km、相对湿度_dist3km、相对湿度_dist4km、相对湿度_dist5km、降水量_dist1km、降水量_dist2km、降水量_dist3km、降水量_dist4km、降水量_dist5km
例如，下图展示了天瀚科技(吴江)有限公司周边1~5km 范围内的年均气温：
处理方法
税调数据我使用的是之前分享的：
2007~2016 年税调企业地理信息数据（含经纬度及其所处的省市区县）：https://rstata.duanshu.com/#/course/76d38022cd004b09b2aa09647936beb0
其中经纬度是使用高德地图地理编码接口进行解析的。
具体处理方法如下：
- 首先分别生成税调企业地址周边 1km、2km、3km、4km、5km 的缓冲区（因为计算量非常大，所以这次没有计算 10km、15km 和 20km 的）。
- 读取各种气象指标的栅格数据。
- 使用缓冲区数据和读取到的栅格数据进行地理计算，就可以得到每个缓冲区内的各种气象指标的均值了。
- 按照这个思路循环各个年份的即可（不过运算量非常大，建议大家不要尝试，非常耗时）。
- 处理之后再整理即可得到分享给大家的
税调周边的年度气象数据.dta了。
附件中也提供了处理代码，感兴趣的小伙伴可以学习～
点击这里跳转到 RStata 短书平台获取附件：旧版本｜2007～2016 年税调周边 1～5km 范围内的年度气象面板数据
评论